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ABSTRACT 
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URLIM , a unified  radome  limitations  computer 
program  has  been  developed  to  aid  the  radome  design 
engineer  by  providing  a definition  of  the  maximum  flight 
performance  capabilities  of  radome  materials.  URLIM 
numerically  determines  the  response  of  the  radome  to 
aerodynamic  heating  and  loading.  It  computes  the  fol- 
lowing as  functions  of  trajectory  time:  thermal  stress; 
radar  boresight  error  slopes;  missile-radome  attachment 
stresses  caused  by  maneuvers,  pressure,  and  drag  forces; 
and  the  onset  of  melting.  The  basic  output  of  the 
program  is  a notation  of  trajectory  time  at  which  the 
radome  reaches  its  design  limitations.  Many  options  are 
available  to  the  user  of  the  URL in  program  that  provide 
a wide  variety  of  analysis  capability.  For  this  reason 
URLIM  may  also  be  considered  as  a general  purpose  aero- 
dynamic heat  transfer  program  as  well  as  a specific 
purpose  radome  limitations  program.  Volume  1 of  this 
report  presents  the  theoretical  background  of  the 
analysis  technigues  used  in  URLIM ; Volume  2 provides  a 
detailed  explanation  cf  how  to  use  URLIM. 
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PREFACE 


The  purpose  of  this  volume  is  to  provide  a detailed 
explanation  of  how  to  use  the  URLIM  program.  Also 
included  herein  is  a discussion  of  the  general  objec- 
tives of  the  UFLIH  program  itself  (i.e.r  what  physical 
conditions  it  can  simulate  and  what  its  essential  out- 
puts are) . This  later  discussion  provides  an  overview 
of  the  information  required  by  the  program  so  that  the 
specific  data  descriptions  can  be  understood  in  their 
proper  context.  Also  embodied  in  this  volume  (specifi- 
cally in  the  Appendixes)  is  a description  of  the  com- 
plete subroutine  library  that  complements  the  URLIM 
program.  These  routines  comprise  the  Standard  Heat 
Transfer  Program  (SHTP)  library  and  are  documented 
herein  for  reference. 
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1.  BACKGROUND  AND  INTRODUCTION 


The  Unified  Radome  Limitations  Computer  Program 
(URLIM)  has  been  developed  to  improve  significantly  the 
accuracy  and  speed  with  which  a radome  designer  may 
define  the  permissible  operating  regime  of  a supersonic 
missile.  Moreover,  the  program  may  have  greater  utility 
as  a preliminary  systems  design  tool  by  affording  the 
ability  to  rapidly  assess  a candidate  radome  material's 
usefulness  for  high  speed  flight.  A secondary  aspect  of 
the  program,  but  one  that  is  no  less  useful,  is  that  the 
structure  being  analyzed  need  not  be  a radome  but  may  be 
any  body  subjected  to  transient  heating.  The  program 
contains  numerous  options  that  are  user-selectable  and 
that  provide  a means  of  controlling  the  amount  and  type 
of  analysis  performed. 

Before  a detailed  discussion  of  the  use  of  the 
URLIM  program  is  presented,  a general  discussion  of 
radome  limitations  will  be  given.  The  discussion  will 
include  a definition  of  the  data  requirements  for  URLIM. 

The  flights  of  current  and  proposed  tactical  mis- 
siles through  the  atmosphere  occur  at  relatively  high 
speeds  so  that  the  target  seeker  equipment  requires  a 
radome  to  protect  it  from  aerodynamic  heating  and  wind 
leading.  Also,  the  radome  provides  aerodynamic  fairing 
to  reduce  drag.  Consequently,  rademes  used  on  missiles 
must  satisfy  a number  of  criteria  to  be  successful. 
First,  they  must  be  dielectric  materials,  (i.e.  trans- 
parent to  radar  frequency  electromagnetic  (EM)  radia- 
tion) . Second,  they  must  be  streamlined,  so  as  to  offer 
low  aerodynamic  drag.  Third,  they  must  be  able  to  with- 
stand the  environment  created  by  the  high  speed  flight 
of  the  missile  through  the  atmosphere. 

In  order  to  satisfy  these  three  basic  design 
criteria  the  radome  designer:  (1)  selects  a plastic  or 
ceramic  material  of  suitable  thickness,  (2)  shapes  it 
aerodynamically , and  (3)  attempts  to  define  how  fast  and 
how  far  the  material  will  fly  before  it  no  longer  satis- 
fies the  first  two  criteria.  This  latter  aspect  of  a 
radome  design  is  seen  to  be  a definition  of  the  radome' s 
flight  limitations  and  is  a complicated  analysis  task 
for  the  following  reasons:  (1)  the  aerodynamic  heating 
and  loading  environment  must  be  defined  reliably. 
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(2)  the  effects  on  radar  transmission  must  be  deter- 
mined, and  (3)  the  response  of  the  material  to  transient 
heating  oust  be  calculated.  To  be  more  specific,  a 
radoae  limitations  analysis  must  examine  the  thermal  and 
mechanical  response  of  the  radome  over  all  possible 
trajectories  commensurate  with  the  missile's  propulsion 
capabilities.  A flight  limit  would  then  be  presented  as 
plots  of  all  of  these  trajectories  with  the  failure 
points  marked  on  each  trajectory.  Since  there  are  many 
combinations  of  trajectory  parameters,  it  must  be 
decided  how  to  "plot"  a trajectory  and  how  to  select  a 
reasonable  number  of  parameters. 

Trajectories  have  three  basic  aspects:  a velocity 
history,  an  altitude  history,  and  a downrange  history 
(we  assume  no  cross-range  components  in  the  trajectory 
with  no  loss  of  generality).  The  flight  regime  of  a 
missile  may  generally  be  depicted  as  an  area  in  an 
altitude-versus-range  plot  (an  elevation  view  of  the 
flight  space).  Such  a picture  is  sketched  in  Fig.  la. 
The  figure  shows  several  flight  paths  for  various  launch 
angles  ^nd  ceiling  altitudes.  For  each  one  of  these 
flight  paths,  a variety  of  velocity  histories  could  be 
used  depending  on  the  propulsion  system  used  by  the 
missile  and  other  considerations.  Figure  1b  shows  some 
typical  velocity  history  possibilities.  To  complete  the 
description  Fig.  1c  shows  the  altitude  histories  that 
would  result  from  the  trajectories  sketched.  To  des- 
cribe the  flight  potential  of  a missile,  it  would  suf- 
fice to  evaluate  each  of  the  velocity  histories  over 
each  of  the  flight  paths.  In  order  to  show  the  flight 
limitations  of  the  missile  radome,  a dot  is  placed  on 
each  velocity  history  and  on  the  appropriate  position  in 
the  altitude-range  plot.  The  velocity  history  plot 
(Fig.  1b)  has  five  such  dots  plotted  and  the  dashed  line 
plotted  through  the  dots  approximates  the  locus  of  all 
failure  points  on  all  similar  velocity  histories.  There 
are  also  five  dots  plotted  on  the  altitude-range  figure 
that  depict  the  flight  limitations.  In  this  way  dif- 
ferent flight  paths  are  flown  with  several  velocity 
histories,  and  in  each  case  a velocity  limit  plot  is 
made.  Corresponding  points  on  the  altitude-range  plot 
can  also  be  drawn.  Combining  these  data  will  then  give 
a "plot"  of  the  particular  radome 's  limitations  in  the 
atmosphere.  For  any  particular  radome  material,  as  dif- 
ferent wall  thicknesses,  or  different  launch  angles,  or 
different  flight  path  shapes  (e.g.  line-of-sight  versus 
boost-cruise)  are  used  the  resulting  velocity  limits  may 
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a.  Hypothetical  mission 
profiles 
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be  significantly  different.  While  it  is  helpful  to  know 
the  permissible  flight  regime  in  as  many  aspects  as 
possible,  the  most  useful  information  tor  radoae  design 
is  the  velocity  limit  on  a particular  flight  path. 
Because  of  this,  the  velocity  limit  curve  (c.f.  Fig.  lb) 
will  receive  the  emphasis  throughout  this  report.  With 
ci  set  of  velocity  limit  curves  for  several  flight  paths 
and  the  downrange  history  information,  the  other  ways  ct 
displaying  the  flight  limits  can  be  generated  easily. 

The  descriptions  and  discussions  that  follow  in 
this  report  will  be  concerned  with  the  following  four 
radome  failure  modes:  (1)  temperatures  in  the  radome 
that  exceed  the  melt  temperature  of  the  radome  material, 
(2)  transient  temperature  gradients  that  cause  thermal 
stresses  sufficient  tc  fracture  the  radome,  (3)  tem- 
perature changes  in  the  radome  wall  that  cause  radar 
aberations  (boresight  errors) , and  (4)  aerodynamic  pres- 
sure and  maneuver  loads  that  cause  stresses  at  the 
radome-missile  attachment  point  sufficient  to  cause 
f racture. 

Another  major  cause  of  radome  failure  is  rain 
erosion  but  this  material  response  is  not  covered  by  the 
U BUM  program.  However,  the  problem  of  rain  erosion  is 
serious  and  difficult  to  predict  and  must  be  considered 
in  any  radome  design.  A radome* s response  to  rain 
erosion  is  generally  characterized  by:  (1)  the  severity 
of  the  rain  field  (droplet  size  and  number/unit  volume), 
(2)  the  missile's  velocity,  and  (3)  the  shape  of  the 
radome.  Moreover,  the  susceptibility  of  any  given 
radome  to  rain  erosion  will  depend  primarily  on  the 
missile's  velocity,  with  only  secondary  dependence  on 
transient  material  responses.  In  other  words,  a rain 
erosion  limit  will  state  a combination  of  rain  field 
intensity,  length  of  rain  field,  and  velocity  beyond 
which  the  missile  cannot  fly. 

The  four  limits  of  melting,  thermal  stress,  aero- 
dynamic load,  and  boresight  error  are  secondary  effects 
that  result  from  the  transient  thermal  response  of  the 
radome  and  as  such  are  very  trajectory-dependent.  For 
example,  an  alumina  radome  flown  on  a trajectory  with 
low  acceleration  could  eventually  fly  at  a very  high 
speed  so  long  as  the  transient  thermal  gradients  were 
low  and  severe  thermal  stresses  were  not  developed. 
Conversely,  an  alumina  radome  flown  at  a high  accelera- 
tion will  fail  due  to  thermal  shock  quite  early  in  the 
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flight.  In  summary,  the  secondary  material  responses 
are  highly  dependent  on  trajectory  parameters  and  must 
be  evaluated  via  trajectory  simulation  to  determine 
their  attendant  flight  limits. 

The  example  cited  above  of  the  alumina  radome  with 
different  flight  limits  on  different  trajectories  points 
out  an  interesting  source  of  confusion  with  regard  to 
thermal  stresses.  An  appropriately  selected  set  of 
trajectories  can  produce  one  set  of  flight  limits  while 
a different  set  of  trajectories  over  the  same  range  of 
velocities  and  altitudes  can  produce  another  set  of 
limits.  This  ambiguity  is  possible  because  thermal 
stresses  arise  from  temperature  gradients  and  simply 
raising  a radome* s temperature  uniformly  (or  with 
moderate  difference)  will  not  produce  appreciable 
stresses.  Therefore,  it  is  not  inconceivable  that  a set 
of  trajectories  could  be  tailored  for  any  given  material 
wherein  thermal  stresses  would  never  be  a limit.  It  is 
also  possible  that,  for  some  materials,  these  trajec- 
tories might  not  be  too  useful  for  a tactical  missile. 

It  will  suffice  here  only  to  mention  these  "stress-free” 
trajectories  and  note  that  the  DRUM  pogram  will  accept 
any  velocity  history  as  input.  A series  of  analyses 
could  be  run  with  URL IM  to  discover  the  stress  failure 
free  paths,  if  desired.  From  this  discussion  we  concl- 
ude that,  when  thermal  stress  is  the  overall  limiting 
factor,  the  definition  of  a flight  limit  for  a radome 
design  will  not  be  unique.  Moreover,  the  region  of 
uncertainty  should  be  explored  by  the  URLIM  user. 

In  order  to  use  URLIM  to  perform  a limit  analysis, 
the  user  must  establish  a model  of  the  radome  flight. 

As  will  become  clear,  this  is  not  always  an  easy  task 
because  each  of  the  four  limiting  factors  mentioned 
above  require  considerable  basic  input  data,  and  the 
URLIM  user  may  spend  considerable  effort  in  gathering 
these  required  data.  The  first  item  that  should  be 
determined  is  the  description  of  the  structure  to  be 
modeled.  Previous  work  (Ref.  1)  has  indicated  that  the 
position  on  a radome  where  the  aerodynamic  flow  changes 
from  laminar  to  turbulent  is  the  most  severe  from  a heat 
transfer  standpoint.  Therefore,  this  station  on  the 
particular  radome  of  interest  should  be  modeled  for 
calculations  of  boresight  error,  melting,  and  thermal 
stress.  It  is  generally  permissible  to  perform  only  a 
one-dimensional  heat  transfer  analysis  at  this  tran- 
sition point,  but  a two-  or  three-dimensional  analysis 
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can  be  done  if  desired.  In  general,  a multidimensional 
analysis  will  significantly  increase  the  amount  of  input 
required  and  has  been  shown  to  offer  a negligible 
increase  in  accuracy  (c.f.  Ref.  1,  Sections  3 and  4)  . 

The  heat  transfer  analysis  section  of  URL  in  requires 
that  the  structure  being  analyzed  be  divided  into  lumped 
masses  or  nodes,  within  which  temperatures  will  be 
considered  constant  at  any  specific  time.  For  the  sake 
of  accuracy,  the  radome  wall  thickness  at  the  transition 
point  should  therefore  be  divided  into  several  small 
nodes.  However,  as  the  size  of  the  nodes  is  decreased, 
the  amount  of  time  required  by  URLIH  for  the  trajectory 
simulation  increases  so  that  a compromise  must  be 
reached  between  accuracy  and  program  execution  time. 
Volume  1 of  this  report  provides  information  on  how  to 
calculate  an  appropriate  node  size. 

To  complete  the  description  of  the  radome,  the 
thermal  transport  and  mechanical  properties  of  the 
radome  material  must  be  provided.  These  properties 
include  thermal  conductivity,  specific  heat,  infrared 
emissivity.  Young's  modulus,  Poisson's  ratio,  and  free 
thermal  expansion.  Whenever  possible  these  properties 
should  be  given  as  functions  of  temperature  up  to  the 
melting  point  of  the  material. 

Figure  2a  is  a sketch  of  the  radome  geometry  selec- 
ted in  the  sample  problem  that  is  discussed  throughout 
this  report.  The  figure  shows  the  specific  data  that 
were  mentioned  generally  above.  The  tigure  shows  the 
general  missile  geometry,  the  specific  location  on  the 
radome  for  the  limit  analysis,  and  the  material  proper- 
ties that  are  required.  All  of  the  data  shown  in  this 
figure  must  be  determined  before  URLIH  is  run;  i.e., 
these  data  are  part  of  the  required  input.  The  radome 
wall  is  shown  divided  into  six  thermal  nodes  with  a 
seventh  node  of  zero  volume  at  the  exterior  surface. 

This  zero-volume  (or  surface)  node  is  a requirement  of 
the  aerodynamic  heating  calculations.  As  will  be  descri- 
bed in  detail  later,  the  definition  of  the  thermal 
interconnections  between  the  nodes  can  be  accomplished 
by  using  standard  geometries  or  by  manual  input  of  the 
conduction  paths.  The  details  of  this  data  entry  will  be 
explained  later.  Figure  2c  shows  some  of  the  detail 
required  in  specifying  the  geometry  of  the  thermal 
models.  The  models  shown  in  this  figure  are  for  the 
sample  case  described  throughout  this  report. 
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Von  Karman  profile  with  fineness  Missile  characteristics 

ratio,  2.1;  base  radius,  rb;  axial  and  CjviMand  maximum 
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Fig.  2a  Sketch  of  Some  Required  Information  for  an  URLIM  Run 
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After  the  geometry  and  material  definitions  of  the 
radome  have  been  decided,  the  trajectories  that  will  be 
used  must  be  specified.  Since  the  specification  of  the 
velocity  limits  is  the  locus  of  points  on  all  possible 
trajectories  at  which  the  radome  fails,  an  infinite 
number  of  trajectories  could  be  used.  It  remains, 
therefore,  for  the  URLIM  user  to  examine  the  flight 
potential  of  the  particular  missile  and  select  a man- 
ageable number  of  trajectories  that  will  represent  the 
range  of  possible  flights.  A single  trajectory  specifi- 
cation requires  the  velocity  and  altitude  as  a function 
of  time.  To  complete  the  trajectory  information  two 
additional  tables  must  be  supplied  --  the  thermal  and 
mass  transport  properties  of  air  as  functions  of  tem- 
perature and  pressure,  and  the  temperature  and  pressure 
of  the  atmosphere  as  a function  of  altitude. 

In  the  sample  problem  in  this  report,  two  line-of- 
sight  trajectories  were  flown  at  five  different  accel- 
lerations.  The  altitudes  achieved  were  governed  by  the 
launch  angle  of  the  trajectory  but  were  limited  to  a 
maximum  of  100  000  ft.  Figure  2b  shows  the  10  different 
trajectories  that  were  used  and  in  so  doing  indicates 
the  type  of  trajectory  information  that  is  required  by 
UBLIM. 

With  the  radome  and  trajectory  completed,  the  user 
must  decide  what  limits  he  wants  to  calculate.  The 
limits  mentioned  earlier  (thermal  stress,  melting, 
boresight  error,  and  aerodynamic  load)  may  be  selected 
in  any  combination  for  evaluation  during  a run.  Also, 
if  more  than  one  limiting  factor  is  being  run  at  one 
time,  the  user  may  select  whether  to  terminate  the 
simulation  when  any  one  limit  is  reached  or  to  continue 
through  the  trajectory  until  all  possible  limits  have 
been  met.  This  latter  mode  of  limit  calculation  will 
provide  all  of  the  possible  information  about  the  radome 
material  but  may  take  considerable  additional  computer 
time  to  execute.  Also,  the  calculations  may  be  meaning- 
less; that  is,  it  would  not  be  helpful  to  find  out  that 
the  boresight  error  rate  exceeded  its  permissible  level 
some  time  after  the  radome  had  melted  through.  More- 
over, deciding  which  limit  mode  to  use  will  depend  on 
several  items.  Some  of  these  items  may  be  known  a 
priori;  e.g. , how  well  the  radome  material  limitations 
are  known  and  how  strictly  these  limits  must  be  adhered 
to.  Other  less  easily  known  items  involve  the  interac- 
tions of  material  capability  and  the  severity  of  the 


22 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

laurel  Maryland 


selected  trajectory.  In  general  the  decision  between 
the  "complete  limit  analysis"  mode  and  the  "lowest 
limit"  mode  will  be  made  through  experience. 

The  general  procedures  outlined  above  are  what  must 
be  followed  to  generate  a radome  material  limits  plot. 
Other  types  of  aerodynamic  analyses  can  be  performed  by 
URLIM  because  the  program  has  been  designed  with  con- 
siderable flexibility.  That  is,  the  general  analysis 
capabilities  of  URLIM  can  be  selected  in  combinations  to 
perform  any  of  a variety  of  tasks.  For  example,  the 
very  simple  problem  of  determining  the  one-dimensional 
transient  response  of  a structure  with  one  boundary 
exposed  to  an  independently  specified  heat  flux  can  be 
easily  modeled  wth  URLIM.  If  the  structure  is  a cylin- 
der, then  it  would  be  possible  to  have  the  thermal 
stresses  associated  with  the  particular  heat  flux 
calaculated  and,  if  desired,  plotted  in  history  form  by 
a CalComp  plotter.  On  the  other  end  of  the  spectrum  of 
capabilities  for  URLIM,  it  is  possible  to  do  more  than 
one  limit  calculation  at  one  time.  As  many  as  nine 
materials  can  be  evaluated  over  any  number  of  trajec- 
tories and  the  resulting  limits  plotted  in  any  number  of 
combinations.  In  practice,  it  has  been  observed  that  the 
limitations  in  URLIM's  use  lie  in  the  ability  of  the 
user  to  formulate  the  problem  rather  than  in  limitations 
of  capacity  within  the  program.  This  general  capability 
of  URLIM  must  be  carefully  considered  by  the  user;  i.e., 
is  it  better  to  try  to  solve  a large  limits  problem  with 
one  run  of  URLIM  or  break  the  problem  down  into  smaller 
units  that  are  run  one  at  a time?  This  question  of  how 
to  apply  URLIM's  capability  will  be  answered  through 
experience. 

In  the  section  that  follows  the  required  inputs  to 
URLIM  will  be  described.  Within  this  section  the 
organization  of  the  basic  analysis  units  of  URLIM  will 
be  shown  along  with  the  essential  data  flow  paths.  This 
will  be  accomplished  via  a general  block  diagram  of  the 
URLIM  program.  Following  this  introductory  description 
will  be  a specific  card-by-card  discussion  of  the  re- 
quired input. 
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2.  PROGRAM  ORGANIZATION 


The  ORLIM  program  can  simulate  missile  flights 
through  the  atmosphere  and  in  doing  so  requires  a con- 
siderable amount  of  basic  input  data.  Figure  3 depicts 
the  overall  internal  structure  of  the  URLIM  program  and 
shows  many  of  its  essential  data  flows.  At  the  top  of 
the  figure  is  the  section  that  deals  with  the  reading  in 
of  the  required  data.  The  first  group  of  required  data 
is  that  which  defines  the  thermal  transport  properties 
and  dielectric  constant  of  the  radome  material.  Since 
the  ORLIM  program  is,  at  its  most  basic  level,  a tran- 
sient thermal  response  program,  the  thermal  properties 
of  the  material  are  essential.  In  order  to  calculate 
thermal  stresses  the  mechanical  properties  of  the 
material  are  required  and  are  shown  as  the  next  block  of 
required  data. 

Since  a flight  through  the  atmosphere  is  being 
simulated,  the  thermal  and  mass  transport  properties  of 
air  must  be  available.  The  third  block  at  the  top  of 
Fig.  3 represents  the  reading  in  of  these  data,  which 
are  assumed  to  be  functions  of  temperature  and  pressure. 
The  next  block  of  data  completes  the  description  of  the 
atmosphere  by  tabulating  the  temperature,  pressure,  and 
speed  of  sound  as  functions  of  altitude.  The  next  block 
of  data  provides  the  means  for  calculating  the  local 
aerodynamic  conditions  that  exist  at  the  precise  loca- 
tion of  interest  on  the  radome  by  tabulating  the  pres- 
sure ratio,  Mach  ratio,  and  velocity  gradient  term  (for 
stagnation  points  only)  as  functions  of  freestream  Mach 
number.  The  following  block  shows  the  input  of  some  of 
the  aerodynamic  characteristics  of  the  radome  and  mis- 
sile. These  values  are  required  by  the  aerodynamic  load 
routine  so  that  maneuver  and  pressure  forces  on  the 
radome  may  be  calculated  and,  in  turn,  the  attachment 
area  stresses  ca'n  be  evaluated.  The  last  block  of  input 
data  on  the  top  line  of  Fig.  3 provides  for  specifying 
the  trajectory  in  terms  of  altitude  and  velocity  versus 
time.  In  this  block  of  data,  there  may  be  other  time- 
dependent  parameters  that  would  be  required  by  the 
particular  run  of  URLIM.  Examples  of  such  inputs  are 
shown  in  the  form  of  externally  supplied  temperatures 
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and/or  heat  fluxes,  along  with  the  maximum  permissible 
angle  of  attack  for  the  missile. 

The  second  half  of  the  input  section  involves  the 
geometry  definitions  that  are  required  to  model  the 
flight  structure.  The  analysis  technique  tor  heat  trans- 
fer requires  that  the  structure  be  divided  into  smaller 
elements  (nodes)  and  that  the  thermal  interconnection 
between  nodes  be  defined.  In  the  cases  where  the  struc- 
ture may  be  described  by  cylindrical,  spherical,  or  two- 
or  three-dimensional  Cartesian  coordinates,  the  appro- 
priately named  data  input  blocks  can  be  optionally  used. 
If  the  geometry  is  arbitrary  and  cannot  be  simply 
described,  then  the  nodal  network  information  may  be 
input  via  the  READCP  input  block.  The  first  block  in  the 
geometry  section  shows  that  the  node  numbers  must  be 
given  as  functions  of  radius  tor  the  purpose  ot  cal- 
culating the  thermal  stresses.  This  information  is 
required  by  the  thermal  stress  routine  that  uses  a thick 
wall  cylindrical  approximation  for  the  radome  shape. 

The  trajectory  simulation  section  of  the  program 
determines  the  temperature  histories  of  ach  of  the 
nodes  described  by  the  geometry  data  via  an  explicit, 
forward  marching,  finite  difference  technique.  The 
boundary  conditions  specified  by  the  input  data  can 
include  the  following  time-dependent  functions:  velocity 
and  altitude  histories,  independently  specified  nodal 
heat  fluxes,  or  temperatures.  All  of  these  optional 
specifications  come  together  in  the  program  at  the  block 
called  CAPCON/STEP  wherein  the  finite  difference  solu- 
tion to  the  heat  conduction  problem  is  obtained.  The 
results  of  this  trajectory  simulation  are  the  nodal 
temperatures  that  vary  with  time.  This  information  is 
combined  with  the  flight  profiles  in  the  Limits  section 
to  provide  the  thermal  stresses,  the  aerodynamic  load 
information,  and  the  boresight  error  rates.  During  a 
run  of  the  URLIM  program  the  four  limits  (including 
melting)  are  monitored;  when  excessive  values  are 
reached  they  are  noted  by  the  LIMIT  block  routines  and 
saved  for  later  printout. 

The  last  section  of  Fig.  3 is  concerned  with  the 
program’s  output.  The  first  output  is  the  error/warning 
messages  that  may  be  generated  during  the  execution  for 
any  of  a variety  of  reasons.  The  next  output  is  the 
printer  generated  plots  of  the  combined  limits  data, 
samples  of  which  are  given  in  Figs.  4 and  5.  These 
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Fig.  5 Example  of  a Total  Limits  Plot  Generated  by  the  URLIM  Program 
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plots,  which  were  run  for  a sample  problem  that  will  be 
discussed  below,  show  the  set  of  trajectories  examined 
plotted  (with  asterisks)  as  velocity  versus  time.  The 
letters  printed  on  the  trajectories  indicate  the  occur- 
rence of  one  of  the  four  limits.  The  locus  of  the 
lowest  of  these  points  delineates  the  permissible 
operating  regime  for  the  particular  radome  modeled 
during  the  run.  Although  this  output  is  somewhat  inac- 
curate due  to  the  method  of  plotting  (line  printer)  it 
is  a very  handy  and  concise  way  of  reading  the  results 
of  the  URLIM  run.  The  output  bloc*  WRITE  prints  all  of 
the  other  data  as  it  is  calculated.  In  the  case  of  a 
very  large  and/or  very  long  run,  this  output  may  be  put 
on  microfiche.  The  final  output  block  provides  the 
capability  for  having  off-line  CalComp  plots  made  of 
selectable  parameters  (e.g.,  temperatures,  stresses, 
boresight  errors,  electrical  wall  thickness  changes,  or 
mechanical  load  stresses)  as  functions  of  time. 

Sufficient  generality  has  been  built  into  URLIM  to 
allow  virtually  any  type  of  thermal  analysis.  For 
example,  as  many  as  nine  different  materials  may  be 
flown  on  any  number  of  user-defined  trajectories  and  any 
combination  of  the  four  basic  limits  may  be  determined. 
Such  a run  could  use  all  of  URLIM's  capability.  On  the 
other  hand  it  may  be  required  simply  to  take  some  ex- 
perimental data,  such  as  front  face  temperature  versus 
time  on  a cylinder,  and  to  determine  the  thermal  stres- 
ses that  result.  Such  a problem  is  easily  modeled  with 
URLIM  and  the  CalComp  plotter  routines  can  provide  the 
required  thermal  stress  history  in  graphical  form.  In 
summary,  the  URLIM  program  can  be  used  to  almost  any 
degree  of  complication  because  its  basic  computational 
units  are  under  the  user's  control  to  be  selected  as 
required. 

The  sample  problem  that  will  be  used  for  illustra- 
tive purposes  in  the  next  section  was  assembled  to 
demonstrate  the  versatility  of  the  URLIM  program.  The 
problem  involves  determining  the  combined  velocity 
limits  of  three  materials  (alumina,  fused  silica,  and 
Pyroceram  9606)  over  the  10  trajectories  mentioned 
earlier  (c.f.  Fig.  2b).  A K-band  radar  frequency  of  30 
GHz  was  used  for  purposes  of  wall  thickness  calcula- 
tions. The  three  lesultant  nodal  breakups  are  shown  in 
Fig.  2c. 
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It  should  be  noted  that  simulating  the  ten  trajec- 
tories for  each  of  the  three  Materials  implies  running  a 
total  of  30  separate  thermal  Models.  It  is  therefore 
necessary  to  create  30  separate  and  uniquely  numbered 
groups  of  thermal  nodes  and  fly  each  of  them  on  its 
respective  trajectory.  With  the  basic  scope  of  the 
sample  problem  now  stated,  a description  of  the  precise 
inputs  required  to  solve  this  problem  is  given. 
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3.  REQUIRED  INPUTS 


Before  an  DRUM  run  can  be  made  the  user  must  have 
(1)  an  IBM  System  360  or  370  computer  available  tc  him 
and  (2)  a copy  of  the  completely  executable  URLIM  code. 
Given  these  two  essentials  the  user  must  assemble  a 
series  of  job  control  language  (JCL)  statements  and  data 
cards  that  will  access  the  URLIM  program  via  the  operat- 
ing system  (OS)  at  the  computer  center.  The  examples  of 
JCL  that  are  included  with  this  report  will  be  ones  that 
apply  to  the  360/91  computer  at  the  Applied  Physics 
Laboratory  using  CS/MVT  under  control  of  ASP  (The  Asym- 
metric Multi-processing  system).  Further  descriptions  of 
OS  and  JCL  are  found  in  "Interface  with  the  Operating 
System",  in  Section  3 of  Volume  1. 

Figure  6 shows  the  sequence  of  JCL  used  to  run  the 
sample  case  that  will  be  referred  to  throughout  -this 
report.  For  the  moment,  we  will  concentrate  only  on  the 
first  few  cards  that  deal  with  accessing  the  URLIM 
program.  The  first  card  is  the  JOB  card  and  contains 
user  and  job  identification  information  along  with 
accounting  data  that  is  generic  tc  the  APL  computer 
center.  The  second  JCL  statement  is  the  execute  (EXEC) 
statement: 


// 

EXEC 

PGM= URLIM, 

// 

REGION=240K,TIME=10, 

// 

PARK* 'ISASIZE  (6  8K)  , REPORT 

These  cards  specify  to  the  operating  system  that  the 
program  URLIM  is  to  be  executed,  that  it  will  require 
240  000  bytes  of  storage  and  will  run  for  no  more  than 
10  minutes.  Additionally  a parameter  is  passed  to  the 
program  on  this  card  specifying  an  Initial  Storage  Area 
(ISASIZE)  of  66  000  bytes  and  that  a report  is  ashed 
for.  These  two  parameters,  which  are  important  to 
efficient  execution  of  the  program,  will  be  described  in 
detail  later  (c.f.  "Dynamic  Storage"  in  Section  3 of 
Volume  1).  The  third  JCL  statement  indicates  the  data 
set  that  contains  the  program  URLIM: 
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//KELLY# 1 

JOB 

(11116,  C,IJ,N)  ,'  KELLY  FRAZE?' , NOTIP Y=EFAZE? 

//  EXEC 

PGH=URLIH, 

// 

REGiON=240K, 

// 

71 M7= 1 C , 

// 

PARH=* ISASTZE  (63K)  , REPORT* 

//STEPLIB 

DL 

CS N=  B B E . 7 R AZ  EP . U R 1 1 ?.  A LL , DI S P = 3 K * 

//SYsrar nt 

i)D 

SY  30UT-A  , 

// 

CCB*  (RECEM=VEA,LR’!CL=-137,BLK3IZ  E = 3 13  3) 

//PLIDUflP 

DD 

S Y SO  (J  T = A 

//SYSIN 

DD 

csn*bbe.fr\zfp.5.m«pl:n,disf=3hf 

// 

DD 

* 

• 

/* 

//FLO RAT 5 

DC 

CSN=!3BE.  FRAZER.  F!li 4 DATA  (VK2 114)  , DISP=SHR 

//FORCOEF 

n 

CSN=PRE.FRAZFR. BHWDATA (AFCCFT) ,EI3P=SHP 

//ATMOS 

DC 

C3N=P3E. FRAZER. HHMHATA (ATWC2) ,DISP=3H3 

//AIR 

DD 

C3N  = EBE.FBAZER.«H«DATA(MB7X2S)  ,DISP=3I>? 

//LIfIPLT 

pr 

CSN=BBE. FRAZER. LI MPLT01 , 

// 

SPACE*  (3124,  (2,0)  , R L 3 E)  , 

// 

DISP=  (NFM, CAT  LO)  , 

// 

UNIT*  SAVE 

//MAXPKNT 

OP 

SYSJUT=A, 

// 

CCP= (RECPM  = V E A , L°ECI=137,ELKSIZF=3135) 

Fig.  6 JCL  Statements  Required  to  Run  Sample  Problem 
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//STEPLIB  DD  DSN=BBE. F FAZEE .URLIMALL , 

//  DISP=SHF 

The  data  set  name  (DSN)  on  this  STEPLIB  card  indicates 
that  the  partioned  data  set  BBE. F RAZEE. URLIHALL  is  the 
data  set  that  must  contain  a member  named  URLIM.  The 
fourth  JCL  statement  describes  the  basic  output  file 
SYSPRINT  and  shows  via  the  SYSOUT=A  specification  that 
the  output  will  be  on  standard  computer  paper.  A 
provision  has  been  made  to  receive  this  output  on 
microfiche  and  is  used  by  simply  specifying  SYSGUT=Q 
rather  than  SYSOOT=A  for  the  particular  print  file  in 
question.  (This  option  is  available  in  this  way  at  the 
APL  computer  center) . The  other  JCL  statements  all 
describe  input  or  output  files  that  are  needed  by  the 
program  during  execution;  their  application  will  become 
evident  as  the  required  input  is  further  described.  In 
general,  the  required  input  data  can  be  read  from  the 
standard  input  file  (SYSIN)  or  from  other  files  of  the 
user's  specification.  (Appendix  A gives  a description  of 
the  data  input  formats  that  are  used  by  the  program  and 
referred  to  below.)  The  first  data  cards  that  are  read 
by  the  UFLIM  program  are  basic  control  values  that  allow 
the  user  to  specify  what  analysis  will  be  made.  These 
values  are  read  from  the  SYSIN  file  only  and  are  expec- 
ted to  be  in  the  PL/I  DATA  format;  that  is,  the  variable 
values  are  given  in  the  form 

variable  name  = value 

in  sequence  seperated  by  a comma  cr  one  or  more  blanks 
until  a semicolon  is  encountered.  All  of  the  possible 
control  variables  have  specific  default  values,  as  shown 
in  Table  1.  The  significance  of  these  variables  and 
their  values  will  become  clear  in  the  subsequent  ex- 
planations. Because  they  must  be  the  first  data  input  to 
the  program  they  are  given  here  with  brief  explanations. 
If  no  value  is  specified  for  one  of  these  variables, 
then  the  default  value  will  apply. 
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Table  1 

UPLIM  Control  Variables 


Variable 

Default 

Value 

Use 

CALL A BO 

0 

Number  of  separate  aerodynamic 
heating  blocks 

CALLLOD 

0 

Number  of  individual  load  cal- 
culations 

CALLSIG 

0 

Number  of  separate  thermal 
stress  calculations 

CALLBSE 

0 

Number  of  separate  boresight 
error  calculations 

CALLCYL 

0 

Number  of  thermal  networks  that 
are  two-dimensional  cylindrical 

CALLCL3 

0 

Number  of  thermal  networks  that 
are  three-dimensional  cylindri- 
cal 

CALLFLO 

0 

Number  of  separate  aerodynamic 
flow  tables  to  be  read  in 

CALLLIN 

0 

Number  of  combined  limit  plots 
to  be  made 

CALLTOD 

0 

Number  of  nodal  networks  that 
are  two-dimensional  Cartesian 

CALLTIH 

0 

Number  of  separate  time  tables 
to  be  read  in 

CALLFLX 

0 

Number  of  heat  flux  boundary 
conditions 

CALLFRC 

0 

Number  of  fixed  temperature 
boundary  conditions 

UANLIN 

0 

Flag  indicating  manual  input  of 
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Variable 

HAXENTS 

HCPROFS 

MINPRNT 

AIRPRNT 

BINRUN 

PLOTS 

THPROPS 

XCAPLIN 

tMACH 

INODES 

#SET 

STARTIH 

STOPTIH 


Default 

Value  Use 

nodal  parameters:  0 = NO  INPUT, 

1 = INPUT  ON  NETFILE 

0 Maximum  number  of  entries  to  the 

mechanical  property  tables 

0 Number  of  mechanical  propel ty 

tables 

0 Flag  for  maximum  print  option: 

1 = minimum  print,  0 = regular 
print,  2 = minimum  and  maximum 
print 

0 Flag  for  printing  the  air 

property  tables:  1 = print  them 
on  SYSPRINT,  0 = no  print 

0 Flag  for  minimum  execution 

option:  1 = minimum  run;  0 = 
maximum  run 

0 Flag  for  CalComp  plots;  0 = no 
plots,  1 = plot  information  will 
be  saved 

1 Number  of  material  thermal 
property  tables 

20  Highest  thermal  node  number 

0 Number  of  Mach  number  entries  to 

aerodynamic  force  tables 

0 Total  number  of  thermal  nodes 

0 Number  of  nodes  to  be  initial- 

ized to  a set  temperature 

0 The  initial  trajectory  simula- 

tion time  (s) 

10  The  final  trajectory  simulation 

time  (s) 
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Variable 

Default 

Value 

(Jse 

SIGLIST 

0 

Flag  for  specifying  input  format 
for  radius-node-temperature  data 
for  thermal  stresses 

TINIT 

0 

The  initial  temperature  for  all 
nodes  except  those  set  specifi- 
cally 

# ALPHA 

0 

Number  of  angle  of  attack  aata 
points  in  aeroforce  tables 

TAIRFACT 

1 

Converstion  factor  for  air  tem- 
perature values  to  convert  to  °R 

HAIRFACT 

1 

Conversion  factor  for  air 
enthalpy  tables  to  convert  to 
Btu/lbm 

MAXSTEP 

1 

Maximum  allowable  time  step  (s) 

HIM  STEP 

0.05 

Minimum  allowable  time  step  (s) 

PALTFACT 

2115.36 

Conversion  factor  for  pressure 
versus  altitude  tables  to  con- 
vert to  lbf/ft2 

PAIRFACT 

2115.36 

Conversion  factor  for  pressure 
variable  in  air  property  tables 
to  convert  to  lbf/ft2 

METFILE 

• SYSIN • 

File  name  for  specified  nodal 
information 

AFCFILE 

• SYSIN* 

File  name  for  aerodynamic  force 
data 

TREF 

536 

Reference  temperature  for  all 
enthalpies  (°R) 

TALTFACT 

1 

Conversion  factor  for  temperat- 

ures  in  the  altitude  tables  to 
convert  to  °F 
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In  the  second  major  input  section  the  material 
properties  are  read  in.  The  thermal  transport  proper- 
ties are  read  in  first  for  each  material  followed  by  the 
mechanical  properties  (elastic  modulus,  Poisson's  ratio, 
free  thermal  expansion).  The  formats  are  described  in 
the  following  paragraphs. 


THERMAL  PROPERTIES 


There  must  be  a number  of  data  groups  equal  to  the 
value  of  the  control  variable  THPROPS  with  the  rollowing 
structure: 


1.  On  the  SYSIN  file  there  must  be  a card  (or 
cards)  in  LIST  format  with: 


a.  a file  name. 


b.  the  number  of  table  entries  to  be 
read,  and 

c.  the  maximum  allowable  boresight 
error  for  this  material. 


2.  On  the  file  named  in  a (above)  there  must  be 
data  values  that  will  be  read  as:  two  comment 
cards  (160  characters)  followed  by  the  number  of 
cards  specified  in  b (above)  , with  the  following 
information  on  each: 


Cols  Value 

1-14  Temperature  (°R)  , F-format; 

15-19  Density  times  specific  heat 

(Btu/ft3-°R)  , F-format; 

30-44  Thermal  conductivity  (Btu/ft-h 
°R)  , F-format; 

45-59  Infrared  emissivity,  F-format; 
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Cols  Value 

60-69  Dielectric  constant,  F-format; 

70-80  Loss  tangent,  F-format. 

Material  code  numbers  starting  with  1 and  running 
serially  to  THPROPS  will  be  associated  with  each  group 
of  thermal  properties.  In  Appendix  B the  sample  input 
shows  that  the  thermal  properties  of  three  materials 
were  read  from  the  SYSIN  file.  The  sample  output  in 
Appendix  C shows  that  the  three  material  codes  were 
assigned  in  order  and  that  all  of  the  input  values  were 
printed  out  and  reformated.  It  is  important  to  note 
that  the  last  temperature  value  listed  in  each  thermal 
property  table  is  used  as  the  melt  temperature  for  the 
particular  material.  Also,  the  temperatures  must  be 
listed  in  monotonically  increasing  order. 

The  thermal  properties  (with  RanJeine  temperature 
scale  dependence)  are  followed  by  the  mechanical  proper- 
ties with  Farenheit  temperature  scale  dependence. 


MECHANICAL  PROPERTIES 


There  must  be  a number  of  data  groups  equal  to  the 
value  of  the  control  variable  MCPROPS  with  the  following 
format : 


1.  On  the  SYSIN  file  a card  (or  cards)  in  LIST 
format  with: 


a.  a file  name, 

b.  the  number  of  table  entries,  and 

c.  the  maximum  allowable  bending 
stress  (the  modulus  of  rupture)  . 

2.  In  the  file  named  in  a (above)  there  must  be  as 
many  values  of  temperature,  modulus,  Poisson's 
ratio,  and  free  thermal  expansion  as  specified 
in  b (above)  with  the  following  format:  two 
comment  cards  (160  characters)  followed  by  data 
cards: 
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Cols 

Va  lue 

1-14 

Values  of 
f orma  t; 

temperature  (°F)  , F- 

15-29 

Values  of 
E-format; 

elastic  modulus  (psi)  , 

30-44 

Values  of 
E-f or  mat ; 

Poisson's  Ratio  (psi)  , 

45-59 

Values  of 

free  thermal  expansion 

(in. /in.)  , E-format.. 


The  sample  input  shows  that  the  three  mechanical 
property  tables  are  read  from  the  SYSIN  file.  The 
sample  output  in  Appendix  C shows  that  the  three  tables 
were  read  in  and  successive  mechanical  property  code 
numbers  were  assigned.  As  before,  the  temperatures  must 
be  in  increasing  monotonic  order. 


CYLINDRICAL  GEORETRIES  --  TWO- DIMENSIONAL 


The  next  several  cards  in  the  sample  input  section 
are  the  control  cards  for  obtaining  automatically 
generated  cylindrical  geometries.  The  control  variable 
CALLCYL  controls  the  number  of  such  nodal  definitions  to 
be  made.  In  the  sample  case,  CALLCYL  was  set  to  30  in 
the  initial  input  so  there  are  30  sets  of  data  on  the 
SYSIN  file  in  this  section,  each  with  the  following  data 
in  LIST-format  (c.f.  Fig.  7): 

Item  Value 


1 NR,  the  number  of  radial  divisions  for  the 
cylindrical  section; 

2 NTH,  the  number  of  circumferential  divisions 
in  the  section; 

3 ROOT,  the  outer  radius  of  the  section  (ft)  ; 
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Side  6 


Notes: 

• All  dimensions  must  be  in  feet  or  degrees 

• Node  numbering  begins  at  upper  left  corner  {0  = 0,  R = Rout)  and 
counts  first  radially  and  then  circumferentially. 

• The  node  centers  (small  dots)  are  used  for  computation  of  the 
thermal  path  lengths  between  nodes  and  are  located  at  the  geometric 
center  of  gravity  except  for  boundary  or  side  nodes.  When  the  control 
valve  is  set  to  1 for  a side  then  the  node  centers  are  considered  to  be 
at  the  external  surface  for  computations  of  thermal  path  length. 

• List  of  angular  and  radial  spacings  can  be  arbitrary  or  equally  spaced 
according  to  values  input  (see  text  for  explanations). 


Fig.  7 Cylindrical  Section  Geometry  Conventions 
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Item  Value 

4 z,  the  depth  of  the  section,  Z-axis  dimension 
(ft)  ; 

5 Thermal  property  code  number  for  the  material 
in  this  section; 

6 The  node  number  for  the  first  node  in  the 
array  of  nodes  (c.f.  Fig.  7 for  the  convention 
used  to  number  the  nodes)  ; 

7-10  Four  numbers  (either  0's  or  I's)  that  indicate 
the  type  of  nodal  connection  to  use  on  each  of 
the  four  sides  of  the  cylindrical  section; 

11  A list  of  NR  numbers  specifying  the  radial 

thickness  tc  be  used  in  the  section  starting 
from  the  surface  and  proceeding  inward  (if  NR 
is  1,  then  the  first  number  input  here  is 
taken  to  be  the  total  radial  thickness  of  the 
section;  it  is  followed  by  a second  number 
that  indicates  the  number  of  evenly  spaced 
radial  divisions  to  be  used  for  the  section; 
length  or  radius  values  are  in  feet) ; and 

Next  A list  of  NTH  numbers  specifying  the  angular 
divisions  to  be  used  in  the  circumferential 
direction  (deg)  (if  NTH  is  set  to  1,  then  the 
first  number  here  is  the  total  angle  of  the 
section  and  the  next  input  value  is  the  number 
of  evenly  spaced  divisions  that  are  required). 

The  sample  input  of  Appendix  B shows’  the  30  groups  of 
data  that  were  input  and  the  sample  output  shows  the  30 
cylindrical  node  geometries  that  were  generated.  Note 
that  the  surface  nodes  are  not  generated  in  this  com- 
putation (they  are  input  later  via  the  manual  nodal 
input).  This  is  not  tc  say  that  these  surface  nodes 
could  not  have  been  input  through  this  section.  To  do 
so  the  nodal  thickness  is  simply  entered  as  zero.  This 
method  was  not  done  in  the  sample  problem  so  as  to 
illustrate  the  use  of  the  manual  input  feature  described 
in  a subsequent  paragraph. 
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CYLINDRICAL  GEOMETRIES  --  THREE-DIMENSIONAL 


If  the  control  variable  CALLCL3  is  set  to  a non- 
zero value  then  that  number  of  data  groups  will  be  read 
from  the  SYSIN  file.  These  data  will  describe  three- 
diaensional  cylindrical  node  geometries.  The  required 
inputs  for  this  automatic  network  generator  are  very 
similar  to  those  for  two-dimensional  networks  as  de- 
scribed above.  Specifically,  for  each  network  there 
must  be  the  following  data  in  list  format  on  the  SYSIN 
file:  (c.f.  Fig.  7) 

Item  Value 


1 NR,  the  number  of  divisions  to  be  made  in  the 
radial  direction; 

2 NTH,  the  number  of  divisions  to  be  made  in  the 
circumferential  direction; 

3 NZ,  the  number  of  divisions  to  be  made  in  the 
axial  direction; 

U ROOT,  the  outer  radius  of  the  secton; 

5 Thermal  property  code  for  the  material  com- 
prising this  section; 

6 The  node  number  of  the  first  node  in  this 
section ; 

7-12  Six  numbers  (either  0's  or  I's)  that  indicate 
the  type  of  conductance  path  to  be  used  on 
each  of  the  section's  six  sides; 

13  A conversion  factor  that  will  be  multiplied 
times  each  input  variable  that  is  a spatial 
dimension  so  as  to  obtain  units  of  feet  as  a 
result; 

14  A list  of  NR  numbers  specifying  the  series  of 
radial  thicknesses  to  be  used  (this  item  is 
interpreted  identically  to  the  list  of  thick- 
nesses used  by  the  two-dimensional  routine 
(above)  ) ; 
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Item  Value 

Next  A list  of  NTH  numbers  specifying  the  circum- 
ferential divisions  to  be  used  in  the  section 
(deg)  (this  list  is  interpreted  the  same  as 
the  list  of  circumferential  divisions  used  by 
the  two-dimensional  cylinder  routine  (above)); 

Next  A list  of  NZ  numbers  specifying  the  axial 

thicknesses  to  be  used  in  this  section  (this 
list  is  interpreted  in  the  same  manner  as  the 
two  previous  lists) . 

The  units  of  the  linear  dimensions  given  above  will 
be  considered  to  be  feet  after  multiplication  by  the 
13th  value  input  (the  conversion  factor). 

The  sample  problem  presented  here  does  not  contain 
three-dimensional  cylindrical  sections  as  none  were 
required  by  the  model.  The  output  generated  by  this 
facility  is  quite  similar  to  the  previous  two-dimen- 
sional cylindrical  facility  for  which  there  is  an  ex- 
ample in  the  sample  output. 


RECTANGULAR  GEOMETRIES 


If  the  control  variable  CALLTOD  is  set  to  a non- 
zero value  then  that  number  of  data  groups  will  be  read 
from  the  SYSIN  file.  These  data  will  describe  two- 
dimensional  Cartesian  node  geometries.  The  required 
inputs  for  this  automatic  network  generator  are  very 
similar  to  those  for  the  cylindrical  geometries. 
Specifically,  for  each  two-dimensional  grid  network 
there  must  be  data  cards  in  the  SYSIN  file  with  the 
following  information:  (c.f.  Fig.  8) 

Item  Value 


NX,  tha  number  of  division  to  be  made  in  the 
x-direction ; 

NY,  the  number  of  divisions  to  be  made  in  the 
y-direction  ; 
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Notes: 

• All  dimensions  must  be  in  feet 

• Node  numbering  begins  at  upper  left  corner  (x  = 0,  y = 0)  and 
counts  first  in  x -direction  and  then  in  y-direction. 

• The  node  centers  (small  dots)  are  used  for  computation  of  the 
thermal  path  lengths  between  nodes  and  are  located  at  the 
geometric  center  of  gravity  except  for  boundary  on  side  nodes: 
when  the  control  value  is  set  to  1 for  a side  then  the  node  center 
is  taken  to  be  at  the  external  edge  of  the  node  for  computation 
of  thermal  path  length. 

• The  list  of  angular  and  radial  spacings  can  be  arbitrary,  or  equally 
spaced  according  to  input  values  (see  text  for  explanations). 


Fig.  8 Two-Dimensional  Cartesian  Node  Network 
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Item  Value 

3 z,  the  depth  of  the  nodal  grid  in  the  z-direc- 
tion  (ft)  ; 

4 Material  property  code  for  thermal  transport 
properties; 

5 The  node  number  of  the  first  node,  (c.f.  Fig. 

8 for  node  numbering  conventions) ; 

6-9  Four  numbers  (0 's  or  I's)  that  select  the  type 
of  nodal  connections  for  boundary  nodes; 

Next  & list  of  NX  numbers  specifying  the  x-direc- 
tion  thickness  to  be  used  by  the  automatic 
geometry  routine  (if  NX  is  1,  then  the  first 
number  in  this  list  is  the  total  length  of  the 
region  in  the  x-direction  and  the  next  number 
is  the  number  of  equally  spaced  nodal  divi- 
sions that  will  be  used)  (ft) ; and 

Next  A list  of  N Y numbers  specifying  the  y-direc- 
tion  thicknesses  to  be  used  by  the  automatic 
geometry  routine  (If  NY  is  1 then  the  first 
number  in  this  list  is  assumed  to  be  the  total 
length  of  the  region  in  the  y-direction  and 
the  next  number  is  the  number  of  equally 
spaced  nodal  divisions  that  will  be  used 
within  the  y-directicn  thickness)  (ft)  . 

In  the  sample  problem  the  control  variable  CALLTOD 
was  not  set,  so  the  default  value  of  zero  was  used;  i.e. 
no  two-dimensional  Cartesian  ncde  regions  were  required 
and  no  two  dimensional  cards  were  read.  However,  had 
there  been  any  cards  read,  tables  similar  to  the  ones 
for  the  cylinder  input  would  be  printed  on  the  output 
file  for  verification  of  the  input  values. 

The  current  version  of  URLIM  does  not  have  the 
spherical  or  three  dimensional  Cartesian  node  geometry 
routines  in  it.  In  a subsequent  release  of  the  URLIM 
program  these  routines  will  be  available  and  the  input 
for  them  will  occur  at  this  location  in  the  SYSIN  file 
and  will  be  controlled  with  the  variables  CALLSPH  and 
CALL3D. 
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GEOMETRY  TABLE 


If  the  control  variable  NANLIN  is  set  to  1 then  the 
file  specified  by  the  NETFILE  variable  will  be  read  for 
thermal  network  data.  This  data  must  be  in  the  follow- 
ing format: 

Two  comment  cards  (160  characters)  followed  by  a 
variable  number  of  network  data  cards  that  describe 
nodal  interconnections  in  the  following  format: 

Columns  Item 


1-4  Node  number,  F-format; 

10-19  Node  volume  (ft3).  E-format; 

20-29  Node  thermal  property  identification  number 
associated  with  thermal  property  tables  read 
in  earlier,  F-format; 

30-39  The  number  of  interconnections  between  this 
node  and  its  immediate  neighbors,  F-tormat; 

40-49  The  node  number  for  the  first  interconnec- 
tion, F-format; 

50-59  The  thermal  path  length  (area/distance) 
between  node  centers  (ft).  E-format; 

60-69  The  node  number  of  the  next  immediately  con- 
necting node  (if  present),  F-format;  and 

70-79  The  thermal  path  length  between  the  node  and 
its  next  connecting  node  (if  any)  (ft)  , E- 
format . 

There  may  be  as  many  entries  to  this  table  as 
reguired  for  the  model  and  they  may  be  in  any  order. 

The  end  of  the  list  of  cards  is  indicated  by  including  a 
card  with  columns  1-4  blank  or  a zero.  Moreover,  if  any 
datafield  is  left  blank,  a zero  will  be  assumed  for  the 
value.  The  sample  input  of  Appendix  B shows  that  the 
surface  nodes  for  all  the  aerodynamic  heating  conditions 
were  read  from  the  SY SIN  file.  These  particular  nodes 
have  zero  volume  (they  are  actually  surfaces)  and  are 
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connected  to  structural  nodes  as  indicated  by  the  data 
entries.  The  sample  output  data  in  Appendix  C shows  a 
listing  of  the  read  in  data  values,  again  slightly 
reformatted. 

The  routine  in  URLIM  that  reads  and  processes  the 
manual  geometry  input  (READCP)  has  the  capability  of 
connecting  nodes  with  radiation  heat  transfer  as  well  as 
joining  of  different  materials  through  a contact  resis- 
tance. Also,  any  one  node  may  be  described  as  having 
connections  to  any  number  of  other  nodes  by  means  of  a 
special  input  procedure.  Details  of  this  capability  are 
found  in  Appendix  F. 


AERODYNAMIC  HEATING  SPECIFICATIONS 


If  any  aerodynamic  heating  calculations  are  to  be 
made  then  the  control  variable  CALLARO  should  be  set  to 
indicate  how  many  different  aerodynamic  surfaces  or 
"blocks"  there  will  be.  (An  aerodynamic  block  is 
actually  a set  of  surface  (zero-volume)  nodes  that  are 
subjected  to  the  same  aerodynamic  condition.  Usually, 
when  there  is  more  than  one  node  to  a block,  the  surface 
nodes  are  part  of  a two-dimensional  thermal  network  that 
is  subjected  to  the  same  aerodynamic  environment  over 
its  entire  surface.)  Having  CALLARO  set  to  a nonzero 
value  will  cause  the  program  to  expect  several  groups  of 
data  to  occur  in  the  SYSIN  file  directly  after  the  last 
geometry  cards.  These  data  in  LIST  format  are: 


Item  Value 

1 A file  name  indicating  where  a table  of 
thermal  and  mass  transport  properties  of  air 
will  be  found; 

2 The  number  of  pressure  entries  in  this  table; 
and 

3 The  number  of  temperature  entries  in  this 
table. 
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If  the  file  named  in  item  1 (above)  is  SYSIN  then  the 
above  values  are  followed  immediately  with  a table  of 
air  properties  in  the  following  format:  two  comment 
cards  (160  characters)  followed  by  data  cards: 

Columns  Item 


I- 10  Value  of  pressure,  F-format; 

II- 20  Value  of  temperature,  F-format; 

21-30  Value  of  enthalpy  at  the  respective  pressure 
and  temperature,  F-format; 

31-40  Value  of  Prandtl  number  at  the  respective 
pressure  and  temperature,  F-format; 

41-50  Value  of  the  ratio  cf  specific  heats  for  air, 
F-format; 

51-60  Value  of  viscosity  (lbm/ft-s)  F-format;  and 

61-70  Value  of  compressibility  for  air,  F-format. 

If  the  file  named  in  1 (above)  is  other  than  SYSIN,  then 
the  two  comment  cards  and  the  tabular  data  cards  should 
be  on  the  file  associated  with  that  file  name,  in  the 
format  stated  above.  The  logical  records  of  that  file 
will  be  expected  to  be  at  least  80  characters.  The 
physical  units  of  the  pressure,  temperature,  ana  en- 
thalpy values  were  not  specified  in  the  definitions 
above  because  they  may  be  variable  and  may  be  chosen  by 
the  user.  However,  it  must  be  remembered  that  tne  con- 
trol variables  PAIRFACT,  TAIRFACT,  and  HAIRFACT  will  be 
multiplied  by  each  of  the  pressure,  temperature  and 
enthalpy  values  actually  read  and  the  resultant  values 
will  be  assumed  to  have  the  following  units:  pressure, 
lbf/ft2;  temperature,  °R;  and  enthalpy,  Ftu/lbm. 

The  total  number  of  air  property  data  cards  read  during 
this  phase  will  be  the  number  cf  pressure  entries  (item 

2 (above))  times  the  number  of  temperature  entries  (item 

3 (above)).  The  pressure  and  temperature  values  must  be 
in  monotonically  increasing  order;  i.e.  all  of  the  data 
cards  at  the  lowest  pressure  must  be  first  and  they  must 
have  temperature  values  in  increasing  order.  The  next 
set  of  pressure  cards  must  contain  the  same  number  of 
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temperature  cards  as  before  and  the  pressure  value  must 
be  higher.  Also,  each  set  of  temperatures  must  be  the 
same;  if  they  are  not,  a warning  message  will  be  printed 
and  the  temperature  values  from  the  first  set  of  pres- 
sure cards  will  be  used. 

The  sample  problem  input  in  Appendix  B shows  that 
the  control  card  for  this  section  of  data  names  the  tile 
"AIR"  as  containing  the  air  property  data  at  7 pressures 
and  28  temperatures.  The  sample  JCL  in  Fig.  6 shows 
that  the  file  named  AIR  is  associated  with  the  data  set 
BBE. FRAZER. MHWDATA  (AI  R7X28)  and  the  data  that  were  read 
are  shown  reprinted  in  the  sample  output  with  the  ap- 
propriate conversion  factors  already  applied. 


Atmosphere  Definitions 


The  next  set  of  input  data  values  on  the  SYSIN  file 
controls  the  reading  of  the  atmosphere  data  for  the 
aerodynamic  calculations.  These  data  are  read  if  CAL- 
LARO  is  greater  than  zero.  The  data  that  will  be  re- 
quired at  this  location  on  the  SYSIN  tile  (in  LIST 
FORMAT)  are: 

1.  A file  name,  and 

2.  The  number  of  tabular  entries  to  the  atmo- 
sphere table. 

If  the  file  name  in  item  1 (above)  is  SYSIN  then  the 
immediately  following  cards  on  the  SYSIN  file  will  be 
the  atmosphere  table  and  should  have  the  following 
format:  two  comment  cards  (160  characters)  followed  by 
data  table  cards: 

Columns  Value 


1-15  Value  of  altitude  (ft),  F-format; 

16-30  Value  of  static  temperature  at  the  correspon- 
ding altitude,  F-format; 

31-45  Value  of  static  pressure  at  the  corresponding 
altitude,  F-format; 
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Columns  Value 

46-60  Value  of  the  speed  of  sound  at  the  correspon- 
ding altitude  (ft/s),  F-format;  and 

61-75  Value  of  the  effective  space  temperature  for 
use  with  radiation  relief  from  the  aerody- 
namically  heated  surfaces  (°R) , F-format. 

If  the  file  nane  mentioned  in  item  1 (above)  is 
other  than  SYSIN  then  the  two  comment  cards  and  the 
tabular  data  cards  must  be  in  a file  associated  with 
this  other  file  name.  The  static  temperature  and  pres- 
sure values  read  by  this  section  of  UBLIM  are  multiplied 
by  the  control  variables  TALTFACT  and  PALTFACT  respec- 
tively and  the  resultant  values  will  be  assumed  to  have 
the  following  units:  temperature,  °R ; and  pressure, 
lbf/ft1 2 

In  the  sample  problem  the  control  card  for  the 
atmosphere  data  names  the  file  ATMOS  as  containing  30 
data  card  entries  describing  the  atmosphere  model  to  be 
used.  The  sample  output  (Appendix  C)  shows  the  data 
read  by  this  section  with  the  conversion  factors  ap- 
plied. The  data  were  read  from  the  data  set  BBE. FRAZER. 
H HKD AT A ( ATMOS)  which  is  associated  with  the  file  name 
ATMOS  in  the  JCL  for  the  sample  job  (c.f.  Fig.  6). 


Aerodynamic  Heating  Parameters 


The  next  data  reguired  by  the  program  are  a series 
of  data  groups  describing  the  aerodynamic  heating 
parameters.  When  the  control  variable  CALLARO  is  non- 
zero then  there  should  be  CALLARO  number  of  these  data 
groups  on  the  SYSIN  file  with  the  following  data  in  LIST 
format  (c.f.  Fig.  9 for  definition  of  terms): 

Item  Value 


1 The  number  of  surface  nodes  in  the  aerody- 

namic heating  block; 
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Freestream  condition 

Density  - 
Viscosity  - Pqo 
Temperature  - Too 
Pressure  - Poo 
Velocity  - V^ 


Shock  wave  for  compressible  supersonic  flows; 
implies  local  fluid  properties  that  will  depend 
.on  coordinates  X,  Y and  body  shape 

(See  detail  below,  left) 

/ External  surface  of  body 


■First  internal  nodes 


Body  coordinates 


Local  flow  condition  ==  Velocity  ( V j ) 


Boundary 


» L v I u A . , 

Temperature  (T.)  Vf '^er  edge 


Local  coordinates 


Boundary  layer  edge 
-Heated  surface  at  T.., 


0.99  V, 


Fluid  velocity 


Boundary 
layer  edge 

M., 


For  a heated  wall 


-For  an  adiabatic  wall 


Two  thermal  nodes  with  similar  local  flow 
at  boundary  layer  edge 

Required  specifications  include; 

• Critical  Reynolds  number  for  flow  transition 

• Trajectory  identification  number 

• Material  identification  number 

• In-depth  nodal  connections 

• Local  flow  ratio  identification  number 

• View  factor  for  radiation  relief 

• Type  of  flow  field  (flat  plate,  cone,  pipe, 
or  stagnation) 

• Surface  area  (ft2) 


Fig.  9 Terminology  for  Aerodynamic  Heating  Boundary  Conditions 


Tw  Tl  Ti 


Fluid  temperature,  Tf 


recovery 
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Item  Value 

2 A code  for  the  type  of  heat  transfer  cor- 
relation to  use.  This  nuaber  has  the  follow- 
ing legal  values: 

1 flat  plate  flow 

2 cone  flow 

3 pipe  flow 

4 stagnation  points,  cylinders  (two- 

diaensional  flow  fields) 

5 stagnation  points,  spheres  (three- 
dinensional  flow  fields)  ; 

3 The  identification  number  for  the  table  of 
trajectory  parameters  (i.e.,  velocity  and 
altitude  versus  tine)  ; these  tables  and  num- 
bers  are  read  in  at  a later  point  under  the 
control  of  variable  CALLTIM ; 

4 The  identification  number  for  the  table  of 
Hach  and  pressure  ratios  that  will  ne  used  to 
define  the  flow  conditions  at  the  edge  of  the 
boundary  layer  near  the  aero  block  (these 
data  are  read  in  later  under  the  control  of 
variable  CALLFLO) ; 

5 The  critical  Reynolds  number  (used  in  deter- 
mining whether  the  boundary  layer  flow  is 
laninar  or  turbulent)  ; 

6 A view  factor  for  radiation  relief  to  the 
space  teaperature; 

7,8  Tolerance  values  used  for  judging  the  ac- 
curacy of  a nuaerical  solution  technique 
within  the  program  (the  first  is  for  tem- 
perature calculations  and  is  expressed  in  °R; 
the  second  is  a heating  rate  value  used  to 
determine  heat  balances  across  the  shock  wave 
and  across  the  heated  surface  and  is  expres- 
sed in  Btu/ft*-s;  typical  values  for  each  are 
1 °R  and  1 Btu/ft2-s) ; 
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Item  Value 

9 An  iteration  limit  to  be  used  for  the  numeri- 
cal soluticns  mentioned  above  (if,  after  the 
specified  number  of  iterations  either  of  the 
convergence  criteria  are  not  met,  then  an 
appropriate  message  will  be  printed  and  the 
execution  will  continue;  descriptions  of  the 
solution  techniques  and  error  messages  are 
found  in  Appendix  E)  ; 

10  For  each  of  the  surface  nodes  in  the  block 
(item  1 above)  the  following  must  be  speci- 
fied: 

a.  the  surface  node  number, 

b.  the  connecting,  in-depth,  node  number, 

c.  the  surface  area  (ft*)  , 

d.  the  reference  length  (ft)  , and 

e.  the  diameter-tc-length  ratio  (this  item 
used  for  laminar  pipe  flow  calculations 
only  but  some  value  must  be  supplied; 
when  other  than  pipe  flow  is  called  for 
by  i tern  2 (above)  , this  value  is  not 
used. 

In  the  sample  input  table  (Appendix  B)  there  are  30 
cards  with  the  above  information.  The  first  such  card 
indicates  that  the  first  aero  block  will  have  one  sur- 
face node,  that  the  flow  will  be  calculated  using  flat- 
plate  correlations,  the  trajectory  information  will  come 
from  time  table  No.  1,  the  local  flow  parameters  will 
come  from  flow  table  No.  1,  the  critical  Reynolds  number 
is  1x10*,  and  the  view  factor  for  radiation  relief  to 
space  is  1.  The  tolerance  levels  for  the  numerical 
solutions  are  2°R  for  temperature  calculations  and  1 
Btu/ft2-s  for  heat  balance  calculations;  ten  iterations 
will  be  allowed.  The  node  number  of  the  single  surface 
node  is  1 and  it  is  connected  to  interior  node  number  2; 
the  surface  area  is  1 ft2,  and  the  reference  length  from 
the  stagnation  point  is  0.133  ft.  The  diameter-to- 
length  ratio  (not  actually  used)  is  specified  as  1. 
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The  sample  output  shows  the  values  that  veit  read 
in  at  this  point  printed  on  SYSPRINT  for  ready  reference 
and  verification. 


L2£&l;-U2«J2££i2iti2ns 


The  next  section  of  input  that  is  read  (if  CALLARO 
is  nonzero)  contains  the  data  that  pertain  to  the  local 
flow  conditions  in  the  vicinity  of  the  aerodynamic 
heating  blocks.  The  number  of  groups  of  data  that  are 
read  in  this  section  is  the  value  of  the  control 
variable  CALLFLO.  There  should  be  CALLFLO  number  of 
groups  of  data  in  the  following  order  on  the  file  SYSIN 
(LIST  format) : 

1.  A file  name  indicating  where  the  data  table 
values  are  located,  and 

2.  The  number  of  tabular  data  points  in  the  flow 
table 

The  next  data  read  are  the  tabular  values  of  the 
flow  parameters;  they  are  read  from  the  file  named  in 
item  1 above.  The  format,  assuming  80-character  records 
(cards),  is:  two  comment  cards  (160  characters)  followed 
by  the  number  of  tabular  data  cards  specified  in  item  2 
above  with  the  following  format: 

Columns  Value 


1-14  Values  of  freestream  Mach  number,  F-format; 

15-29  Values  for  the  ratio  of  local  body  station 
Hach  number  to  freestream  Hach  number,  F- 
format; 

30-44  Values  for  the  ratio  of  local  body  station 
static  pressure  to  freestream  static  pres- 
sure, F-format;  and 

45-59  Values  of  the  velocity  gradient  term  used 

only  in  stagnation  heat  transfer  analyses  (if 
this  space  is  left  blank,  a value  of  1 will 
be  assumed),  F-format. 
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The  sample  input  shows  the  file  'FLORATS'  to  con- 
tain 17  flow  table  entries  and  the  values  read  (from  the 
card  file  FLORATS)  are  shown  reformatted  in  the  sample 
output. 


AERODYNAMIC  FORCE  COEFFICIENTS 


The  aerodynamic  load  information  for  the  missile 
and  the  radome  are  next  read  in  if  the  control  variables 
# ALPHA  and  #M ACH  are  both  nonzero.  The  data  will  be 
read  from  the  file  identified  by  the  control  variable 
AFCFILE  (default  is  SYSIN) . The  data  should  have  the 
following  in  EDIT-f or  mat:  two  comment  cards  (160  charac- 
ters) followed  by  # ALPHA  groups  of  tabular  data  cards 
each  containing  iKACH  cards  in  the  following  format: 

Columns  Value 


1-13  Value  of  angle  of  attack,  of  the  missile 
(deg)  , F-f cimat; 

14-26  Value  of  freestream  Mach  number,  F-format; 

27-39  Value  of  the  normal  force  coefficient  of  the 
complete  missile  (CNH)  at  the  corresponding 
angle  of  attack  and  Mach  number,  F-format; 

40-52  Value  of  the  normal  force  coefficient  of  the 
radome  only  (CNR),  F-format; 

53-65  Value  of  the  axial  force  (drag)  coefticient 
of  the  radome  only  (CAP),  F-format;  and 

66-78  Value  of  the  center  of  pressure  location, 
measured  from  the  tip  of  the  radome  and 
normalized  to  the  length  of  the  radome,  F- 
format . 

The  sample  input  shows  no  input  for  this  data 
section  from  the  SYSIN  file.  Since  the  control  variable 
AFCFILE  was  set  to  ' FORCOEF' , this  file  is  read  for  the 
tabular  data.  The  sample  output  shows  the  table  of  data 
items  printed  out  for  inspection.  Note  that  the  indepen- 
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dent  variables  (angle  of  attack  and  Mach  number)  are 
both  in  monotonically  increasing  order. 


AERODYNAMIC  LOAD  PARAMETERS 


The  next  data  required  by  URLIM  is  controlled  by 
the  variable  CALLLOD,  which  tells  how  many  aerodynamic 
load  cases  will  be  evaluated.  The  calculations  using 
these  data  determine  the  missile-radome  attachment 
stresses  due  to  maneuver  and  pressure  forces.  On  the 
file  SYSIN  there  should  be  CALLLOD  groups  of  data  in  the 
following  order  in  LIST  format  (refer  to  Fig.  10  for 
nomenclature  conventions) : 

Item  Value 


1 Total  length  of  radcme  (in.); 

2 Radius  of  the  radome  base  (in.); 

3 Thickness  of  radome  wall  at  the  attachment 
ring  (in.)  ; 

4 Location  of  the  center  of  mass  of  the  whole 
missile,  measured  from  station  zero  (the  tip 
of  the  radome)  (in.); 

5 Density  of  the  radome  material  (lbm/it3)  ; 

6 Total  weight  of  missile  (lbm)  ; 

7 Maximum  permissible  lateral  acceleration 
(ft/s*) ; 

8 Maximum  permissible  rate  of  change  of  angle 
of  attack  (deg/s)  ; 

9 Maximum  permissible  acceleration  of  angle  of 
attack  (deg/s*) ; 

10  An  identification  number  for  the  trajectory 
table  used  for  this  flight  simulation  case 
(the  trajectory  table  identification  numbers 
are  assigned  when  the  time*dependent  data  are 
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X',Y'~  Coordinate  system  attached  to  ground 

X,Y  — Coordinate  system  along  and  perpendicular 

to  velocity  vector  of  missile  at  center  of  mass  of 


m m,  mr  are  the  mass  of  the  whole  missile  and  radome  only. 

RB  is  the  radome  base  radius 
CP  is  the  center  of  pressure  on  radome 
o is  the  missile  angle  of  attack 
g is  the  gravitational  acceleration 

F|\lr-  Fa,  are  the  normal  and  axial  forces  on  the  radome  due  to 
the  dynamic  pressure  of  the  air  stream 

and  represent  first  and  second  derivatives  with  respect 
to  time 

© Centroid  of  radome  or  missile 


Fig.  10  Coordinate  System  and  Forces  for  Radome  Mechanical  Load  Analysis 
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Item  Value 

read  in,  under  the  control  of  variable  CALL- 
TIM.) ; and 

11  A psuedo-node  number  for  obtaining  history 

plots  of  the  maximum  tensile  aerodynamic  load 
(if  this  value  is  zero  then  no  action  will 
take  place;  if  nonzero,  the  value  must  not 
correspond  to  an  actual  thermal  network  node 
number) . 

The  sample  problem  input  shows  the  30  LIST  format- 
ted cards  that  were  read  under  the  control  of  CALLLOD, 
which  was  set  to  30  earlier.  The  values  read  in  are 
reformatted  and  printed  on  SYSFRINT  in  the  sample  out- 
put. 


TIME-DEPENDENT  PARAMETERS 


If  the  control  variable  CALLTIM  is  nonzero  then  the 
next  data  in  the  SYSIN  file  will  pertain  to  the  reading 
in  of  tables  of  certain  variables  versus  time;  these 
variables  will  be  explained  subsequent  to  the  format 
description.  The  control  values  that  must  be  on  the 
SYSIN  file  are  (in  LIST  format): 

1.  A file  name  specifying  what  file  will  contain 
the  table  cf  data  and 

2.  The  number  of  tabular  data  values  to  be  read. 

On  the  file  named  in  1 (above)  there  should  be  card 
images  (80  character  records)  with  the  following  data:  2 
comment  cards  (160  characters)  followed  by  the  specified 
number  of  tabular  data  cards  in  the  following  EDIT 
format : 

Columns  Value 


11-9  Value  of  time  (s)  , F-format; 

10-19  First  dependent  variable,  F-format; 
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Columns  Value 

20-29  Second  dependent  variable,  F-format; 

30-39  Third  dependent  variable,  F-format; 

40-49  Fourth  dependent  variable,  F-format;  and 

50-59  Fifth  dependent  variable,  F-format. 

Each  of  the  CALL TIN  tables  will  be  assigned  a 
serial  identification  number  starting  with  1.  These 
identification  numbers  are  then  input  by  the  user  at 
other  points  in  the  data  for  the  various  routines  re- 
quiring the  particular  data.  For  aerodynamic  heating 
the  required  column  format  is: 

Columns  Required  Variable 

1-9  Time  variable  (s)  ; 

10-19  Altitude  values  (ft); 

20-29  Velocity  (ft/s) ; and 

30-39  Values  of  externally  applied  infrared  radia- 
tion to  the  aerodynamic  surface  (Btu/ft2-s) . 

In  0FLIM  runs  where  the  aerodynamic  loading  is 
calculated  (CALLLOD  * 0)  the  values  of  the  maximum 
permissible  angle  of  attack  must  also  be  provided;  e.g.: 

Columns  Required  Variable 

40-49  Maximum  permissible  angle  of  attack  (deg)  . 

In  the  case  where  boundary  conditions  of  tempera- 
ture or  heat  flux  histories  must  be  input  then  the 
specific  number  of  the  column  positions  is  not  critical 
but  the  specific  dependent  variable  number  will  be 
needed  by  other  parts  of  the  program.  Further  explana- 
tion of  this  aspect  of  the  time-dependent  data  is  given 
in  the  next  section,  where  input  for  the  temperature  and 
heat  flux  boundary  conditions  is  described. 

The  sample  problem  input  shows  10  trajectory  tables 
were  read  in  as  time- dependent  input  for  the  aerodynamic 
heating  calculations  (CALLTIM  was  set  to  20.)  The 
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control  cards  indicate  that  the  tables  are  on  the  SYSIN 
file  so  that  each  table  follows  immediately  each  control 
card.  (Note  that  the  dependent  variable  ordering  con- 
vention is  adhered  to  and  that  the  maximum  permissible 
angle  of  attack  is  set  in  all  cases  to  20°. ) The  sample 
output  shows  each  of  the  tables  printed  out  and  refor- 
matted slightly.  The  serial  identification  numbers  have 
also  been  assigned. 


HEAT  FLUX  BOUNDARY  CONDITIONS 


The  next  input  to  the  URLIH  program  is  for  extern 
ally  supplied  heat  fluxes  and  is  controlled  by  the 
variable  CALLFLX.  In  the  SYSIN  file  there  should  be 
CALLFLX  number  of  data  groups  with  the  following  items 
in  LIST  format: 

Item  Value 


1.  The  identification  number  of  the  time  table 
where  the  values  of  externally  specified  heat 
flux  (Btu/ft2-h)  are  fcund; 

2.  The  dependent  variable  identification  number 
within  the  above  named  table; 

3.  The  identification  number  for  the  thermal 
property  table  of  the  material  to  which  the  flux 
is  applied; 

4.  The  number  of  thermal  nodes  to  which  the  flux  is 
applied;  and 

5.  For  each  of  the  nodes  to  which  the  flux  is 
applied,  the  following  three  values; 

a.  the  thermal  node  number; 

b.  the  surface  area  (ft2)  ; and 

c.  the  view  factor  to  the  radiation  source. 

In  the  sample  problem,  there  are  no  heat  flux 
boundary  conditions  (i.e.,  CALLFLX  3 0)  and  so  there  is 
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no  input  to  this  section  in  the  sample  input.  In  the 
cases  where  heat  flux  data  is  input,  the  values  are 
reformatted  and  printed  on  SYSFRINT  for  reference. 


FORCED  TEMPERATURE  BOUNDARY  CONDITION 


If  the  control  variable  CALLFRC  is  set  to  a 
nonzero  value  then  that  number  of  sets  of  data  values 
(in  LIST  format)  must  appear  in  the  SYSIN  data  stream. 
These  values  are: 

1.  The  identification  number  of  the  time  table 
where  the  values  of  externally  specified  tem- 
peratures are  found; 

2.  The  dependent  variable  identification  number 
within  the  above  table; 

3.  The  number  of  thermal  nodes  that  will  be  forced 
to  follow  the  specified  temperature  history;  and 

4.  The  list  of  thermal  nodes 

Since  the  variable  CALLFRC  was  set  to  zero,  there  are  no 
data  in  the  sample  input.  In  the  cases  where  tempera- 
ture boundary  conditions  are  input,  the  values  are 
reformatted  and  printed  on  SYSFRINT  for  reference. 


THERMAL  STRESS  PARAMETERS 


If  thermal  stresses  are  to  be  calculated,  then  the 
control  variable  CALLSIG  should  be  set  to  indicate  the 
number  of  thermal  stress  cases  that  will  be  needed.  For 
each  of  the  CALLSIG  number  of  cases  the  following  data 
will  be  expected  on  the  SYSIN  file,  in  LIST  format  (for 
a definition  of  terms  refer  to  Fig.  11  and  Appendix  K) : 
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Notes: 

• Total  wall  thickness  (bnf  j hj ) is  divided  into  n arbitrarily  sized  regions, 
creating  n+1  interfaces. 

• Temperatures  aie  provided  at  m arbitrarily  spaced  radius  values  and  are  constant 
with  respect  to  the  axial  and  circumferential  directions. 

• Radial  displacements  from  an  initial  isothermal  state  occur  for  each  interface,  Uj 

• Mechanical  properties  are  assumed  constant  within  regions  at  values  that  are 
temperature  averaged  with  respect  to  the  radial  temperature  variance. 

• Each  region  may  have  different  mechnical  properties 

Fig.  11  Thermal  Stress  Geometry  and  Nomenclature 
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Item  Value 

f l T 

1.  A file  name  for  indicating  where  the  node  versus 
radius  data  is  stored; 

2.  The  number  of  positions  through  the  wall  at 
which  temperatures  will  be  specified; 

3.  The  number  of  subdivisions  through  the  wall  for 
the  finite  difference  solution  to  the  thermal 
stress  problem; 

4.  The  pressure  at  the  outside  surface  of  the 
cylinder  (psi) ; 

5.  The  pressure  at  the  inside  surface  of  the  cylin- 
der (psi)  ; 

6.  The  temperature  at  which  the  cylinder  will  be 
assumed  stress-free  (°F)  ; 

7.  The  angle  of  incidence  of  the  radar  mainbeam 
(deg)  ; 

8.  A flag  for  the  type  of  boundary  condition  for 
the  ends  of  the  cylinder  (1  for  fully  constrai- 
ned ends;  0 fcr  unconstrained  ends); 

9.  A pseudo-node  number  used  for  saving  the  maximum 
tensile  stress  values  versus  time  for  plotting 
in  another  section  of  the  program  (the  plotting 
is  controlled  by  the  variable  PLOTS,  whose  value 
is  set  in  the  ctrol  section;  the  value  input 
here  will  be  a thermal  "node"  number  for  the 
stress  value  and  therefore  must  not  be  an  actual 
thermal  node  number  in  the  thermal  network). 

After  these  values  have  been  read  then  the  file 
named  in  1 (above)  will  be  read  for  the  following  data: 

2 comment  cards  (160  characters)  followed  by  tabular 
data  with  one  of  the  following  two  types  of  entries: 

a.  In  columns  40  to  55  a stress  boundary  card  that 
gives  a radial  dimension  (in.)  for  the  interface 
between  two  adjacent  stress  regions,  F-format; 
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b.  A thermal  information  card  that  gives  the  fol- 
lowing : 

(1)  In  columns  1 to  15,  the  radial  position  of 
the  thermal  node,  (in.),  F-format; 

(2)  In  columns  16  to  25,  the  thermal  network 
node  number  whose  node  center  is  at  the 
corresponding  radial  position,  F-for*at;  and 

(3)  In  columns  26  to  35,  the  mechanical  property 
identification  number  (used  for  accessing 
the  table  of  temperature  dependent  mechani- 
cal properties) . 

These  data  (items  a or  b)  must  not  occur  on  the 
same  card.  Moreover,  the  separate  cards  are  to 
be  arranged  in  a manner  that  resembles  tne 
geometry  being  modeled;  that  is,  beginning  with 
the  innermost  radius  (the  inner  surface)  the 
data  cards  of  type  a or  b are  interleaved  in 
increasing  radial  order,  ending  with  the  outer 
surface.  Note  that  the  inner  and  outer  surface 
of  the  wall  section  being  modeled  are  cards  of 
type  a and  must  be  both  the  first  and  tne  last 
cards  of  the  sequence.  Thermal  node  and  inter- 
face cards  may  have  the  same  value  of  radius. 

The  EDlT-format  input  described  above  for  the 
radius-temperature  information  may  be  circumvented  if 
desired  by  an  alternative,  free-form  input  style  through 
the  setting  of  control  variable  SIGLIST.  If  SIGLIST  is 
set  to  1,  then  the  radius-temperatur e-node  information 
for  each  of  the  CALLSIG  cases  may  be  read  in  LIST  format 
from  the  file  named  above.  If  SIGLIST  is  set  to  1,  then 
the  following  data  will  be  read  from  the  file  named  in  1 
(above)  in  LIST  format  for  each  of  the  CALLSIG  cases: 

Item  Value 


1.  For  each  position  through  the  wall  at  which 

temperatures  are  known,  a triplet  (in  increasing 
radial  order)  of: 

a.  the  thermal  node  number, 

b.  the  mechanical  property  code  numbers,  and 
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Item  Value 

c.  the  radius  (in.). 

2.  For  each  radial  subdivision  interface: 
a.  the  radial  position  (in.)* 

The  sample  problem  has  the  control  variable  CALLSIG 
set  to  30  and  SIGLIST  set  to  0r  so  there  are  30  data 
groups  shown  in  the  sample  input  (Appendix  B) . In  all 
of  the  cases,  the  data  tables  of  nodal  positions  are  on 
the  file  SYSIN.  For  the  first  case,  the  control  data 
indicate  that  the  wall  will  have  three  nodal,  temper- 
atures specified  and  the  wall  will  be  divided  into  only 
two  subsections  for  the  stress  solution.  The  inner  and 
outer  surface  pressures  are  both  zero,  the  no-stress 
condition  occurs  at  59.3  °F,  the  radar  incidence  angle 
at  this  station  is  65°,  and  the  end  condition  is 
specified  as  being  unconstrained.  The  last  zero  in  the 
control  data  specification  is  a null  input  because  the 
plot  feature,  controlled  by  the  variable  PLOTS,  was 
turned  off  (PLOTS  = 0). 

The  data  tables  for  each  case  are  shown  in  the 
sample  input  and  they  follow  the  formats  as  described 
here  and  show  much  repitition  of  data.  It  must  be 
noted,  however,  that  each  thermal  node  number  is  dif- 
ferent throughout  all  this  input  so  each  case  is  unique. 

In  the  sample  output  there  is  a reprint  of  each  of 
the  geometry  tables  as  read  in,  with  sequential  case 
numbers  assigned. 


BORESIGHT  ERROR  SLOPES 


For  calculating  the  boresight  error  slopes  several 
parameters  must  be  input  for  each  case  and  they  should 
occur  at  the  next  position  in  the  SYSIN  file.  If  the 
control  variable  CALLBSE  is  nonzero  then  that  number  of 
data  groups  must  appear  in  the  SYSIN  file  in  the  follow- 
ing order,  LIST  format: 

1.  A case  number  specifying  which  thermal  stress 

case  is  to  be  used  for  this  boresight  error  case 
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(it  is  often  convenient  to  use  a one-to-one 
correspondence  between  thermal  stress  and 
boresight  error  cases  but  this  is  not  required) ; 

2.  The  radar  frequency  to  be  used  for  the  boresight 
calculation  (Hz)  ; 

3.  The  phase  center  separation  for  the  radar 
antenna  (in.)  ; and 

4.  A psuedo-node  number  for  obtaining  off-line 
CalCoap  plots  of  boresight  error  and  electrical 
thickness  change  histories  (if  this  value  is 
zero,  then  no  action  will  take  place;  if  this 
value  is  nonzero,  then  this  value  and  the  next 
sequential  value  will  be  used  to  store  boresight 
error  and  electrical  thickness  change,  respec- 
tively, at  each  time  point  indicated  for  these 
calculations) . 

In  the  sample  input  data  (directly  after  the  last  ther- 
mal stress  geometry  input)  there  are  30  triplets  of  data 
for  the  30  boresight  error  cases  to  be  made.  Tne  data 
show  the  one-to-one  correspondence  between  boresight 
error  case  number  and  thermal  stress  case  number.  Also 
each  calculation  is  done  at  a K-band  frequency  (30  GHz) 
with  the  same  antenna  aperture  (3.89  in.).  The  sample 
output  shows  the  data  reprinted  on  the  SYSPRINT  file  for 
reference. 


VELOCITY  LIMITS 


The  next  section  of  data  to  be  found  on  the  SYSIN 
file  concerns  the  specification  of  how  the  velocity 
limits  will  be  calculated.  In  this  section  the  selected 
parameters  that  will  be  potentially  limiting  are  com- 
bined with  the  appropriate  trajectories  and  design 
allowable  properties  so  that  the  URLIM  program  can 
generate  the  velocity  limits  plots.  For  each  plot  that 
is  called  for,  several  identifying  numbers  must  be  sup- 
plied in  a particular  order  so  as  to  fully  specify  the 
limit  plot.  These  numbers  are  to  be  on  the  SYSIN  file 
and  in  LIST  format: 

1.  The  number  of  velocity  versus  time  tables  (tra- 
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jectories)  that  are  to  be  plotted  and 

2.  The  number  of  capacitors  that  will  be  checked 
for  melting  on  each  of  the  trajectories. 

Beginning  at  the  first  column  of  the  next  logical 
record : 

3.  A title  card  of  80  characters  or  less, 

4.  The  thermal  property  material  code  number  of  the 
material  for  which  the  plct  is  being  made, 

5.  The  mechanical  property  identification  number  of 
the  material  for  which  the  plot  is  being  made, 

6.  For  each  of  the  number  of  trajectories  specified 
in  item  1 above,  the  following: 

a.  The  node  numbers  to  be  checked  for  melting 
(there  must  be  as  many  entries  here  as 
specified  in  item  2 above  for  each  trajec- 
tory) , 

b.  The  thermal  stress  case  number  that  will  be 
examined  for  excessive  thermal  stress, 

c.  The  aerodynamic  load  case  number  that  will 
be  examined  for  excessive  attachment  stres- 
ses, and 

d.  The  trajectory  identification  number  for 
this  trajectory  (i.e.,  the  number  of  the 
time  table  to  be  used  for  the  altitude  and 
velocity  versus  time  data)  . 

Since  the  variable  CA LLLI M in  the  sample  problem  was  set 
to  6,  there  are  six  groups  of  data  in  the  sample  input 
that  follow  the  above  format.  The  first  one  specifies 
that  there  are  to  be  five  trajectories  examined  and  that 
on  each  of  them  3 nodes  will  be  examined  for  melting. 

The  comment  card  tells  what  material  and  what  flight 
path  the  limit  plot  is  for.  The  thermal  and  mechanical 
material  identification  numbers  are  both  1 (and  should 
therefore  correspond  to  alumina) . The  nodes  that  will 
be  subject  to  melting  on  the  first  trajectory  are  nodes 
No.  1,  No.  2,  and  No.  3;  the  thermal  stresses,  aeroioad, 
and  trajectory  data  are  those  cases  numbered  1.  The 
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second  trajectory  has  melt  nodes  No.  4,  No.  5,  and  No.  6 
and  the  stress  and  load  data  are  to  come  from  those 
cases  numbered  2.  The  second  time  table  is  to  be  used 
for  velocity  versus  time  data.  The  rest  of  the  data  for 
the  first  limit  plot  follow  in  the  same  way  for  the 
other  trajectories.  The  other  five  sets  of  data  that 
specify  the  other  limit  plots  follow  immediately  in  the 
SYSIN  file.  All  of  this  bookkeeping  data  for  the  limits 
plots  is  reformatted  and  printed  on  the  output  file. 

The  sample  output  shows  the  six  resulting  tables,  which 
can  be  very  useful  later  when  the  analytical  data  are 
received. 


CALCOMP  PLOTS 


Following  the  limits  data  there  are  two  sections 
that  deal  with  the  analytical  output  of  the  program. 

The  first  of  these  deals  with  the  CalComp  plot  facility 
and  is  controlled  by  the  variable  PLOTS.  If  PLOTS  is 
set  to  any  nonzero  value,  then  immediately  after  the 
limits  data  on  the  SYSIN  file  the  following  items  should 
be  input  in  LIST  format: 

1.  A file  name  identifying  a file  onto  which  will 
be  written  the  information  required  to  make  the 
CalComp  plots.  There  must  be  a JCL  card  descri- 
bing this  file  included  with  the  execution  step 
of  the  UR L IH  run.  An  example  of  such  a card 
would  be: 

//filename  DD  DSN=group. userid. xxx, 

//  DIS P= (N EW , C ATLG) , 

//  UNIT=SA  VE,SPACE={3500,  (10,10)  , RLSE) 

The  plots  are  not  made  during  the  URLIH  run,  but 
the  required  data  are  saved.  After  the  URLIH 
program  execution,  when  the  printed  data  have 
been  inspected,  a job  may  be  run  to  plot  selec- 
ted variables  as  required  (c.f.  Appendix  P for 
details  of  the  execution  of  this  secondary 
program)  . 

2.  Two  40-character  strings  (each  enclosed  in 
single  quotes)  to  be  used  as  titles  in  the  plot- 
ting. 
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3.  The  number  of  times  at  which  the  data  points  are 
to  be  recorded. 

4.  The  list  of  times  at  which  the  temperature  data 
will  be  written  on  the  file  named  in  item  1 
(above)  (s)  . 

If  the  number  of  times  (item  3)  is  given  as  1 and 
the  single  value  input  (item  4)  is  negative,  then  the 
temperatures  will  be  written  out  to  the  file  (item  1)  at 
every  time  step  of  the  UELIM  program  execution.  If  the 
number  of  items  is  1 and  the  value  input  is  positive 
then  this  value  will  be  used  successively  as  an  incre- 
ment of  time  from  the  starting  time  at  which  the  values 
will  be  written  out. 

In  the  sample  problem  input  there  are  no  data  for 
this  option  because  the  control  variable  PLOTS  was  set 
to  zero.  When  data  are  input  they  are  echoed  directly 
to  the  SYSPRINT  file  for  reference. 


PRINTED  RESULTS 


During  the  execution  of  UELIM  the  various  analy- 
tical subroutines  can  be  called  upon  to  print  their 
results  at  intervals  specified  by  a list  of  times.  This 
facility  is  not  selectable;  i.e. , the  following  data 
items  must  be  next  in  the  SYSIN  data  stream  in  LIST 
format. 

1.  The  number  of  times  at  which  data  are  to  be 
printed  and 

2.  The  list  of  times  for  printout  (s) . 

Note:  If  the  value  for  item  1 is  1 and  the  single 
value  for  item  2 is  negative  then  data  will  be  printed 
at  every  time  step  during  the  execution  (use  with 
caution  --  considerable  output  can  be  made  this  way)  . If 
the  value  for  item  1 is  1 and  the  single  value  for  item 
2 is  positive,  then  the  item  2 value  will  be  used  as  an 
increment  in  determining  successive  times  from  the 
initial  time  at  which  output  will  be  printed.  If  the 
first  value  in  the  list  of  times  to  print  (item  2)  is 
999  then  the  format  for  the  printed  temperatures  will  be 
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changed  to  a table  that  will  include  the  node  numbers, 
the  temperatures,  the  total  net  conductances,  and  the 
net  heat  flows  for  each  node. 

The  type  of  printed  output  can  be  controlled  by  the 
control  variable  MINPFNT.  By  default  hINPRNT  has  the 
value  0 but  it  may  be  set  equal  to  the  values  1 or  2. 
These  values  have  the  following  effects: 

0 The  printed  output  is  directed  in  its  entirety 
to  the  SYSPRI NT  file.  The  individual  analytical 
routines  print,  in  the  order  of  their  execution, 
detailed  reports  of  the  calculations  involved. 

1 The  printed  output  is  directed  in  its  entirety 
to  the  SYSPRINT  file.  The  individual  routines, 
however,  print  only  a much  abbreviated  (one- 
line)  report  of  their  calculations. 

2 The  analytical  routines  print  their  abbreviated 
output  on  the  file  SYSPRINT.  The  analytical 
routines  also  print  their  extended  output 
reports  to  the  file  MAXPFNT,  which  must  have  a 
JCL  card  describing  its  format  for  the  operating 
system. 

The  sample  problem  was  run  with  MINPFNT  set  to  2; 
the  JCL  statements  in  Fig.  6 will  show  how  the  extensive 
output  data  file  MAXPFNT  was  described.  There  is  a 
facility  available  at  APL  for  obtaining  microfiche 
output  from  print  files  like  MAXPFNT.  This  is  accom- 
plished by  specifying  "SYSOUT  = Q"  on  the  DD  card. 
Appendix  C contains  the  results  of  the  minimum  printing 
and  Appendix  D shows  a sample  of  the  extensive  output 
that  was  written  on  microfiche. 

There  is  a significance  to  the  list  of  printout 
times  in  addition  to  that  described  so  far:  the  list  of 
times  indicates  how  frequently  the  analytical  routines 
for  the  limits  calculations  are  executed.  Specifically 
these  routines  are  the  ones  that  deal  with  thermal 
stress,  boresight  error,  and  aerodynamic  load.  Some 
consideration  must  therefore  be  given  to  how  frequently 
these  parameters  will  be  evaluated  during  the  flight 
simulation  so  that  the  respective  parameter  histories 
are  not  too  coarse. 

In  the  sample  input  27  time  points  are  specified 
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for  printout  ranging  from  0 to  45  s.  There  is  no  direct 
indication  in  the  sample  output  that  these  values  were 
read;  however,  the  output  from  the  analytical  routines 
is,  of  course,  printed  at  the  times  that  were  specified. 


FINAL  INPUT 


The  last  items  that  may  be  input  prior  to  the 
flight  simulation  include  the  values  listed  in  Table  1. 
Any  other  variable  used  in  the  URLIM  main  program  (of 
which  there  are  many)  may  be  given  values  at  this  point 
in  the  SYSIN  file  by  using  the  DATA  format.  These 
variables  will  have  had  values  assigned  to  them  during 
the  earlier  input  but  they  may  have  their  values  changed 
in  this  section.  This  feature  is  mainly  for  the  user 
with  intimate  knowlege  of  the  URLIM  code  and  is  not 
generally  usefull.  However,  whether  or  not  any  varia- 
bles have  their  values  changed,  this  section  must  end 
with  a semicolon. 

In  the  sample  input  there  were  no  DATA-format 
entries  so  only  a semicolon  appears.  The  semicolon  also 
appears  on  the  sample  output  as  the  DATA-format  entries 
would,  if  there  were  any. 

At  this  point  the  input  is  complete  and  the  program 
will  execute  according  to  the  various  data  it  has  been 
given. 
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4.  PhOGRAM  OUTPUT 


In  addition  to  the  system -generated  messages  (JCL, 
data  set  allocations,  and  disposition)  generic  to  the 
particular  operating  system,  the  output  received  from  ar. 
URLIM  run  will  be  in  two  general  parts.  The  first  part 
will  echo  some  of  the  input  data  so  that  it  can  be 
verified  and  later  referred  to  by  the  user.  The  second 
part  contains  the  results  of  the  calculations  by  the 
various  analytical  subroutines  and  is  printed  as  the 
time  variable  in  the  program  advances.  This  second 
section  of  the  output  can  be  thought  of  as  having  two 
parts:  the  first  is  the  normal  output  of  the  various 
program  subroutines  and  the  other  part  is  a set  of 
warning  or  error  messages  that  may  occur  during  execu- 
tion. Appendix  E lists  these  warning  messages.  Appen- 
dix C contains  a part  of  the  output  from  the  sample 
problem  that  was  printed  on  SYSPRINT  and  Appendix  D 
contains  a sample  of  the  output  that  was  written  on 
microfiche  from  the  file  MAXPRNT.  These  sample  results 
will  be  used  to  illustrate  the  general  discussion  that 
will  follow  the  program's  output. 


SECTION  I 


This  output  is  a partial  echo  of  the  input  data  and 
has  been  described  in  detail  in  Section  3 of  this 
report.  The  following  is  a summary  list  of  the  types  of 
data  that  may  be  printed  in  this  section: 

1.  A direct  echo  of  the  control  variable  inputs; 

2.  The  thermal  property  tables  of  the  program 
materials  (if  THPEOPS  > 0); 

3.  The  mechanical  property  tables  (if  any)  of  the 
program  materials,  (NCEROPS  > 0); 

4.  The  manually  input  thermal  network  information 
(if  MANLIN  > 0)  ; 
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5.  The  cylindrical  network  data,  if  any  (CALLCYL 
> 0)  ; 

6.  The  2-dimensional  Cartesian  network  data,  if  any 
(CALLTOD  > 0)  ; 

7.  The  air  property  data,  if  used  (CA HARO  > 0)  ; 

8.  The  altitude  tables,  if  used  (CALLARO  > 0)  ; 

9.  The  local  flow  ratio  tables  if  used  (CALLFLO 
> 0)  ; 

10.  The  missile  and  radome  force  coefficients,  if 
used  (CALLAFC  > 0); 

11.  The  time-dependent  tables,  if  any  (CALLT1M  > 0)  ; 

12.  The  thermal  stress  geometry  definitions,  it  any 
(CALLSIG  > 0)  ; 

13.  The  limits  directory  tables,  if  any  (CA1LLIM 
> 0)  . 


SECTION  II 


The  analytical  output  is  printed  by  each  subrou- 
tine, in  the  order  of  execution,  at  every  time  interval 
specified  by  the  input  data  for  printing.  The  sample 
output  for  Section  II  is  shown  in  Appendix  C and  is  the 
4-s  output  from  the  sample  problem.  When  the  output  for 
any  one  time  interval  is  complete,  a line  of  asterisks 
is  printed  as  a visual  aid  to  separating  the  outputs  of 
different  times.  The  output  information  may  contain  the 
following  data: 

1.  The  current  time  in  the  program  (s)  (always  on 
SYSPRINT) . 

2.  The  thermal  capacitor  temperatures  (°F)  . Only 
the  thermal  nodes  whose  temperatures  are  dif- 
ferent from  the  initial  temperature  are 
printed.  Also,  as  a group  of  thermal  nodes 
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reaches  either  a liait  or  the  end  of  rheir 
trajectory,  they  will  be  dropped  froa  the  cal- 
culations. An  appropriate  warning  nessage  is 
printed  when  this  happens  (this  output  is 
always  on  SYSPBINT). 

3.  If  aerodynaaic  heating  is  being  used  (CALLABO 
> 0)  then  the  subroutine  that  calculates  the 
local  flow  conditions  and  the  surface  tempera- 
ture  will  print  its  data  (according  to  the 
value  of  MIN  PENT)  next.  The  aininua  output, 
which  is  always  on  SYSPBINT,  is  only  one  line 
and  specifies:  (1)  the  case  number,  (2)  the 
Mach  number,  and  (3)  the  iterations  reguired  to 
reach  the  solution  to  the  surface  energy 
balance.  Appendix  C shows  this  output  tor  the 
sample  problem.  If  the  maximum  output  is 
called  for,  then  more  inforaation  will  be 
printed.  As  shown  in  Appendix  D,  the  extended 
output  from  the  surface  aerodynamic  calcula- 
tions includes:  the  block  identification  num- 
ber; the  freestreaa  and  local  conditions  of 
temperature,  pressure,  velocity,  Mach  number, 
and  enthalpy;  and  data  pertinent  to  the  local 
flow  calculation  accuracy.  For  each  surface 
node  in  the  aero  block  the  following  is  prin- 
ted: the  node  number  and  its  newly  calculated 
temperature;  the  enthalpy-based  heat  transfer 
coefficient;  the  net  convective  heat  flux  to 
the  hot  surface  (Btu/ft2-s)  ; the  Reynolds 
number;  the  reference  temperature,  the  adia- 
batic wall  teaperature  and  enthalpy;  the  heat 
to  and  from  the  wall  via  radiation  (Btu/ft*-s); 
the  heat  conducted  from  the  surface  into  the 
structure  (Btu/ft*-s)  ; the  enthalpy  of  the 
stream  at  the  wall  temperature;  the  convergence 
parameters  indicating  the  accuracy  of  the  cal- 
culations; and  the  compressibility  of  the 
boundary  layer  gas  at  the  reference  condition. 

4.  If  aerodynaaic  loads  are  being  calculated 
(CALLLOD  > 0)  then  the  lead  routine  will  print 
its  values  at  the  current  time.  Appendix  C 
shows  samples  of  this  output  in  abbreviated 
fora  with  (1)  the  value  of  the  maximum  stress 
in  the  attachment  area  (2)  the  case  number,  and 
(3)  the  angle  of  attack  of  the  aissile.  (If 
the  angle  of  attack  is  followed  by  an  asterisk 
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then  the  printed  angle  of  attack  is  the  maximum 
peraissible  and  the  aissile  is  not  capable  of 
executing  the  required  maneuver;  otherwise  the 
angle  of  attack  is  that  which  is  required  to 
pull  the  full  maneuver  acceleration.)  Appendix 
D shows  some  examples  of  the  extended  output 
for  the  aero-load  routine.  In  this  output  the 
case  number  is  identified  and  the  time  is 
printed,  the  net  normal  and  axial  forces  on  the 
radome  are  printed  (lbf) , the  net  bending 
aoaent  at  the  attachaent  is  listed  (ft-lbf)  , 
the  windward  and  leeward  stresses  at  the  at- 
tachment point  are  shown  (psi),  the  maxinum 
allowable  and  the  current  angles  of  attack  are 
given  (deg)  , and  the  altitude  (ft)  and  Mach 
number  are  listed. 

5.  If  thermal  stresses  are  calculated  (CALLSIG 

> 0)  then  the  thermal  stress  routine  will  print 
its  abbreviated  output  (if  required)  in  the 
format  shown  in  Appendix  C.  This  information 
is  simply  the  case  number  and  the  maximum 
tensile  stress  in  the  wall.  The  extended 
output,  if  requested,  contains  a full  Descrip- 
tion of  the  thermal  stress  analysis  ana  is 
shown  in  the  example  in  Appendix  D.  The  output 
identifies  the  case  number  and  the  current  time 
of  the  trajectory  simulation.  The  maximum 
tensile  stress  is  also  given.  A table  is  then 
printed  that  shows  the  stresses  and  temperat- 
ures as  a function  of  radius  along  with  the 
effective  elastic  modulus  and  the  net  radial 
expansion.  Thermal  node  numbers  (under  the 
"INDEX"  column)  are  also  shown  at  their  respec- 
tive radii.  The  other  temperatures  given  at 
other  radii  are  interpolated  temperatures.  The 
hoop  and  axial  stresses  are  printed  at  each 
external  surface  and  at  either  side  of  each 
internal  stress  boundary.  These  stresses  may 
vary  at  each  interface  because  of  the  variance 
of  the  elastic  modulus  between  adjacent  re- 
gions. The  radial  stresses  at  the  surface  are 
siaply  the  externally  applied  pressures. 

Within  the  wall  the  radial  stresses  exist 
across  each  boundary.  The  last  line  from  the 
thermal  stress  routine  gives  the  average  axial 
expansion  (in. /in.)  and  also  gives  the  net 
electrical  thickness  change  in  terms  of  the 
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insertion  phase  delay  through  the  wall  thick- 
ness. 

6.  If  boresight  error  slopes  are  being  calculated 
during  the  0 KLIM  run,  they  will  be  printed  on 
the  SYSPRINT  or  HAXPRNT  output  files  (as  re- 
quired) in  the  same  format;  i.e.,  the  total 
output  is  simply  the  foresight  error  case 
number  and  the  boresight  error  slope  (degrees 
of  error  per  degree  of  lcok  angle). 

7.  If  flight  limits  are  being  monitored  and  if  a 
flight  limit  has  been  reached,  a one-line  mes- 
sage is  printed  out  to  the  SYSPRINT  file 
indicating  what  limit  was  exceeded  and  which 
limit  case  number  is  involved.  Also  the  exact 
case  number  of  the  particular  limiting  parame- 
ter is  listed  along  with  the  time  of  failure. 

8.  When  all  the  output  from  a single  time  has  been 
printed,  a statement  regarding  the  time  step 
and  the  thertal  node  that  is  limiting  the  value 
of  the  time  step  is  printed.  The  node  is  the 
"CRITICAL  INDEX  NO."  and  can  be  reviewed  by  the 
user  to  determine  if  thermal  model  changes  need 
to  be  made  to  raise  the  time  step. 

9.  Nhen  the  URLIH  program  has  finished  the  flight 
simulation  aspect  of  its  execution  and  if 
limits  were  monitored,  then  the  velocity  limits 
plots  will  be  printed  out  on  the  SYSPRINT  file. 

10.  After  the  limits  plots,  if  any,  a final  message 
will  be  printed  indicating  the  final  time  in 
the  simulation  and  the  total  number  of  time 
steps  taken  during  the  simulation.  This  mes- 
sage will  be  the  last  output  from  a successful 
execution  of  the  URLIH  program. 

The  only  other  output  will  be  operating  system  messages; 
the  URLIH  output  is  completed. 

Hention  has  been  made  throughout  this  report  of  a 
sample  problem  and  a library  of  subroutines  used  by 
URLIH  in  its  calculations.  The  appendixes  that  follow 
contain  the  detailed  descriptions  of  these  parts  of  the 
URLIH  code. 
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Appendix  A 
DATA  INPUT  FORMATS 


Data  that  are  prepared  for  reading  into  the  URLIM 
program  will  be  character  or  stream  type  data.  Most 
sections  of  URLIM  will  expect  the  input  data  files  to 
have  fixed  length  record  formats  (RECFM)  of  80  charac- 
ters (bytes)  each.  That  is,  the  Data  Control  Block 
(DCB)  of  the  Data  Definitions  (DD)  statement  should  be 

//filename  DD  DCB=  (RECFM=FB,LRECL=80,BLKSIZE=x*80) 

where  the  value  x*80  represents  a decimal  integer  that 
is  some  multiple  of  80.  The  data  on  these  cards  will  be 
in  one  of  three  basic  formats  according  to  what  the 
program  reguires.  These  three  formats  are: 

DATA  wherein  specific  program  variables  are  assig- 
ned values  by  name  in  the  data  stream;  e.g.: 

(variable  name)  = value 

this  format  is  free  (i.e.,  is  not  dependent  on 
column  position  on  the  card) . 

LIST  wherein  the  data  values  are  written  in  any 

decimal  format  in  the  specific  order  reques- 
ted. This  format  is  also  column  free  out  the 
values  must  be  in  the  precise  order  required. 

EDIT  wherein  the  data  values  must  be  in  specified 
formats  in  specified  column  fields  on  the 
cards.  In  EDIT  directed  input  it  will  be 
required  to  use  1-format,  F-format,  or  E- 
format  numerical  representations;  that  is: 

F-format  is  fixed  decimal  point  and  requires 
a decimal  point  to  be  used  (e.g.,  -49.662); 

I-format  where  integers  are  to  be  read 
(e.g.,  -287); 

E-format  where  a floating  point  (scientific 
notation)  number  is  used  (e.g.,  -836.49E-6). 
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In  this  example  the  E -6  represents  a power 
of  10  so  that  the  number  would  be: 

-836.49  x 10-*  = -0.00083649. 

In  all  of  the  above  examples  if  a sign  is  not  included, 
positive  values  will  be  assumed.  Also,  in  the  E-foraat 
there  may  be  no  spaces  between  ♦ or  - signs  and  the 
exponents  (E)  . 
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Appendix  B 
SAMPLE  INPOT 


The  following  JCL  and  data  cards  are  similar  to  the 
ones  used  to  execute  the  saaple  problem  referred  to 
throughout  this  report.  A few  of  the  salient  features 
include: 


• The  inclusion  of  a "//FLIDUMP  DD  " card  so  that 
the  storage  report  nay  be  printed  thereon. 

• The  SYSIN  file  is  a concatenated  file;  that  is, 
it  is  made  up  of  two  data  sets,  the  first  being  BBE. 
FRAZER. MHNDATA  (SAHPLXN)  and  the  other  a brief  card  file. 
The  contents  of  the  BEE.FRAZER. MHNDATA  (SAMPLIN)  dataset 
are  listed  directly  after  the  JCL  (below). 

• The  MAXPRNT  file  is  a regular  print  file 
(SYSOUT=  A)  but  could  just  as  easily  have  been  microfiche 
(SISOT=Q) . 


Contents  of  Dataset 
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CALLARO  = 30 
CALLLOD  = 30 
CALLSIG  = 30 
CALLBSE  = 30 
CA I ICY L = 30 
CALLFLO  = 1 
CALLLIfl  = 6 
CALLTIH  = 10 
BA  N LI  N = 1 

HCPROPS  = 3 HA  X E NTS  = 18 

TH  PROPS  = 3 

XC A PL IH  = 185 

•NODES  = 185 

AFC  FILE  = • FORCOEF ' 

TINIT  = 519 

•ALPHA  = i *HACH  = 1U 
NETPILE  = 'S  YSI N • 

HINSTEP  = .019 
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3600. 

41.1 

.86 

.690 

4.40 

SY3T  N • 18 

.010 

THERMAL  PROPERTIES 

PYROCERAM  - TG 

865 

TEMP  DEG  R 

RHO  CP 

CONDUCT* 

Y EHISS 

DIEL  CONST 

460. 

26.7 

2.22 

.850 

5.65 

560. 

30.8 

2.  16 

.846 

5.65 

660. 

33.4 

2.  10 

.842 

5.67 

760. 

35.6 

2.05 

.8  37 

5.67 

560. 

37.  3 

2.00 

.832 

5.68 

460. 

38.8 

1.96 

.822 

5.69 

1060. 

39.9 

1.93 

.810 

5.70 

1260. 

42.1 

1.87 

.765 

5.71 

1460. 

43.  7 

1.82 

.718 

5.75 

1660. 

45.4 

1.79 

.698 

5.79 

1860. 

47.0 

1. 76 

.678 

5.82 

20  60. 

48.6 

1.74 

.654 

5.88 

22  60. 

50.2 

1.73 

.630 

5.91 

2460. 

51.8 

1.72 

.617 

6.01 

2660. 

53.5 

1.71 

.595 

6.10 

2860. 

55.9 

1.70 

.5  87 

6.  13 

2°20. 

56.7 

1.70 

.580 

6.15 

24  30. 

56.7 

1.70 

.580 

6.  15 

•SYSIN'  11 

1 8700 

MECHANICAL 

PROPERTIES  - 

ALUMINA  TG  865 

TEMP  DEG  F 

MODULUS 

POIS 

EXPANSION 

0. 

50. CE6 

.28 

0.0  EO 

400. 

49.5E6 

.28 

1.3E-3 

800. 

49.CE6 

.28 

2.9E-3 

1200. 

47.8F6 

.28 

4. 8E- 3 

1600. 

46. CE6 

.28 

7.0F.-3 

2000. 

43.0E6 

.29 

9.2E-3 

2400. 

38. CF6 

.30 

11.4E-3 

2800. 

30. 0E6 

.30 

13.8E-3 

3200. 

18.CE6 

.30 

16.2E-3 

3600. 

0. 0E0 

.30 

1R.6E-3 

3610. 

0.  CEO 

.30 

18.7E-3 

•SYSIN'  11 

4000.0 

MECHANICAL 

PROPERTIES  FUSED  SILICA  - TG  865 

TEMP  DEG  F 

HCDfJLHS 

POIS 

EXPANSION 

0. 

8.  CF6 

.15 

0.0  EO 

400. 

9.0E6 

.16 

0. 18E-3 

800. 

1C.CE6 

.17 

0.42E-3 

1000. 

10.0E6 

.20 

0.46E-3 

1200. 

1C.  1E6 

.18 

0.41E-3 

1600. 

9.9E6 

.19 

0.48E-3 

2000. 

8.5P6 

.195 

0.60E-3 

24  00. 

6. 0E6 

.20 

0. 7 IE-3 

2800. 

1 . 5E6 

.20 

0.80F-3 

3140. 

0.0E0 

.20 

0. 8 IE-3 

3600. 

C.CEO 

.20 

1. OOF-3 

•SYSIN'  IP 

22500.0 

MECHANICAL  PROPERTIES  PYRCCERAH 

TG  865 

TEMP  DEG  F 

MODULUS 

POIS 

EXPANSION 

0.0 

15.CE6 

.25 

0.0  EO 

100. 

16. 5E6 

.25 

0.1E-3 

200. 

16.5E6 

.25 

. 5E-3 

400. 

16.  5F.6 

.25 

1.  IE-  3 

600. 

16.6E6 

.25 

1.6E-3 

800. 

16.7E6 

.25 

2.0E-3 

1000. 

16. 9E6 

.25 

2.3E-3 
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1200. 

17. CF6 

.25 

2.8E-3 

1400. 

16. 9F6 

.25 

3.3E-3 

1600. 

16. 6E6 

.25 

3.9E-3 

1000. 

15.0F6 

.25 

4.4E-3 

2000. 

14.3E6 

.25 

5.0E-3 

2100. 

13.0E6 

.25 

5.3E-3 

2200. 

11 . 5E6 

.25 

5.6E-3 

2)00. 

8.  5E6 

.25 

5.8E-3 

2400. 

5.5E6 

.25 

6.2E-3 

2462. 

0.  0E6 

. 25 

6.3E-3 

2470. 

C . CF6 

.25 

6.4E-3 

1 

1 

.07102 

2.216  1 

2 

0 0 1 

c 

5.4667E-3  2 360  1 

1 

1 

.07182 

2.216  1 

5 

0 0 1 

0 

5.4667S-3  2 360  1 

1 

1 

.07182 

2.2  16  1 

8 

0 0 1 

0 

5.4667E-3  2 360  1 

1 

1 

.07182 

2.216  1 

11 

C 0 1 

0 

5.4667E-3  2 360  1 

1 

1 

.07182 

2.216  1 

14 

0 0 1 

0 

5.46  67E- 3 2 360  1 

1 

1 

.07182 

2.216  2 

22 

0 0 1 

0 

9.2167E-3  6 

360 

1 

1 

1 

.07182 

2. 216  2 

29 

0 0 1 

0 

9.2167E-3  6 

360 

1 

1 

1 

.07182 

2.216  2 

36 

0 0 1 

0 

9.2167E-3  6 

360 

1 

1 

1 

.07182 

2.2  16  2 

43 

0 0 1 

0 

9.2167E- 3 6 

360 

1 

1 

1 

.07182 

2.216  2 

50 

0 0 1 

0 

9.2167E-3  6 

360 

1 

1 

1 

.07182 

2.216  3 

62 

0 0 1 

0 

7.04  1 67F-3  4 

360 

1 

1 

1 

.07182 

2.2  16  3 

67 

0 0 1 

0 

7 . 04  1 67  E- 3 4 

360 

1 

1 

1 

.07182 

2.216  3 

72 

0 C 1 

0 

7.04 167E-3  4 

360 

1 

1 

1 

.07182 

2.216  3 

77 

0 0 1 

0 

7. 04  167  E-3  4 

360 

1 

1 

1 

.07182 

2.216  3 

82 

0 0 1 

0 

7.04167E-3  4 

360 

1 

1 

1 

.07182 

2.216  1 

102 

0 0 1 

0 

5.4667E-3  2 360  1 

1 

1 

.07182 

2.216  1 

105 

0 0 1 

0 

5.4667E-3  2 360  1 

1 

1 

.07182 

2.216  1 

108 

0 0 1 

0 

5.46  67E-3  2 360  1 

1 

1 

.07182 

2.216  1 

111 

0 0 1 

0 

5.4667E-3  2 360  1 

1 

1 

. 07182 

2.216  1 

114 

0 0 1 

0 

5.4667E-3  2 360  1 

1 

1 

.07182 

2.216  2 

122 

0 0 1 

0 

9.2167E-3  6 

360 

1 

1 

1 

.07182 

2.216  2 

129 

0 0 1 

0 

9.2167E-3  6 

360 

1 

1 

1 

.07102 

2.216  2 

136 

0 0 1 

0 

9.2167E-3  6 

360 

1 

1 

1 

.07182 

2.216  2 

143 

0 0 1 

0 

9.2167F.-3  6 

360 

1 

1 

1 

.07182 

2.2  16  2 

150 

0 0 1 

0 

9.2167E-3  6 

360 

1 

1 

1 

.07182 

2.216  3 

162 

0 0 1 

0 

7 . 04  1 67  F- 3 4 

360 

1 

1 

1 

.07182 

2.216  3 

167 

0 0 1 

0 

7. 04  1 67 E-  3 4 

360 

1 

1 

1 

.07132 

2.216  3 

172 

0 0 1 

0 

7.04  167  E-3  4 

360 

1 

1 

1 

.07182 

2.216  3 

177 

0 0 1 

0 

7.04  167  E- 3 4 

360 

1 

1 

1 

.07182 

2.216  3 

182 

0 0 1 

0 

7 . 04  1 67  E ~3  4 

360 

1 

GFCMETFY  DEFINITIONS  FOR 

SURFACE  NODES 

MODE 

VOLUME 

MATE  • 

L 

CONS 

FIRST  A/L 

1 

0.0  E0 

1. 

1. 

2.  743. OFO 

4 

0.0E0 

1. 

1. 

5.  743. OEO 

7 

0.CE0 

1. 

1. 

8.  743. OEO 

10 

0.  0EC 

1. 

1. 

11.  743. OEO 

13 

0.0E0 

1. 

1. 

14.  743. OEO 

21 

0.  OEO 

2. 

1. 

22.  1294. OEO 

23 

0.  CEO 

2. 

1. 

29.  1294. OEO 

35 

C.C'O 

2. 

1. 

36.  1294. OEO 

42 

O.OEC 

2. 

1. 

43.  1294. OEO 

49 

0.0E0 

2. 

1. 

50.  1294. OEO 

61 

O.OEO 

3. 

1. 

62.  1135. 8E0 

66 

0.0E0 

3. 

1. 

67.  1135. 8E0 

71 

O.OEO 

3. 

1. 

72.  1135. 8E0 

76 

O.OEC 

3. 

1. 

77.  1135. 8E0 

81 

0.  CEO 

3. 

1. 

82.  1135. 8E0 

101 

O.OEO 

1. 

1. 

102.  743. OEO 

104 

O.OEO 

1. 

1. 

105.  743. OEO 

. 
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107 

0.0E0 

1. 

1. 

108. 

743. GFO 

110 

0.  CRO 

1 . 

1. 

111. 

743. OEO 

1 13 

0.0E0 

1. 

1. 

114. 

743. OEO 

121 

0.  02C 

2. 

1. 

122. 

1294. CFO 

129 

0.0E0 

2 • 

1. 

129. 

1294. OEO 

135 

0.0E0 

2. 

1. 

136. 

1294. OEO 

142 

0.0E0 

2. 

1. 

14  3. 

1294. OEO 

149 

O.OEC 

2. 

1. 

150. 

1294 . 0 EO 

161 

O.CEO 

3. 

1. 

162. 

1135. 8E0 

166 

O.CFC 

3. 

1. 

167. 

1135. 8E0 

171 

O.CEO 

3. 

1. 

172. 

1135. 8E0 

176 

0.  OF  0 

3. 

1. 

177. 

1135. 8 EO 

181 

0.  C EO 

3. 

1. 

162. 

1135. REO 

•AIR' 

7 28 

•ATflOS'  30 

1 1 

1 1 

1.0E6  1 

2 

1 10 

1 

2 

1 

. 13333 

1 

1 1 

2 1 

1.0E6  1 

2 

1 10 

4 

5 

1 

.13333 

1 

1 1 

3 1 

I.OFfe  1 

2 

1 10 

7 

8 

1 

. 1 3333 

1 

1 1 

4 1 

1 . OF  6 1 

2 

1 10 

10 

11 

1 

.13333 

1 

1 1 

5 1 

1.0E6  1 

2 

1 10 

13 

14 

1 

. 1 3333 

1 

1 1 

1 1 

1.0E6  1 

2 

1 10 

21 

22 

1 

.13333 

1 

1 1 

2 1 

1 . 0 E6  1 

2 

1 10 

28 

29 

1 

.13333 

1 

1 1 

3 1 

I.OEo  1 

2 

1 10 

35 

36 

1 

.13333 

1 

1 1 

4 1 

1.0E6  1 

2 

1 10 

42 

43 

1 

.13333 

1 

1 1 

S 1 

1 . OF  6 1 

2 

1 10 

49 

50 

1 

.13333 

1 

1 1 

1 1 

I.OEb  1 

2 

1 10 

61 

62 

1 

. 1 33  33 

1 

1 1 

2 1 

1.0F6  1 

2 

1 10 

66 

67 

1 

.13333 

1 

1 1 

3 1 

1.0E6  1 

2 

1 10 

71 

72 

1 

.13333 

1 

1 1 

4 1 

1 . OF  6 1 

2 

1 10 

76 

77 

1 

.13333 

1 

1 1 

5 1 

1.0E6  1 

2 

1 10 

81 

82 

1 

.13333 

1 

1 1 

6 1 

1.0F6  1 

2 

1 10 

101 

102 

1 

.13333 

1 

1 1 

7 1 

1.0E6  1 

2 

1 10 

104 

105 

1 

.13333 

1 

1 1 

8 1 

1.0E6  1 

2 

1 10 

107 

108 

1 

.13333 

1 

1 1 

9 1 

1.0F6  1 

2 

1 10 

110 

111 

1 

.13333 

1 

1 1 

10  1 

I.OEb  1 

2 

1 10 

113 

114 

1 

.13333 

1 

1 1 

6 1 

1.0E6  1 

2 

1 10 

121 

122 

1 

.13333 

1 

1 1 

7 1 

1 . OR  6 1 

2 

1 10 

128 

129 

1 

.13333 

1 

1 1 

8 1 

1.0E6  1 

2 

1 10 

135 

136 

1 

.13333 

1 

1 1 

9 1 

1.0E6  1 

2 

1 10 

142 

143 

1 

.13333 

1 

1 1 

10  1 

1.0F6  1 

2 

1 10 

149 

150 

1 

.13333 

1 

1 1 

6 1 

1 . OE  6 1 

2 

1 10 

161 

162 

1 

.13333 

1 

1 1 

7 1 

1.0F6  1 

2 

1 10 

166 

167 

1 

.13333 

1 

1 1 

8 1 

1.CE6  1 

2 

1 10 

171 

172 

1 

.13333 

1 

1 1 

9 1 

1.0E6  1 

2 

1 10 

176 

177 

1 

.13333 

1 

1 1 

10  1 

1 . OE  6 1 

2 

1 10 

181 

182 

1 

.13333 

1 

•FLORATS*  17 

28.1 

6.75 

.0656 

110 

248 

1500 

966 

0 

650 

1 

0 

28.3 

6.75 

.0656 

110 

248 

1500 

966 

0 

650 

2 

0 

28.3 

6.75 

. 0656 

110 

248 

1500 

966 

0 

650 

3 

0 

28.3 

6.75 

. 0656 

110 

248 

1500 

966 

0 

650 

4 

0 

28.3 

6.75 

. 0656 

110 

248 

1500 

966 

0 

650 

5 

0 

28.3 

6.75 

. 1 106 

110 

120 

1500 

966 

0 

650 

1 

0 

28.3 

6.75 

. 1106 

110 

120 

1500 

966 

0 

650 

2 

0 

28.3 

6.75 

. 1106 

110 

120 

1500 

966 

0 

650 

3 

0 

28.  3 

6.75 

. 1106 

110 

120 

1500 

96  6 

0 

650 

4 

0 

28.3 

6.75 

. 1106 

110 

120 

1500 

966 

c 

650 

5 

0 

28.3 

6.75 

.oeus 

110 

162 

1500 

96  6 

0 

650 

1 

0 

28.3 

6.75 

. 0845 

110 

162 

1500 

966 

0 

650 

2 

0 

28.3 

6.75 

. 0845 

110 

162 

1500 

96  6 

0 

650 

3 

0 
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28.3 

6.75 

.0845  110 

162 

1500 

966 

0 

650 

4 

0 

28.3 

6.75 

.0845  110 

162 

150  0 

966 

0 

650 

5 

0 

28.  3 

6.75 

.0656  110 

248 

1500 

966 

0 

650 

6 

0 

28.3 

6.75 

.0656  110 

248 

1500 

966 

C 

650 

7 

0 

28.3 

6.75 

.0656  110 

248 

1500 

96  6 

0 

650 

8 

0 

28.3 

6.75 

.0656  110 

248 

1500 

966 

0 

65C 

9 

0 

28.3 

6.75 

.0656  110 

240 

1500 

966 

0 

650 

10 

0 

28.  3 

6.75 

.1106  110 

120 

1500 

966 

0 

650 

6 

0 

28.3 

6.75 

.1106  110 

120 

1500 

96  6 

0 

650 

7 

0 

28.3 

6.75 

.1106  110 

120 

1500 

96  6 

0 

650 

8 

0 

28.3 

6.75 

.1106  110 

120 

1500 

96  6 

0 

650 

9 

0 

28.3 

6.75 

.1106  110 

120 

1500 

96  6 

C 

650 

10 

0 

28.3 

6.75 

.0845  110 

162 

1500 

966 

0 

650 

6 

0 

28.3 

6.75 

.0845  110 

162 

1500 

96  6 

0 

650 

7 

0 

28.3 

6.75 

.0845  110 

162 

1500 

966 

0 

650 

8 

0 

28.3 

6.75 

.0845  110 

162 

1500 

966 

0 

65  C 

0 

0 

28.3 

6.75 

.0845  110 

162 

1500 

966 

0 

650 

10 

0 

• SY  SIN ' 

2 

TIME  TABLE  FOB  VI  TRAJ 

a qe  = 

20 

TIME 

ALT 

VELOCITY 

QRAt  IN 

RAX 

ALEHA 

0. 

0. 

1 . 

20. 

U5. 

92345.440 

1 2000. 

20. 

•3YSIN  • 

2 

TIRE  TABLE  FOB  V2  TRAJ 

a QE  = 

20 

hre 

ALT 

VELOCITY 

QR AD  IN 

RAX 

ALEHA 

0. 

0. 

1. 

20. 

30. 

61563.6 

12C00. 

20. 

•SYSTN' 

2 

TIRE  TABL 

E FOR  V3  TRAJ 

a QE  = 

20 

tire 

ALT 

VELOCITY 

QRAC  IN 

RAX 

ALPHA 

0. 

0. 

1. 

20. 

20. 

4 1 0 4 2 . V 

12000. 

20. 

•SYSIN  • 

2 

TIRE  TAELF  FOP  V4  TRAJ 

a QE  = 

20 

TIRE 

ALT 

VELOCITY 

QR AD  IN 

RAX 

ALEHA 

0. 

0. 

1. 

20. 

10. 

20521.2 

12C00. 

20. 

•SYSTN* 

2 

TIME  TABLE  FOR  V5  TRAJ 

a QE  = 

20 

TIRE 

AIT 

VELOCITY 

QRAC  IN 

RAX 

ALFHA 

0. 

0. 

1. 

20. 

5. 

10260.0 

12000. 

20. 

•SYSTN  • 

3 

TIRE  TABLE  FOR  VI  TRAJ 

a QE  = 

80 

TIBB 

ALT 

VELOCITY 

QBAD  IN 

RaX 

ALEHA 

0. 

0. 

1. 

20. 

21.  89 

1OCOC0. 

7358.5 

20. 

*s. 

10C000. 

12000. 

20. 

3 

TIRE  TABLE  FOR  V2  TRAJ 

a QE  = 

80 

* » • t 

ALT 

VELOCITY 

QRAD  IN 

RAX 

ALEHA 

9. 

0. 

1. 

20. 

1OCOC0. 

9012. 

20. 

100000. 

12000. 

20. 

-?•!  TABLF  FOR  V3  TRAJ 

« QE  * 

80 

Alt 

VELOCITY 

QRAD  IN 

RAX 

ALEHA 

1.  20. 

1102».5  20. 

* 2000.  20. 
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I 


•SYSIN ' 

2 

TIHE  TABLE  FOR  V4 

TP  AJ 

3 

OF 

= 80 

TIHE 

ALT  VELOC 

IT  Y 

QR  AD 

IN  1 A X 

ALPHA 

0. 

0.  1. 

20. 

10. 

59088.  12000. 

20. 

•S7STN  • 

2 

TIHE  TABLE  FOR  V5 

TP  AO 

3 

OF 

= 80 

TIKE 

ALT  VELOCITY 

QR  AD 

IN  HA  X 

ALEHA 

0. 

0.  1. 

20. 

5. 

29544. C 12000. 

20. 

SYS  I N' 

3 2 0 0 

54. 

3 

65  0 

0 

THERHAL  STRESS  GECHETRY  DATA 

FOR 

LIHIT 

STUDY 

RADIUS 

NODE  HAT 

CODE 

INTERFACE  RADIOS 

.79624 

.79624 

3.  1 

• 

.82904 

.84544 

2.  1 

. 

.86184 

1.  1 

• 

.86184 

SYSIN' 

3 2 0 0 

59. 

3 

65  0 

0 

.79624 

.79624 

6. 

1. 

.82904 

. 84  54  4 

5. 

1. 

.86184 

4. 

1 . 

.86184 

SYSIN'  3 

2 0 

0 59.3 

65  0 

0 

.79624 

.79624 

0, 

1. 

.82904 

.84  54  4 

8. 

1. 

.86184 

7. 

1. 

. 86184 

V 

SYSIN'  .1 

2 

0 

0 

59.3 

65  0 

0 

.79624 

.79624 

12. 

1. 

.82904 

.84544 

11. 

1. 

. 86184 

1C. 

1. 

.86184 

SYSIN'  3 

2 

0 

0 

59.3 

65  0 

0 

.79624 

.79624 

15. 

1. 

.82904 

.84  54  4 

14. 

1. 

.86194 

13. 

1. 

.86184 

•SYSIN' 

7 4 

0 

0 

59.3  65 

0 

0 

.75124 

.75124 

27. 

2. 
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77887 

26. 

2. 

.7881 

79723 

25. 

2. 

81575 

24. 

2. 

.82496 

83418 

23. 

2. 

.84339 

85261 

22. 

2. 

86184 

21. 

2. 

.86184 

•SYS  IN' 

7 4 

0 0 

59.3  65 

75124 

34. 

2. 

.75124 

77889 

33. 

2. 

79723 

22. 

2. 

.7881 

81575 

31. 

2. 

33418 

30. 

2. 

.82496 

85261 

29. 

2. 

. 84339 

86184 

28. 

2. 

•SYSIN' 

7 4 

C 0 

.86184 
59.3  65 

75124 

41. 

2. 

.75124 

77887 

40. 

2. 

79723 

39. 

2. 

.7881 

91575 

33. 

2. 

8 1418 

37. 

2. 

.82496 

95261 

36. 

2. 

. B4339 

96184 

35. 

2. 

•SYSTN' 

7 4 

0 0 

.86194 
59.3  65 

75124 

48. 

2. 

.75124 

77889 

47. 

2. 

79723 

46. 

2. 

.7881 

81575 

45. 

2. 

93418 

44. 

2. 

.82496 

35261 

4 3. 

2. 

.84339 

36104 

42. 

2. 

• SYSTN  • 

7 4 

0 0 

.86194 
59.3  65 

.75124 

.75124  55.  2. 
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,77889 

54. 

2. 

. 7881 

,79723 

53. 

2. 

,91575 

52. 

2. 

.82496 

.83419 

51. 

2. 

.84339 

,95261 

50. 

2. 

,86184 

49. 

2. 

. 86184 

•SYSIN* 

c 

4 

0 

0 

59.3  6 

.77784 

,77784 

65. 

3. 

.79884 

80  0 34 

64. 

3. 

.81984 

,83034 

63. 

3. 

. 8 4 C 84 

,95134 

62. 

3. 

,86184 

61. 

3. 

.86184 

• STSIN • 

5 

4 

0 

0 

59.3  6 

.77784 

77784 

o 

• 

3. 

. 79894 

80934 

69. 

3. 

.81994 

93034 

68. 

3. 

.84084 

85134 

67. 

3. 

86184 

6 £ . 

3. 

.86194 

• S YS  T N * 

5 

4 

0 

0 

59.3  6 1 

.77794 

77784 

75. 

3. 

. 79P84 

B0°34 

74. 

3. 

.81984 

83034 

73. 

3. 

.84084 

35134 

72. 

3. 

86184 

71. 

3. 

.86194 

• SYS  IN • 

5 

4 

0 

0 

59.3  6' 

.77784 

77784 

80. 

1. 

.79884 

90934 

79. 

3. 

.81994 

33034 

78. 

3. 

.84084 

851  34 

77. 

3. 

89< 
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3. 
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0 
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6 
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0 

7 

3E1  0 
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8 

3E10 
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0 

9 

3S  10 

3.89  0 

10 

3E10 

3.89 

0 

11 

3F10 

3.89  0 

12 

3 E 10 

3.  89 

0 

13 

3E 10 

3.89  0 

14 

3E10 

3.89 

0 

15 

3510 

3.89  0 

16 

3E10 

3.89 

0 

17 

3E  10 

3.89  0 

18 

3510 

3.89 

0 

19 

3 El  0 

3.  89  0 

20 

3E10 

3.89 

0 

21 

3E10 

3.89  0 

22 

3510 

3.89 

0 

23 

3E1  0 

3.89  0 

24 

3E1  0 

3.89 

0 

25 

3E10 

3.89  0 

26 

3510 

3.89 

0 

27 

3E10 

3.89  0 

28 

3E10 

3.89 

0 

29 

3E10 

3.89  0 

30 

3E10 

3.89 

0 

5 

3 

ALUMINA  AT  20 

DEGREE  LAUNCH  ANGLE 

1 1 1 

2 3 

1 1 

1 

4 

5 6 

2 2 

2 

7 

e s 

3 3 

3 

10 

1 1 

12 

4 4 

4 

13 

14 

1 5 

5 5 

5 

5 

7 

FUSED  SILICA  AT 

20  DEGREE 

LAUNCH  ANGLE 

2 

2 21 

22 

23 

24 

25 

26 

27 

6 

6 

1 

28 

29 

30 

31 

32 

33 

34 

7 7 

2 

35 

36 

37 

38 

39 

40 

41 

8 

8 

3 

42 

43 

44 

45 

46 

47 

4 8 

9 9 

4 

49 

50 

51 

52 

53 

54 

55 

10  10 

5 

5 

5 

PYRCCERAH  AT 

20  DEGREE 

LAUNCH  ANGLF 

3 

3 

61 

62 

63 

64 

65 

1 1 

11 

1 

66 

67 

68 

6 9 70 

12 

12 

2 

71 

72 

73 

74 

75 

13 

13 

76 

77 

78 

79 

80 

14 

14 

4 

81 

82 

83 

84  85 

15 

15 

5 

ALUMINA  A 1 80  DEGRFE  LAUNCH  ANGLE 

1 1 101  102  103  16  16  6 

104  105  106  17  17  7 

107  108  109  18  18  8 

110  111  112  IS  19  9 

113  114  115  20  20  10 

5 7 

FUSED  SILICA  AT  80  DEGREF  LAONCH  ANGLE 
2 2 121  122  123  124  125  126  127  21  21  6 

128  129  130  131  132  133  134  22  22  7 
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135  136  137  138  139  140  141  23  23  8 

142  193  144  145  146  147  148  24  24  9 

150  151  152  153  154  155  25  25  10 

PYSCCE  8 AH  AT  80  DEGREE  LAUNCH  ANGIE 

3 3 161  162  163  164  165  26  26  6 

166  167  168  169  170  27  27  7 

171  172  173  174  175  28  28  8 

176  177  178  179  180  29  29  9 

181  182  183  184  185  30  30  10 

1 9 

4 5 6 8 10  12  15  18  20  22  25  27  30  32  35  37  40  42  45 
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Appendix  C 


SAMPLE  OUTPUT 


The  listing  included  herein  was  produced  during  the 
sample  problem  execution  and  was  written  on  file  SYS- 
PRINT.  It  contains,  in  addition  to  the  system  generated 
JCL  and  disposition  lines  the  partial  echo  of  the  input 
data  and  a sample  of  the  minimum  print  output. 
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Appendix  £ 

ERROR  AND  WARNING  MESSAGES 

In  addition  to  the  PL/I  error  messages  that  can 
occur  due  to  faulty  execution,  several  messages  have 
been  incorporated  into  the  URLIM  code  that  are  generic 
to  the  radome  limitations  analyses.  The  PL/I  messages 
that  may  occur  are  described  in  Ref.  2.  The  OKLIM  com- 
putational messages  are  listed  below. 

1.  1 SUBROUTINE  CYL2D:  CALCULATED  NODE  I LARGER 

THAN  ALLOWED  (LASCAP),  EXECUTION  HAY  BE 
FAULTY.  CONTROL  RETURNED  TO  MAIN  PROGRAM' 

Cause 

A control  value  has  been  given  to  the  two 
dimensional  cylindrical  network  subroutine 
(CYL2D)  for  specifying  the  number  of  nodes 
within  the  section  and  a resultant  node  number 
has  been  computed  that  is  larger  than  the 
control  variable  INODES. 

Effect 

The  subroutine  is  halted  from  further  com- 
putation and  the  URLIH  program  continues.  The 
thermal  network  that  results  may  not  be  the 
one  intended. 

Solution 

Either  make  INODES  larger  so  as  to  accommodate 
the  number  of  nodes  actually  needed  or  reduce 
the  control  parameters  to  the  two  dimensional 
cylindrical  subroutine. 

2.  'SUBROUTINE  TWOD:  CALCULATED  NODE  # LARGER  THAN 

ALLOWED  (LASCAP),  EXECUTION  MAY  BE  FAULTY. 
CONTROL  RETURNED  TO  MAIN  PROGRAM' 

Cause 

Same  as  for  message  1 above 
Effect 

Same  as  for  message  1 above 
Solution 

Same  as  for  message  1 above 
-141- 
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3.  'THE  FOLLOWING  TEMPERATURE  ENTRY  WAS  INCONSIS- 

TENT WITH  THE  FIRST  SET  OF  PRESSURE  ENTRIES. 
CHECK  THE  TABULAR  INPUT'  --  from  subroutine 
READAIR. 

Cause 

In  subroutine  READAIR  one  set  of  temperature 
values  input  does  not  match  the  first  set  of 
temperatures  (it  is  required  that  each  set  of 
temperatures  input  at  each  value  of  pressure 
be  identical). 

Solution 

Examine  the  tabular  input  data  and  be  sure  it 
is  consistent  with  the  requirements. 

4.  'THE  FOLLOWING  MACH  # ENTRY  WAS  INCONSISTENT 

WITH  THE  FIRST  SET  OF  ANGLE-OF- ATTACK  ENTRIES. 
CHECK  THE  TABULAR  INPUT'  --  from  subroutine 
READAFC. 

Cause 

In  subroutine  READAFC  one  of  the  Mach  number 
entries  did  not  match  the  first  set  of  Bach 
number  entries. 

Solution 

Examine  the  tabular  input  data  and  be  sure  it 
is  consistent  with  the  requirement.;. 

5.  'ITERATION  FOR  T LOCAL  IN  SAERO  BLOCK  XXX  HAS 

EXCEEDED  xix.  T LOCAL  = XXX,  DELTA  TL  = xxx, 
AND  ENERGY  IMBALANCE  = xxx*  — from  subroutine 
SAERO. 

Cause 

The  solution  technique  for  calculating  the 
local  temperature  behind  the  shock  wave  did 
not  converge  to  the  required  limits.  The 
value  of  TL  that  is  printed  out  is  used  and 
the  program  continues. 

Solution 

Increase  the  number  of  allowed  iterations 
and/or  the  tolerance  values  for  the  conver- 
gence. Check  the  velocity  and  altitude  tables 
for  anomolies  and  compare  with  expected 
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values.  Check  the  Mach  ratio  and  pressure 
ratio  input  tables. 

6.  • ITERATION  ON  SUBFACE  NODE  # XXX  HAS  EXCEEDED 

XXX,  TSURF  = XXX  AND  DELTA  TSURF  = xxx.  HEAT 
IMBALANCE  = xxx ' --  from  subroutine  SAEBO. 

Cause 

The  surface  energy  balance  routine  has  failed 
to  find  a solution  within  the  thermal  and  heat 
flux  tolerance  levels  input.  The  printed 
value  of  TSURF  will  be  used. 

Solution 

Raise  either  the  number  of  iterations  or  the 
convergence  criteria,  check  input  and  cal- 
culated data  for  erroneous  or  excessively  high 
heating  rates,  and  check  for  inconsistent  air, 
altitude  and/or  material  property  data.  Check 
also  that  the  temperatures  are  not  beyond  the 
extent  of  the  air  property  tables. 

7.  'AN  ERROR  IN  THE  COMCCN  DATA  HAS  BEEN  ENCOUN- 

TERED: CAPACITOR  XXX  OR  CAPACITOR  XXX  HAS  BEEN 
CONNECTED  TO  ITSELF  THE  COMCON  DATA  IS: 

I I A (I)  IB  ( I)  CT1  (I)  CT  2 (I) 

XX  XX  XX  XX  XX 

THE  SEARCH  FOR  A TEMFER ATURE  FOE  RADIUS  = XXX 
HAS  CONTINUED*  — from  subroutine  SIGMA 

Cause 

The  thermal  network  data  contains  an  inconsis- 
tent nodal  interconnection,  namely  a node  has 
been  specified  as  connected  to  itself.  The 
error  was  detected  while  determining  the 
contact  temperature  between  nodes  of  different 
materials  by  the  thermal  stress  routine. 

Solution 

Inspect  the  diagnostic  information  printed  and 
correct  the  thermal  model  information.  IA  and 
IB  are  lists  of  the  capacitor  numbers  that  are 
connected  through  a contact  resistance  and  CT1 
and  CT2  are  the  respective  temperatures  at 
either  side  of  the  contact  surface. 

8.  'THE  COMCON  DATA  HAS  BEEN  SEARCHED  & A TEM- 
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PER ATURE  VALUE  FOR  RADIUS  = XXX  HAS  NOT  BEEN 
FOUND.  CONTROL  HAS  BEEN  RETURNED  TO  THE  SHTP 
FROM  SUBROUTINE  SIGMA' 

Cause 

In  the  nodal  network  data  a contact  resistance 
connection  has  been  incorrectly  specified  or 
the  input  data  (node  number  versus  radius)  is 
in  error.  Inspect  these  input  data  for  proper 
assignments  of  interface  locations. 

Effect 

The  thermal  stress  calculation  for  this  case 
is  terminated  and  the  main  program  continues. 

9.  'ALL  MODULUS  COMPUTATIONS  WERE  ZERO,  SIGMA 

TERMINATED  ABNORMALLY'  — from  subroutine 
SIGMA. 

Cause 

The  elastic  modulus  for  the  entire  wail  thick- 
ness of  a thermal  stress  case  was  computed  to 
be  zero.  Either  the  temperatures  are  hrgh 
enough  to  produce  melting  or  the  Modulus 
versus  Temperature  tables  are  in  error. 

Solution 

In  the  cases  where  the  temperatures  are  all 
past  melting,  this  is  actually  a normal  return 
from  the  thermal  stress  routine.  Otherwise, 
inspect  the  Modulus  versus  Temperature  input 
data  for  errors. 

Effect 

Control  is  given  back  to  the  main  program  and 
execution  continues. 

10.  ****THE  ABOVE  GRADIENT  IS  NOT  MONOTONIC.  IT 
IS  POSSIBLE  THAT  THE  INTERPOLATION  ROUTINE  HAY 
HAVE  TRUNCATED  A PEAK***'  --  from  subroutine 
SIGMA. 

Cause: 

The  distribution  of  temperature  versus  radius 
through  a wall  section  is  not  monotonic;  the 
interpolation  routine  (AIT)  may  have  truncated 
the  distribution. 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

Laurel  Maryland 


Solution: 

Inspect  the  tabulated  temperature  versus 
radius  data  for  such  a peak. 

Effect : 

Execution  continues,  the  possibly  truncated 
function  is  used. 

11.  'INTEGRATION  ROUTINE  HAS  EXCEEDED  20  ITERA- 
TIONS. THE  RATIO  = xxx ' --  from  subroutine 
SIGHA. 

Cause 

In  calculating  either  elastic  modulus  or 
thermal  strain  values  in  the  thermal  stress 
routine,  successive  approximations  to  the 
value  of  an  integral  have  not  converged  to 
within  0.5*  in  subroutine  TRAP. 

Solution: 

The  ratio  of  successive  approximations  is 
printed  and  can  be  examined  for  accuracy.  If 
the  ratio  is  unsatisfactory,  the  free  thermal 
expansion  and/or  elastic  modulus  tables  can  be 
expanded  or  made  smoother  by  adding  entries. 

Ef  f ect : 

The  program  continues  and  uses  the  inaccurate 
value  for  the  integral. 

12.  'TIME  STEP  LESS  THAN  MINIMUM  ALLOWED  CRITICAL 
INDEX  NO.  = xx  DELTIM  = xxx'  — from  subrou- 
tine STEP. 

Cause: 

The  stability  time  step  was  calculated  to  be 
less  than  the  minimum  allowed  by  control 
variable  HINSTEP. 

Solution: 

The  thermal  model  must  be  altered  so  that  the 
ratio  of  total  capacitance  to  total  conduc- 
tance for  the  critical  node  is  larger  than  the 
minimum  allowable.  Alternatively,  a smaller 
allowable  time  step  can  be  specified. 

Effect : 

The  URLIM  execution  is  halted. 
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13.  * STEADY  STATE  SOLUTION  HAS  EXCEEDED  xxx 

ITERATIONS. 

THE  NOT  YET  STEADY  TEMPERATURES  ARE: 

[ A list  of  nodes,  temperatures  and  heat 
rates  ] 

— from  subroutine  STEADY. 


Cause: 

In  the  routine  STEADY,  the  convergence 
criteria  have  not  been  met.  In  the  message, 
the  node  number  and  current  temperature  are 
listed  and  the  last  correction  to  the  tem- 
perature and  the  net  heat  flow  for  the  node 
are  also  shown,  only  the  not  yet  converged 
nodes  are  printed  (i.e.,  those  nodes  satis- 
fying the  convergence  criteria  are  not  printed 
in  the  above  table) . If  PUNCH  * 1 then  all  of 
the  temperatures  will  be  saved  (as  described 
in  Appendix  N)  in  spite  of  the  incorrect 
solution.  Before  stopping  the  execution  of 
the  subroutine,  whether  convergence  is  met  or 
not,  a call  to  subroutine  WRITE  is  maae  so 
that  a printing  of  all  of  the  temperatures  is 
obtained. 

It  is  possible  that  the  solution  may  be 
intrinsically  unstable  with  the  initial  guess 
as  supplied.  If  this  is  true  a random  oscil- 
lation of  the  temperatures  may  result  (with 
the  iteration  limit  exceeded)  or  the  tempera- 
tures will  increase  (or  decrease)  in  magnitude 
until  one  of  the  temperatures  becomes  larger 
(or  smaller)  than  can  be  represented  by  the 
computer  on  which  the  program  is  running. 

When  this  occurs  an  overflow  or  underflow 
error  message  will  be  printed  and  execution 
will  be  stopped  without  a printout  of  tem- 
perature. These  underflow  or  overflow  mes- 
sages are  not  a part  of  the  URLIM  program. 

Effect : 

The  program  execution  is  halted. 
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Appendix  F 

THE  READ  ROUTINES 


The  READ  routines  are  used  by  the  URLIM  program  to 
read  tabular  data  in  EDIT  format  from  card  files.  Each 
card  file  has  the  general  format  requirement  of  two 
comment  cards  followed  by  the  tabular  data.  In  many  of 
the  routines  the  identification  number,  the  sequence  of 
the  arguments  and  the  order  of  the  data  on  the  cards  are 
important  in  specifying  required  data  to  other  sections 
of  the  program.  The  data  are  required  to  be  in  specific 
column  fields  and  are  either  in  F-  or  E-format.  The 
specific  formats  will  be  listed  with  each  subroutine 
description.  Each  READ  routine  prints  the  data  it  reads 
on  the  output  file  SYSPRINT  in  a slightly  reformatted 
manner.  The  first  line  of  output  serves  to  identify  the 
specific  subroutine  and  the  identification  number  as- 
sociated with  the  particular  group  of  data,  when  appro- 
priate. The  next  two  lines  simply  echo  the  two  comment 
cards  that  are  input  at  the  beginning  of  each  table. 

The  data  table  then  fellows. 

READRK 

Subroutine  READRK  is  used  in  the  main  program  of 
URLIM  to  read  thermal  transport  and  dielectric  constant 
material  properties.  The  identification  number  will 
associate  the  properties  read  in  by  READRK  with  the 
property  code  numbers  used  elsewhere  in  the  program. 

The  list  of  arguments  for  READRK  must  be  in  a specific 
order;  i.e.,  READRK  must  be  called  as  follows: 

CALL  READRK  (ID#,  RKTEMP,  EHOCP,  K,  EMISS,  DIEL, 

NPTS , FILENAME)  ; 


is  the  material  code  for  the  properties 
which  follow  in  the  calling  sequence 
(FIXED  BIN)  ; 

is  the  vector  of  Rankine  temperatures  and 
must  be  the  second  argument  ((*)  FLOAT 


where: 

ID# 

RKTEMP 
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BIN),  of  dimension  > NPTS  ♦ 1; 

RHOCP  is  the  vector  of  density  times  specific 

heat  values  and  must  be  the  third  argument 
((*)  PLOAT  BIN)  , of  dimension  > NPTS; 

K is  the  vector  of  thermal  conductivity  data 

and  must  be  the  fourth  argument  ((*)  FLOAT 
BIN),  of  dimension  > NPTS; 

EM IS S is  the  vector  of  IR  emissivity  factor 

values  and  must  be  the  fifth  argument  ( (*) 
FLOAT  BIN),  of  dimension  > NPTS; 

DIEL  is  a vector  of  dielectric  constant  values 

((*)  FLOAT  BIN),  of  dimension  > NPTS; 

NPTS  is  the  number  of  tabular  data  caras  read 

in  by  RADRK  (FIXED  BIN)  (all  the  READ 
routines  must  have  2 comment  cards  ahead 
of  the  tabular  cards);  and 

FILENAME  is  the  name  of  the  file  from  which  the 
card  images  will  be  read  (FILE) . 

The  data  on  the  card  file  named  by  FILENAME  must  be  in 
the  following  format: 

2 comment  cards  (160  characters)  followed  by; 


Columns 


Value 


1-14 

15-29 

30-44 

45-69 

60-69 

70-80 


Values  of  temperature  (°R)  , F-format; 

Values  of  density  times  specific  heat 
(Btu/ft 3-°R)  , F-format; 

Values  of  thermal  conductivity  (Btu/ft-h- 
°R) , F-format; 

Values  of  IR  emissivity,  F-format; 

Values  of  RF  dielectric  constant,  F-for- 
mat ; and 

Values  of  RF  loss  tangent,  F-format. 
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If  it  is  desirable  to  have  the  values  of  the  ther- 
mal property  variables  set  by  some  other  method  than 
reading  the  above  described  table,  a call  to  subroutine 
READRK  must  still  be  made.  In  such  a situation,  the 
number  of  data  points,  NPTS,  is  set  to  zero  and  the  call 
is  made;  this  mill  serve  to  associate  the  variables  in 
the  calling  sequence  with  the  material  identification 
number  but  no  cards  will  be  read  from  the  file  FILENAME. 

It  should  be  noted  that  READRK  vill  store  the  value 
NPTS  as  the  first  entry  to  the  independent  variable, 
RKTEMP  (unless  NPTS  = 0)  ; this  is  consistent  with  the 
data  management  requirements  of  the  interpolating 
routines  used  to  extract  values  from  the  tables,  (c.f. 
"Interpolation  Techniques,"  Section  2 of  Volume  1). 


READFH 


Subroutine  READFH  reads  mechanical  property  data 
that  is  dependent  on  the  Farenheit  scale.  The  calling 
sequence  and  data  table  format  are  as  follows: 

CALL  READFH  (ID,  FHTEHP,  EMOD,  NU,  EXPANS,  ENI.1, 
NPTS,  FILENAME)  ; 


where  the  variable  names  are  associated  with  the  data 
table  entries  on  the  file  FILENAME  with  the  following 
format: 

Item  Columns  Parameter  and  Format 


ID 

FHTEMP 


EMOD 


NU 


* Table  identifying  number; 

1-14  Independent  variable  temperature 

(®F) , (*)  FLOAT  BIN,  dimensioned  > 
NPTS  ♦ 1,  F-format; 

15-29  Values  of  elastic  modulus  (psi) , 

(*)  FLOAT  BIN,  dimensioned  £ NPTS, 
E-format ; 

30-44  Values  of  Poisson's  ratio,  (*) 
FLOAT  BIN,  dimensioned  > NPTS, 
F-format ; 
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Item 


Columns  Parameter  and  Format 


EXPANS  45-59  Values  of  free  thermal  expansion 

(in. /in.)  FLOAT  BIN,  dimensioned  > 
NPTS , E-fcrmat; 

ENU  - Values  of  EMOD/ (1-NU)  calculated  by 

the  subroutine  (i.e.,  not  read)  and 
printed  on  the  output,  FLOAT  BIN, 
dimensioned  > NPTS; 


NPTS  - The  number  of  tabular  data  cards 

read  (recall  that  all  the  READ 
routines  require  two  additional 
comment  cards  ahead  of  the  data 
cards  that  are  not  counted  as  data 
cards)  ; and 

FILENAME  - The  name  of  a file  that  contains 

the  NPTS  ♦ 2 card  images  read  by 
FEADFH. 


READFH  stores  NPTS  as  the  first  value  of  FHTEMP, 
and  then  stores  the  tabular  values  of  temperature; 
therefore  FHTEMP  must  be  dimensioned  2:  NPTS  ♦ 1. 


READCP 


Subroutine  READCP  reads  data  from  a file  that 
describes  part  or  all  of  the  thermal  interconnections 
between  nodes.  These  interconnections  may  be  via 
Fourier  conduction  within  the  same  material,  Fourier 
conduction  between  nodes  of  different  materials,  or 
radiation  heat  transfer.  Any  number  of  interconnections 
can  be  specified  for  any  one  node.  The  calling  sequence 
for  READCP  is: 

CALL  READCP  (FILENAME)  ; 

where  FILENAME  is  the  name  of  the  file  where  the  network 
data  are  written  and  has  the  FILE  attribute.  The  data  on 
the  file  FILENAME  must  be  in  EDIT  format  as  follows: 
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Item 

Columns 

; Value 

1 

1-4 

Thermal  node  number,  F-format; 

2 

10-19 

Thermal  node  volume  (ft3) , E-for- 
mat; 

3 

20-29 

Thermal  material  identification 
number,  F-format; 

4 

30-39 

The  number  of  thermal  interconnec- 
tions for  this  node,  F-format; 

5 

40-49 

The  node  number  of  the  first  inter- 
connecting node,  F-format; 

6 

50-59 

The  area-to-length  ratio  of  the 
thermal  path  between  the  centers  of 
the  two  nodes  (items  1 and  5 above) 
(ft)  , E-fcrmat; 

7 

60-69 

The  node  number  of  the  second 
interconnecting  node,  F-foLmat; 

8 

70-79 

The  area-to-length  ratio  of  the 
thermal  path  between  the  centers  of 
the  two  nodes  (items  1 and  7 above) 
(ft)  , E-fcrmat. 

The  number  of  cards  read  by  FEADCP  is  not  restric- 
ted but  the  end  of  the  input  is  signalled  by  including  a 
blank  card  at  the  end  of  the  data. 

The  following  example  illustrates  the  use  of  the 
eight  item  data  field  in  defining  a node  (No.  23)  con- 
nected to  two  other  nodes  (Nos.  24  and  206) : 

23  4.8E-2  5.  2.  24.  2.301E0  206.  1.8E-1 

This  example  entry  shows  node  23  as  having  a volume 
equal  to  0.048  ft3,  being  composed  of  material  number  5, 
and  having  two  interconnections:  one  with  an  A/L  ratio 
to  node  24  of  2.301  ft  and  an  A/L  ratio  to  node  206  of 
0.18  ft. 

For  radiation  heat  transfer  the  interconnecting 
node  numbers  (item  5 cr  7)  are  input  with  minus  signs 
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and  the  A/L  ratios  (either  item  6 or  8)  are  interpreted 
as  being  the  area  tiles  view  factor  between  the  nodes. 

In  order  to  describe  a coaposite  thermal  intercon- 
nection (i.e.,  nodes  of  differing  material  identifica- 
tion numbers  with  a contact  resistance  at  the  interface) 
an  alteration  to  the  input  conventions  is  required.  In 
this  case  the  iteas  that  are  different  than  before  are 
as  follows: 

Itea  Columns  Value 


3 20-29  A composite  material  code  tnat  is  a 

three  digit  number  identifying  the 
material  from  the  node  (item  1)  to 
the  other  node  (item  5)  with  the 
contact  resistance  values  represen- 
ted as  the  center  digit,  F-fofmat; 

5 40-49  The  node  number  of  the  compositely 

interconnected  node,  F-format; 

6 50-59  The  area -to- length  ratio  from  the 

primary  node  center  to  the  inter- 
face (ft)  , E-tormat; 

7 60-69  The  area  cf  the  interface  (ft2)  , 

F-format;  and 

8 70-79  The  area  to  length  ratio  from  the 

interface  to  the  center  of  the 
connecting  node  (ft)  , E-format. 

Note  that  the  contact  resistance  values  are  assumed  to 
have  a unique  material  property  code  and  must  be  input 
via  a call  to  READRK.  The  values  read  in  should  be 
conductances  (Btu/f t* -h-°R)  so  that  when  multiplied  by 
the  contact  area  (Itea  7 above)  a true  conductance 
(Btu/h-°R)  is  achieved.  (In  some  cases  the  contact 
resistance  is  assuaed  small  and  a very  large  nuaber 
(1.0E9)  can  be  input  in  Item  7 along  with  any  arbitrary 
aaterial  code  for  the  interface  so  that  a large  conduc- 
tance for  the  interface  is  achieved.) 


. 

* 
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The  following  example  will  show  the  form  of  a com- 
posite connection  description: 

23  4.8E-2  452.  1.  24.  2.3E0  20.6  2.3S0 

This  shows  node  23  (the  first  item)  connected  to  node  24 
(the  fifth  item)  through  a composite  or  contact  connec- 
tion (by  virtue  of  the  three  digit  third  item).  The 
first  digit  (4)  of  this  composite  material  code  is  for 
node  23,  the  second  digit  is  for  tne  conductance  of  the 
contact,  and  the  third  digit  is  for  the  compositely 
connected  node,  24.  The  path  length  from  node  23  to  the 
contact  is  item  6 (2.3  ft),  the  contact  area  is  item  7 
(20.6  ft2),  and  the  path  length  to  node  24  from  the 
contact  is  item  8 (2.3  ft). 

In  the  general  case  where  one  node  must  be  connec- 
ted to  more  than  two  ether  nodes,  or  must  be  connected 
to  several  nodes  via  composite  connections,  a one-card 
entry  offers  insufficient  room  for  all  the  data.  In 
such  cases  the  number  of  connections  entry  (Item  4)  is 
set  to  some  number  larger  than  the  total  number  of 
active  thermal  nodes.  The  rest  of  the  entries  are  then 
set  as  necessary  according  to  the  conventions  mentioned 
above.  A second  card  entry  will  then  be  required  whose 
first  entry  (Item  1)  will  be  the  same  value  as  the 
previous  number-of -connections  entry.  This  convention 
is  then  seen  to  be  an  indirect  addressing  scheme  wherein 
the  number-of -connections  entry  refers  to  a card  number 
where  more  interconnection  data  may  be  found.  This 
technique  may  be  used  to  any  number  of  levels  and  is 
obviously  terminated  when  a standard  n umber-of-connec- 
tions  is  encountered  (i.e.,  a value  of  1 or  2)  . It  is 
very  important  that  the  indirectly  addressed  cards  have 
Item  1 entries  larger  than  the  highest  active  capacitor 
number.  To  insure  proper  allocation  of  storage  for  all 
of  the  network  information  the  proper  call  should  be 
made  to  subroutine  STOEE  (c.f.  Appendix  I).  To  illus- 
trate the  use  of  the  indirect  addressing  scheme  consider 
the  following  two  entries: 

23  4.8E-2  452.  101.  24.  2.3E0  20.6  2.3E0 

101  4.  2.  22.  1.6E0  32.  3.6E0 

In  this  example  the  card  for  node  23  has  its  num- 
ber-of -connections  entry  set  to  101;  the  next  card  has 
for  its  "node”  number  the  value  101  and  will  be  unders- 
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tood  to  contain  further  connection  data  for  node  23, 
namely  that  nodes  22  and  32  connect  via  material  4 to 
node  23  with  the  appropriate  conduction  paths. 

Another  situation  that  often  arises  is  when  a 
thermally  anisotropic  material  is  used.  In  such  cases 
the  thermal  conductivity  is  different  in  different 
spatial  directions.  The  URL1M  program  considers  only 
two-dimensional  anistropy  and  therefore  two  sets  of 
thermal  conductivity  data  are  needed.  To  indicate  this 
the  material  property  entry  (Item  3)  will  be  a two-digit 
number,  the  first  of  which  indicates  the  material  code 
for  the  connection  described  by  Items  5 and  6 and  the 
second  digit  will  apply  to  the  connection  described  by 
Items  7 and  8.  In  this  way  two  thermal  conductivities 
can  be  applied  for  the  two  different  thermal  path  direc- 
tions. 

Any  of  the  above  entry  conventions  for  network 
description  can  be  combined  in  any  manner  on  indirect 
address  cards. 


READAIR 


Subroutine  READAIB  must  be  used  to  read  in  the 
temperature-  and  pressure-dependent  properties  of  air 
when  aerodynamic  heating  is  to  be  used.  The  subroutine 
is  invoked  by  the  following  CALL  sequence  to  read  values 
of  pressure,  temperature,  enthalpy,  Prandtl  number, 
specific  heat  ratio,  viscosity,  and  compressibility: 

CALL  READAIR  (#PRES,  #TEH PS,  CPRES,  CTEHP,  CENTH, 
FILENAHE) ; 

where: 

iPRES  is  the  number  of  different  pressure  values 
for  the  bivariant  data  (FIXED  BIN)  , 

9TENPS  is  the  number  of  different  temperature 

values  for  the  bivariant  data  (FIXED  BIN) , 

CPRES  is  a conversion  factor  that  will  be  mul- 

tiplied with  each  pressure  value  read 
(this  number  represents  the  conversion 
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required  to  obtain  lbf/ft2  from  whatever 
units  the  pressure  values  are  supplied  in) 
(FLOAT  BIN)  , 

CTEMP  is  a similar  factor  for  obtaining  °B  from 

the  the  temperature  values  as  supplied 
(FLOAT  BIN) , 

CENTH  is  a similar  conversion  factor  for  obtain- 
ing Btu/lbm-°B  from  the  enthalpy  values  as 
supplied  in  the  data  (FLOAT  BIN)  , 

FILENAME  is  the  variable  naming  the  file  where  the 
data  are  to  be  found  (FILE)  . 

The  cards  read  by  BEADAIE  must  have  the  following  data 
in  the  following  format: 

Two  cards  used  for  comment  followed  by  data  caLds: 


Columns  Variable  and  Format 


I- 10  Pressure,  F-format; 

II- 20  Temperature,  F-fcrmat; 

21-30  Enthalpy,  F-format; 

31-40  Prandtl  number,  F-format; 

41-50  Specific  heat  ratio,  F-format; 

51-60  Viscosity  (lbm/ft-s);  and 

61-70  Compressibility , F-format. 

The  two  independent  variables  must  be  in  monotonically 
increasing  order;  i.e.,  all  values  at  the  lowest  pres- 
sure must  be  first  followed  by  the  next  largest  pres- 
sure, etc. , and  all  cards  for  the  same  pressure  value 
must  have  the  temperatures  in  increasing  order.  The 
routine  will  expect  (#PBES  times  #TEMPS)  plus  two  cards 
to  be  in  the  file  FILENAME. 

The  values  read  are  stored  in  several  external 
variables  whose  management  is  cf  no  direct  concern  to 
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the  user.  The  only  requirement  is  that  if  the  routine 
SAERO  is  used  for  aerodynamic  heating  then  READA1R  must 
be  called  to  read  the  air  properties. 


READAL 


READAL  is  a subroutine  that  reads  from  a file  the 
altitude-dependent  properties  of  the  atmosphere,  namely 
the  ambient  temperature,  pressure,  speed  of  sound,  and 
effective  space  temperature  as  functions  of  altitude 
(expressed  in  feet) . READAL  is  called  as  follows: 

CALL  READAL (#ENTS,  CPPES,  CTEHP,  FILENAME); 
where: 


#ENTS  is  the  number  of  entries  to  he  read  in. 


and  the  other  values  have  the  same  definitions  as  for 
READAIR  but  apply  to  the  data  read  by  READAL.  The  data 
format  for  READAL  is  tabular  Fixed  Point  (F-format) , 15 
column  widths  preceded  by  two  comment  cards: 


Columns 

Value 

1- 

15 

Values 

of 

a ltit 

variable  for  th 

16 

-30 

Values 

of 

a mbie 

31 

-45 

Values 

of 

a mbie 

46 

-60 

Values 

of 

the  s 

61 

-75 

Values 

cf 

the  e 

space  : 

Eor 

r adia 

ude,  the  independent 
is  table  (ft)  ; 

nt  temperature; 

nt  static  pressure; 

peed  of  sound  (ft/s)  ; and 

ffective  temperature  of 
tion  relief  (°R) . 


The  data  read  must  be  in  increasing  order  with 
respect  to  altitude  and  will  be  stored  in  external 
variables  outside  of  the  user's  awareness.  The  later 
use  of  subroutine  SAERO  requires  the  reading  of  the 
altitude  table  by  routine  READAL. 
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READFL 


Subroutine  READFl  has  the  function  of  reading  in 
the  Mach  and  pressure  ratio  data,  as  functions  of  Mach 
number.  These  data  are  required  by  SAERO  to  determine 
the  local  flow  conditions  at  the  station  of  interest  on 
the  flight  body  so  that  heat  transfer  can  be  calculated. 
The  READFL  subroutine  is  invoked  by  the  following 
statement: 


CALL  READFL  (ID#,  MACH#,  MACHPAT,  PRAT,  VGRAD,  NPTS, 
FILENAME) ; 


where: 

ID# 

MACH# 


identifies  the  table  of  data  that  are  read 
(FIXED  EIN)  ; 

is  the  list  of  Mach  numbers  (the  indepen- 
dent variable)  ( (*)  FLOAT  BIN)  ; 


MACHPAT  is  the  list  of  Mach  ratios  that  exist  from 
the  free  stream  to  the  local  condition  on 
the  body  ((*)  FLOAT  BIN); 


PRAT  is  the  list  of  pressure  ratios  that  exist 

from  the  free  stream  to  the  local  con- 
dition on  the  body  ((*)  FLOAT  BIN); 


VGRAD  is  the  list  of  velocity  gradient  terms 
that  exist  at  the  stagnation  point  of  a 
body  in  flight  (in  cases  where  stagnation 
point  heating  is  not  required,  these 
values,  though  read  in,  will  not  be  used 
((*)  FLOAT  BIN)  ; and 

NPTS  is  the  number  of  data  points  in  the  table 

to  be  read  (FIXED  BIN). 

The  data  on  the  file  FILENAME  should  have  two  com- 
ment cards  followed  by  NPTS  tabular  data  cards: 


-157- 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREk  MAHVLANO 


Column 

Variable  and  Format 

1-14 

MACH#,  F-format; 

15-29 

HACHRAT,  F-format; 

30-44 

PRAT,  F-format;  and 

45-59 

VGRAD,  F-format. 

If  any  VGRAD  data  entry  is  left  blank,  a value  of  1 
will  be  assumed  for  that  entry. 


READAFC 


The  READAFC  routine  reads  the  bivariant  aerodynamic 
force  coefficient  data  when  aerodynamic  loads  are  to  be 
calculated  by  subroutine  AERLOAD.  The  data  reaa  are  the 
Mach  number  and  angle  of  attack  dependent  values  of 
missile  normal  force  coefficient,  radome  axial  and 
normal  force  coefficient,  and  axial  location  of  center 
of  pressure.  The  routine  is  invoked  as  follows: 

CALL  READAFC  (#ALPHA,  iMACH,  FILENAME); 


w here : 

tALPHA  is  the  number  of  angle  of  attack  entries 
in  the  tabular  data  (FIXED  BIN)  ; 

#MACH  is  the  number  of  Mach  number  entries  in 
the  tabular  data  (FIXED  BIN) ; and 

FILENAME  is  the  name  of  the  file  where  the  data 
will  be  found. 

The  data  on  file  FILENAME  should  include  two  comment 
cards  followed  by  the  data  cards  in  monotonically 
increasing  order  of  angle-of -attack  and  Mach  number  as 
follows: 

Columns  Variable  and  Format 
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Columns  Variable  and  Format 

1-13  Angle  of  attack  (deg),  F-format; 

14-26  Mach  number,  F-format; 

27-39  Normal  force  coefficient  of  the  missile, 

F-format ; 

40-52  Normal  force  coefficient  of  the  radooe 

alone,  F-format; 


53-65  Axial  force  coefficient  of  the  raaome 

alone  (i.e.,  without  base  drag  effects), 
F-format; 


66-78  Axial  location  of  the  center  of  pressure, 

normalized  to  the  total  length  of  the 
radome,  F-format. 


The  data  read  are  stored  in  external  variables  whose 
management  is  of  no  direct  concern  to  the  user. 


READTH 


The  READTH  subroutine  can  be  used  to  read  in  tables 
of  time-dependent  data.  The  routine,  like  many  of  the 
other  READ  routines,  has  general  utility  for  reading 
arbitrary  data  but  must  be  used  in  a specific  form  if 
the  aerodynamic  heating  routine  SAERO  is  used.  For  such 
cases  the  routine  REAETH  is  called  as  shown: 

CALL  READTH  (ID,  TIHE , ALT,  VEL,  QRAD,  DUM,  DUH, 

NPTS  , FILENAHE)  ; 

where: 


ID 

is  the  trajectory  number  which  will 
used  later  by  SAERO  (FIXED  BIN)  ; 

be 

TIHE 

the  list  of  times  during  the  flight 
((*)  FLOAT  BIN)  ; 

(s) 
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ALT  is  the  list  of  altitude  values  (ft)  ((*) 

FLOAT  BIN)  ; 

VEL  is  the  list  of  freestream  velocities 

(ft/s)  ( (*)  FLOAT  BIN)  ; 

QR AD  is  a list  of  independently  supplied  heat 

flux  values  that  Bay  be  applied  to  the 
SAEKO  surface  (Btu/ft*)  ((*)  FLOAT  EIN) ; 

DUN  is  a suitably  dimensioned  dummy  variable 

or  a valid  variatle  for  use  by  another 
subroutine  called  in  the  main  proyrain;  and 

NPTS , are  as  before  except  that  they  apply  to 

FILENAME  the  data  read  by  READTM. 

The  data  on  file  FILENAME  should  be  of  the  following 
form : 

T vo  comment  cards  followed  by  NPTS  tabular  data  cards 
with  the  following  format: 


Columns  Variable  and  Format 


1-9  TIME,  F-format; 

10-19  ALT,  F-format; 

20-29  VEL,  F-format; 

40-49  URAB,  F-format;  and 

50-59  DUN  or  any  other  variable,  F-format. 

The  first  variable,  TIME,  must  be  a list  of  at 
least  NPTS  ♦ 1 length  because,  as  with  other  READ 
routines,  the  READTM  routine  places  the  value  NPTS  as 
the  first  entry  to  TIME  and  then  the  actual  time  values 
follow. 

If  it  is  desired  not  to  read  the  trajectory  data 
for  a subsequent  SAERO  call  via  READTM  then  a call  to 
READTM  must  still  be  made  but  NPTS  is  set  to  zero.  This 
null  CALL  statement  will  cause  the  identification  number 
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to  be  correctly  associated  with  the  variable  names  in 
the  argument  list  but  no  access  will  be  attempted  to  the 
file  FILENAME. 


READSG 


The  READSG  routine  reads  data  that  describes  the 
cylindrical  geometry  and  temperature  information  re- 
quired by  the  cylindrical  thermal  stress  routines,  SIGMA 
and  SIGMET.  The  information  is  in  the  form  of  tables  of 
thermal  node  number  versus  radius  plus  an  indication  of 
the  finite  difference  mesh  size  for  the  thermal  stress 
solution.  The  subroutine  is  called  as  follows: 


CALL  READSG  (ID#,  RADII,  XINTF,  IND,  MAT,  NREGS, 

NHTP  , FILENAME)  ; 

where  (c.f.  Fig.  11  fcr  nomenclature  conventions): 

ID#  identifies  the  table  of  values  (FIXED 

BIN)  ; 

RADII  is  a list  of  radii  (in.)  at  which  tem- 

peratures will  also  be  given  ( (*)  FLOAT 
BIN)  ; 

X I NT  R is  a list  of  radii  (in.)  that  define  the 
boundaries  of  the  finite  difference  mesh 
through  the  cylinder's  thickness  (there 
should  be  NREGS  ♦ 1 such  radii)  ((*)  FLOAT 
BIN)  ; 

IND  is  a list  of  thermal  node  numbers  that 

correspond  one-to-one  with  the  entries  in 
RADII  above  ((*)  FIXED  BIN); 

MAT  is  a list  of  mechanical  property  code 

numbers  that  are  used  by  subroutine  SIGMA 
and  correspond  one-to-one  with  the  radius 
entries  of  RADII  above  (the  numbers  serve 
to  indicate  which  material  properties 
exist  through  the  thickness  of  the  cylin- 
der) ((*)  FIXED  BIN); 
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NREGS  is  the  number  of  subdivisions  that  are  to 
be  made  through  the  thickness  of  the 
cylinder  (FIXED  El  N)  ; 

NHTP  is  a flag  for  indicating  whether  or  not 

the  data  read  are  to  be  used  (by  the 
routine  SIGMA)  in  conjunction  with  a heat 
transfer  program  run  (subroutine  SIGMA  may 
be  executed  either  with  temperatures  sup- 
plied by  the  user  (NHTP=0)  or  with  tem- 
peratures supplied  by  the  heat  transfer 
calculations  of  the  ShTP  (NHTP=1))  (FIXED 
BIN)  ; and 


FILENAME  is  the  file  that  the  data  will  be  read 
from  (FILE). 


For  NHTP=0  (indicating  user  supplied  temperatures  and 
SIGMA  run  independently  of  the  SHTP)  the  following  data 
must  be  on  the  file  FILENAME: 


Two  comment  cards  followed  by  a group  of  data  cards 
arranged  in  increasing  order  of  radius.  These  cards  are 
of  two  types,  either  region  boundary  or  temperature 
(all  values  are  in  F- format). 

Region  boundary  cards  have  this  format: 


Columns 

45-54 

55-64 

Temperat  ure 
Columns 
1-15 
16-30 
31-40 


Value  and  Variable  Name 

Region  boundary  radius  (in.)  (XINIR) ; and 

Region  boundary  temperature  (°F)  . 

cards  have  the  following  format: 

Value  and  Variable  Name 

Radius  value  (in.)  (RADII)); 

Temperature  value  (°F)  ; and 

Mechanical  property  identification  number 
(MAT)  . 


The  first  and  last  cards  must  be  Region  boundary 
cards.  Region  boundary  and  temperature  cards  may  have 
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the  sane  value  of  radius  but  the  two  types  of  data  must 
be  on  separate  cards. 

For  NHTP= 1 (indicating  routine  SIGMA  being  used  in 
conjunction  with  SHTP  routines  that  will  calculate  tem- 
peratures) the  tabular  data  entries,  arranged  in  increa- 
sing order  of  Radius  values,  will  be  of  the  following 
two  types: 

Region  Boundary: 


Columns  Value  and  Variable 

40-54  Region  boundary  radius  (in)  (XINTR) . 


Temperature: 

Columns 

1-15 

16-25 

26-35 


Value  and  Variable 

Value  of  radius  (in.)  (RADII); 

Thermal  node  number  from  the  thermal 
network  (IND)  ; and 

Mechanical  property  identification  number 
(MAT)  . 


As  before,  the  first  and  last  entries  must  be 
region  boundary  cards.  Temperature  and  region  boundary 
cards  may  have  the  same  value  cf  radius  but  the  two 
types  of  information  must  not  appear  on  the  same  card. 


The  arangement  of  cards  just  described  is  designed 
to  allow  the  geometry  data  to  be  input  in  a way  that 
resembles  the  geometry  being  modeled;  i.e.,  the  cards 
defining  stress  region  boundaries  and  those  defining 
temperatures  at  various  radii  are  interleaved  in  increa- 
sing radial  order.  The  only  other  convention  is  that  the 
first  and  last  cards  must  be  stress  boundary  cards  even 
if  temperatures  are  to  be  defined  at  these  radii  as 
well.  Recall  that  region  boundary  cards  and  temperature 
cards  may  have  the  same  value  cf  radius,  if  desired. 
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Appendix  G 


UNIFORM  GEOMETRY  ROUTINES 


Subroutines  are  available  in  the  SHTP  library  that 
may  be  used  to  calculate  automatically  the  nodal 
volumes,  numbers,  and  interconnection  information  for 
uniform  thermal  grids.  The  routines  can  be  used  for 
one-  or  two-dimensional  grids  in  rectilinear  (Cartesian) 
or  cylindrical  coordinates. 


T NOD 


Subroutine  TWOD  is  used  for  describing  Cartesian 
nodal  networks  and  is  called  as  follows  (refer  to  Fig.  6 
for  symbol  conventions): 

CALL  TWOD  (ID# , Z,  XS  , NX,  YS  , NY,  IPROP,  FRSTNOD, 
SIDE1,  SIDE2,  SIDE3 , SIDE 4)  ; 


where: 

ID#  is  an  identification  number  for  the  two- 

dimensional  region  being  generated  (FIXED 
BIN)  ; 

Z is  the  depth  dimension  to  the  nodal  grid, 

the  z-direction  (ft)  (FLOAT  BIN); 

XS  is  a list  of  x-direction  thicknesses  to  be 

used  in  the  nodal  break-up  ( (*)  FLOAT  BIN) 
(if  XS  is  dimensioned  to  1,  then  the 
single  value  will  be  taken  as  the  total 
x-direction  thickness  and  the  variable  NX 
will  be  taken  as  the  number  of  evenly 
spaced  divisions  that  will  be  made  within 
the  x-direction  thickness)  ; 

NX  is  the  number  of  x-direction  divisions 

that  are  to  be  used  (FIXED  BIN)  ; 
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YS  is  the  list  of  y-direction  thicknesses  to 

be  used  in  the  nodal  break-up  ( (*)  FLOAT 
BIN)  (if  YS  is  dimensioned  to  1 then  the 
single  value  will  be  taken  to  indicate  the 
total  extent  of  the  region  in  the  y-direc- 
tion and  NY  will  be  understood  to  be  the 
number  of  equally  spaced  divisions  that 
will  be  made)  ; 


NY  is  the  number  of  y-direction  divisions 

that  are  to  be  used  (FIXED  BIN)  ; 

IPROP  is  the  thermal  property  identification 

number  for  the  single  material  that  will 
be  used  in  the  two  dimensional  region 
(FIXED  BIN)  ; 

FRSTNOD  is  the  node  number  of  the  first  node  in 
this  region  (FIXED  BIN)  , 

SIDE1-  are  flags,  either  0 or  1 that  determine 

SIDE4  the  type  of  thermal  path  length  that  will 

be  used  for  the  four  surfaces  of  the  2-D 
region  (FIXED  BIN). 


The  routine  TWOD  is  seen  to  have  two  modes  of 
operation,  one  to  allow  the  specification  of  an  ar- 
bitrary mesh  break-up  or  the  ether  to  allow  a uniform 
one.  The  choice  is  determined  by  the  dimensions  of  the 
two  lists  XS  and  YS.  All  length  values  are  to  be  in 
feet.  The  routine  prints  out  a table  of  results  that 
shows  node  number,  volume,  material  identification,  the 
number  of  interconnecting  nodes,  and  the  node  numbers 
and  thermal  path  lengths  for  the  interconnecting  nodes. 
The  output  is  annotated  as  being  from  the  TWOD  routine 
and  gives  the  region  identification  number. 


CYL2D 


The  CYL2D  routine  calculates  the  volumes  and  ap- 
propriate A/L  path  lengths  for  two-dimensional  cylin- 
drical regions  of  arbitrary  depth.  The  routine  is 
called  as  follows  (refer  to  Fig.  7 for  nomenclature) : 

CALL  CYL2D (ID#,  Z,  RADII,  NT,  ROUT,  THETAS,  NTH, 
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where : 

ID# 

Z 

RADII 

NR 

ROOT 

THETAS 

NTH 

IPROP 

PRSTNOD 

SIDE  1 , 
SIDE4 


IPROP,  FRSTNOD,  SIDE1,  SIDE2,  S1DE3, 
SIDE 4)  ; 


is  the  region  identification  number  (FIXED 
BIN)  ; 

is  the  region  depth  (rt)  (FLOAT  blN)  ; 

is  a list  of  radial  spacings  that  will 
determine  the  node  thicknesses,  measured 
from  the  outer  radius  inward  (ft)  (if 
RADII  is  dimensioned  to  1 then  the  single 
value  will  be  taken  to  be  the  total  radius 
thickness  and  NR  will  be  the  number  of 
evenly  spaced  divisions  that  will  be  used 
in  the  nodal  descriptions)  ( (*)  FLOAT 
BIN)  ; 

is  the  number  of  radial  divisions  (FIXED 
BIN)  ; 

is  the  cuter  radius  of  the  cylindrical 
region  (ft)  (FLOAT  BIN)  ; 

is  the  list  of  angular  (circumferential) 
divisions  that  will  be  used  to  describe 
the  nodal  network  (deg)  (if  THETAS  is 
dimensioned  to  only  1,  then  the  single 
value  will  be  taken  to  be  the  total 
angular  extent  of  the  region  and  NTH  will 
specify  the  number  of  equally  spaced 
divisions  that  will  be  used)  ( (*)  FLOAT 
BIN)  ; 

the  number  of  circumferential  divisions  in 
the  regions  (FIXED  EIN)  ; 

is  the  thermal  property  identification 
number  for  the  single  material  in  the 
region  (FIXED  BIN)  ; 

is  the  number  of  the  first  node  in  the 
region  (FIXED  BIN)  ; and 

are  flags  (either  0 or  1)  indicating  the 
type  of  thermal  path  calculations  that 
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will  be  made  at  the  four  region  surfaces 
(FIXED  BIN) . 

The  routine  prints  out  the  results  of  its  calculations 
in  the  form  of  a table  with  an  entry  for  each  node,  as 
follows:  node  number,  volume,  material  identification, 
the  number  of  interconnecting  nodes,  and  the  node  number 
and  thermal  path  lengths. 


CYL3D 


The  CYL3D  routine  is  used  much  the  same  as  CYL2D 
except  that  the  region  being  modeled  is  assumed  to  be 
three-dimensional;  i.e.,  the  depth  dimension,  Z,  is 
allowed  to  be  divided  into  sections.  The  calling  se- 
quence is: 


CALL  CYL3D  (NREG , Z,  NZ,  RADII,  NR,  ROUT,  THETAS, 
NTH,  IPROP,  FRSTNOD , SIDE1,  SIDE2,  SID 
SIDE 5 , SIDE6  , XTOFEET) ; 


where : 

Z is  a list  of  axial  (z-direction)  spacings 

to  be  used  in  defining  the  grid  (as  with 
the  other  lists  cf  spacings  this  list  may 
contain  NZ  numbers  that  will  specify  the 
z-direction  grid  spacings  or,  if  dimen- 
sioned to  1 , will  specify  the  total  axial 
extent  of  the  region,  and  NZ  will  be  taken 
to  be  the  required  number  of  evenly  spaced 
divisions; 


NZ  is  the  number  of  z direction  spacings; 

SIDES,  are  the  code  numbers  (either  0 or  1)  for 

SIDE 6 determining  the  type  of  connection  length 

to  be  used  for  the  additional  surfaces  (in 
reference  to  Fig.  7,  these  surfaces  are 
the  front  (SIDE5)  and  back  (SIDE6) ) ; and 

XTOFEET  is  a conversion  factor  for  obtaining  the 
units  of  feet  from  the  units  as  supplied 
for  Z,  ROUT,  and  RADII  (this  term  will  be 
multiplied  times  each  value  of  Z,  ROUT, 
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and  R AD II  and  the  results  will  be  assumed 
to  be  in  feet) . 

The  CYL3D  routine  will,  as  required,  store  connec- 
tion information  for  nodes  with  three  interconnections 
via  the  indirect  address  method  discussed  in  Appendix  F, 
READCP.  In  this  technique,  the  interconnection  matrix, 
XCAP , must  be  dimensioned  large  enough  to  allow  the 
appropriate  number  of  extra  entries  (refer  to  the 
description  of  subroutine  STORE,  Appendix  H)  . 

The  parameter  XTCFEET  represents  an  improvement 
over  CYL2D,  and  will  be  included  therein  upon  the  next 
revision  to  that  routine.  Further,  the  conversion 
implied  by  this  term  will  not  change  the  numbers  sup- 
plied as  arguments;  i.e. , the  conversions  will  be  stric- 
tly internal  to  the  routine  CY13D  and,  upon  return  from 
the  routine,  the  values  of  Z,  ROUT,  and  RADII  will  not 
have  been  altered. 

The  results  of  CYL3D's  calculations  will  be  printed 
in  a manner  similar  to  that  of  CYL2D  with  the  audition 
that  the  indirect  address  entries  to  XCAP  will  be 
included. 
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Appendix  H 

INITIALIZATION  ROUTINES 


In  order  to  execute  a heat  transfer  program  using 
the  SHTP  library,  certain  values  must  be  set  and  the 
proper  storage  must  be  allocated  for  certain  variables. 
Because  the  dynamic  storage  allocation  feature  of  the 
PL/I  language  has  been  widely  incorporated  into  the  SHTP 
library,  a specification  of  the  storage  requirements 
must  be  made.  Essentially,  the  amount  of  storage  needed 
for  variables  is  a direct  function  of  the  number  of 
nodes  in  the  thermal  network.  The  routines  SET  and 
STORE  accomplish  the  allocation  and  initialization  of 
many  of  the  needed  variables. 

STORE 


Subroutine  STORE  allocates  storage  according  to  its 
arguments  and  is  called  by  the  following  statement: 

CALL  STORE  ( #CA PS , XCAPLIM,  TYPRUN); 
where : 

♦CAPS  is  greater  than  or  equal  to  the  number  of 

thermal  nodes  in  the  current  problem 
(FIXED  BIN) ; 

XCAPLIM  is  the  number  of  locations  to  be  reserved 
in  the  XCAP  array;  note  that  XCAPLIM 
should  be  > #CAPS,  (FIXED  BIN)  (in 
problems  where  multiple  interconnections 
between  nodes  are  specified  via  indirect 
card  entries  to  the  READCP  routine, 

XCAPLIM  must  be  set  to  the  value  of  the 
largest  such  indirect  card  number  (c.f. 
READCP,  Appendix  F)); 

TYPRUN  indicates  the  type  of  run  being  made 
according  to  the  following  values: 

1 Normal  run,  transient  heat  trans- 
fer problems; 
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2 Steady  state  analyses,  iterative 
approach ; 

3 Special  run,  used  for  certain 
ablation  subroutines; 

4 Steady-state  analyses  via  a matrix 
solution  technique. 

The  TYPRUN  parameter  is  normally  set  to  1 but  must 
be  set  to  the  other  values  according  to  the  type  of 
analysis  desired.  Appendix  N cf  this  volume  explains 
the  usage  of  the  two  steady  state  analyses,  and  the 
other  code  (TYPRUN=3)  is  for  a program  described  in  Ret. 
4 and  will  not  be  applicable  to  any  of  the  discussions 
in  this  report.  All  cf  the  variables  allocated  by  STORE 
are  initialized  to  values  of  zero. 


SET 


Subroutine  SET  is  used  to  initialize  the  nodal 
temperatures  and  to  set  the  starting  and  ending  times 
for  the  transient  heat  transfer  solution.  The  routine 
is  called  as  follows: 

CALL  SET  (STAPTIM,  STOPTIM,  TINIT,  #SET,  5ETN0DES, 
SETTEMPS) ; 

where : 

STARTIM  is  the  starting  time  for  the  transient 
heat  transfer  calculations  (s)  (FLOAT 
BIN)  ; 

STOPTIN  is  the  final  time  for  the  transient  cal- 
culations (s)  (the  heat  transfer  program 
will  stop  when  this  time  is  reached) 

(FLOAT  BIN)  ; 

TINIT  is  the  temperature  that  all  of  the  thermal 

nodes  will  be  set  to  (°R)  (FLOAT  blN) ; 

•SET  is  the  number  of  nodes  that  will  be  set  to 

temperatures  other  than  the  initial  tem- 
perature (TINIT)  (FIXED  BIN) ; 
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SETNODES  is  the  list  of  node  numbers  that  will  be 
set  to  specific  initial  temperatures  ((*) 
FIXED  BIN)  ; and 

SETTER  PS  is  the  corresponding  list  of  specific 

initial  temperatures  for  the  list  SETNODES 
((*)  FLOAT  BIN). 

All  of  the  nodes  listed  in  SETNODES  mill  be  initialized 
to  the  corresponding  values  in  SETTEMPS.  Zero  is  a 
valid  value  for  #SET,  in  which  case  the  arguments  SET- 
NODES  and  SETTEMPS  may  be  null  variables  (dimensioned  to 
0). 


-171- 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAU«EL  MARYLAND 


Appendix  I 

BOUNDARY  CONDITION  ROUTINES 


After  the  physical  model  of  the  structure  being 
analyzed  is  established,  the  boundary  conditions  the 
body  is  subjected  to  must  be  specified  as  a function  of 
time.  There  is  provision  in  the  SHTP  library  to  allow 
specified  temperature  and  heat  flux  boundary  conditions 
as  well  as  the  effects  of  aerodynamic  heating.  These 
routines  are  discussed  below. 


FORCER 


Subroutine  FORCER  allows  for  the  specification  of 
temperature  histories  for  selected  thermal  nodes.  The 
routine  is  called  as  follows: 

CALL  FORCER  (NODES,  TEMPS,  TIMES); 
w here : 

NODES  is  the  list  of  thernal  node  numbers  that 

will  be  constrained  to  the  specified  tem- 
perature history  ( (*)  FIXED  BIN);  and 

TEMPS,  are  the  schedule  of  temperatures  and 

TIMES  times,  respectively,  that  will  be  applied 

to  the  thermal  nodes  ( (*)  FLOAT  DIN). 

The  TIMES  array  must  have  as  its  first  entry  the  uumber 
of  time- temperature  pairs  that  there  are.  Also,  a value 
of  zero  for  any  entry  to  NODES  will  cause  any  subsequent 
entries  not  to  be  forced  to  follow  the  specified  tem- 
perature history. 
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FLUX 


Subroutine  FLUX  is  used  tc  impose  a specified  heat 
flux  to  a list  of  thermal  nodes.  The  routine  is  called 
as  follows: 

CALL  FLUX  (N , NODES,  AREAS,  VIEWS,  ALPHA,  TEMP, 

QDOT  , TIME)  ; 

where: 

N is  the  number  of  nodes  in  the  list  NODES 

to  be  affected  by  the  heat  flux  (FIXED 
EIN)  ; 

NODES  is  the  list  of  thermal  nodes  that  will  be 

affected  ((*)  FIXED  BIN); 

AREAS  is  the  list  of  surface  areas  that  will  be 

exposed  to  the  heat  flux  (ft2)  in  one-to- 
one  correspondence  with  the  entries  in 
NODES,  ( (*)  FLOAT  BIN)  , 

VIEWS  is  the  list  of  view  factors  that  will  be 

applied  to  the  surfaces  listed  iu  AREAS 
((*)  FLOAT  BIN), 

ALPHA,  are  the  lists  of  absorptivities  and  tem- 

TEMP  peratures  ( °F)  that  describe  the  absorp- 

tive properties  of  the  surfaces  (each  (*) 
FLOAT  BIN)  ; and 

QDOT,  are  the  lists  of  heat  flux  and  time  values 

TIME  that  define  the  heat  flux  history  to  be 

applied  to  each  node  in  NODES  (Btu/ft2-h) 
(each  (*)  FLOAT  EIN). 

The  value  of  heating  rate  that  will  be  applied  to 
each  of  the  nodes  in  NODES  will  be  the  product  of  the 
following  terms:  the  surface  area  (AREAS) ; the  view 
factor  (VIEWS) ; the  absorptivity  of  the  surface,  evalua- 
ted at  the  surface  temperature;  and  the  heat  flux  at  the 
current  time.  Note  that  the  end  result  of  this  series 
of  multiplications  is  a heat  rate  term  expressed  in 
Btu/h.  With  appropriate  input  by  the  user,  the  flux 
routine  can  be  used  to  simulate  another  boundary  con- 
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dition,  namely  internal  heat  generation.  In  this  case 
the  QDOT  values  should  be  in  units  of  Btu/ft*-h  and  the 
AREAS  values  should  be  the  volumes  of  the  nodes  listed 
in  NODES.  The  net  result  in  the  heat  transfer  program 
is  the  same:  a heat  rate  (Btu/h)  is  added  to  the  appro- 
priate node. 

The  first  entries  to  the  independent  variables  TEMP 
and  TIME  must  be  the  number  of  total  data  entries  in  the 
lists  ALPHA  and  QDOT,  respectively.  This  is  the  stan- 
dard convention  used  for  the  READ  routines  and  the 
interpolation  routines  that  must  use  the  data. 


SAERO 


The  subroutine  SAERO  is  used  to  determine  the  sur- 
face temperature  at  any  point  on  a structure  exposed  to 
high  speed  fluid  flow.  The  heat  transfer  calculations 
use  the  Reynolds  analogy  and  the  reference  enthalpy 
method  of  Eckert  to  derive  heat  transfer  coefficients. 
(Refer  to  "Heat  Transfer  --  Convection"  in  Section  1 of 
Volume  1 for  a complete  discussion  of  the  techniques 
used) . It  will  suffice  at  this  point  to  say  that  the 
SAERO  routine  performs  a heat  balance  across  the  heated 
surface  and  numerically  determines  the  surface  tempera- 
ture. This  temperature  is  then  provided  to  the  rest  of 
the  heat  transfer  program  essentially  as  a boundary 
condition.  The  routine  is  called  as  follows  (Refer  to 
Fig.  9 for  nomenclature)  : 

CALL  SAERO  (ID,  #CAPS,  SNODES,  RNODES,  AREA,  REFLF.N, 
SHAPE,  D L,  RECRIT,  VIEW,  TRAdlD, 

FLOID,  T LIM , QLIM,  ITIIM)  ; 


is  the  identifying  number  for  the  case  in 
point  (FIXED  BIN) ; 

is  the  number  of  individual  surface  nodes 
that  are  in  this  aero  block  and  are  there- 
fore affected  by  the  same  body  station 
parameters  (FIXED  BIN)  ; 

are,  respectively,  the  lists  of  surface 


where: 

ID 

iCAPS 

SNODES, 
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RHODES 

AREA 

REFLEN 

SHAPE 


D L 


RECRIT 


VIEW 


TRAJID 


node  numbers  and  connecting  capacitor 
numbers  used  in  modeling  the  structure 
( (*)  FIXED  BIN)  , 

is  a list  of  surface  areas  (ft*)  for  each 
of  the  SNODES  ((*)  FLOAT  BIN); 

is  a list  of  reference  lengths  (used  in 
the  Reynolds  number  calculation)  for  each 
of  the  SNODES  (ft)  ((*)  FLOAT  BIN)  ; 

is  a code  for  determining  the  kind  of 
aerodynamic  heating  calculations  that  will 
be  done.  The  following  codes  apply: 

1 flat  plate  correlations, 

2 cone  flow  correlations, 

3 pipe  flow  correlations, 

4 stagnation  point  (two-dimension- 
al) , and 

5 stagnation  point  (three-dimen- 
sional) (FIXED  BIN)  ; 

is  a list  of  values  used  in  determining 
the  heating  rate  for  pipe  flow  and  is  the 
ratio  of  the  diameter  to  the  length-down- 
stream- from-opening  for  each  of  the 
SNODES;  if  SHAPE  is  not  equal  to  3,  then 
these  values  have  no  effect  ((*)  FLOAT 
BIN)  (when  considering  pipe  flow  it  should 
be  remembered  that  REFLEN  must  be  set 
equal  to  the  hydraulic  diameter)  ; 

is  the  value  of  the  transition  Reynolds 
number  between  laminar  and  turbulent  flow 
(FLOAT  BIN)  ; 

is  the  single  view  factor  to  be  applied  to 
each  SNODE  equally  for  the  input  radia- 
tion, Q RAD , described  above  for  READTH, 
(FLOAT  BIN)  , 

is  an  identifying  number  that  numerically 
matches  an  identification  number  tor  some 
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FLOID 


TLIM  , 
QLIM  , 
ITLIM 


READTH  statement  that  will  contain  values 
of  velocity  and  altitude  versus  tine 
(FIXED  BIN)  , 

is  an  identifying  number  that  corresponds 
numerically  with  an  identification  number 
for  some  READFL  statement  wherein  the 
appropriate  values  of  local  to  freestream 
Mach  No.  ratio,  local  to  freestream  pres- 
sure ratio,  and  velocity  gradient  factor 
were  read  in  versus  freestream  Mach  No. 
(FIXED  BIN)  ; 

are  the  parameters  for  judging  the  con- 
vergence of  the  two  iteration  processes 
carried  out  by  SAERO  (in  each  of  these 
iterations  the  calculated  temperature  must 
be  within  TLIM  of  the  true  temperature 
solution  and  the  resulting  energy  balance 
must  have  a residual  of  not  more  than 
QLIM.  Further,  no  more  than  ITLIM  itera- 
tions will  be  allowed.  If  any  of  the 
convergence  criteria  are  not  met,  a per- 
tinent message  is  printed  and  the  execu- 
tion will  continue). 
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Appendix  J 

HEAT  TRANSFER  CALCULATIONS 


The  heat  transfer  program  calculates  the  transient 
thermal  response  of  a heated  structure  via  a forward 
■arching,  finite  difference  technique.  The  details  of 
the  finite  difference  solutions  can  be  found  in  Section 
2 of  Volume  1.  The  subroutines  that  perform  this  task 
are  described  in  the  following  paragraphs. 


CAPCON 


The  routine  CAPCON  uses  the  nodal  network  infor- 
mation stored  in  the  external  array  XCAP  to  calculate 
the  temperature  dependent  values  of  thermal  capacity 
( Btu/f  t3  -°R)  for  each  node  and  the  sum  of  all  the  con- 
ductances for  each  node  (dtu/°R-h).  These  values  are  the 
essential  coefficients  in  the  finite  difference  equa- 
tions and  are  used  later  by  subroutine  STEP.  CAPCON  is 
called  as  follows: 

CALL  CAPCON; 

There  are  no  arguments  to  CAPCON  because  all  the  essen- 
tial parameters  (material  identification  numbers,  con- 
duction paths,  volumes,  etc.)  have  been  organized 
previously  and  stored  in  the  appropriate  external  varia- 
bles via  calls  to  FEADRK,  FEADCP,  CYL2D,  THOD,  etc. 


STEP 


Subroutine  STEP  has  two  essential  functions:  to 
determine  the  stability  time  step  and  calculate  the 
change  in  nodal  temperatures  from  one  time  point  to  the 
next.  Refer  to  "Transient  Heat  Conduction"  in  Section  2 
of  Volume  1 for  a discussion  of  the  finite  difference 
equations  and  the  stability  criterion.  Subroutine  STEP 
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is  called  as  follows: 

CALL  STEP  ( MINSTEP,  HAXSTEP) ; 

where: 

MINSTEP  is  the  smallest  permissible  time  step  for 
the  calculations  (FLOAT  BIN) , and 

HAXSTEP  is  the  maximum  permissible  time  step  that 
will  be  used  (FLOAT  BIN)  . 

If  a calculated  time  step  is  less  than  MINSTEP  the 
program  will  terminate  with  an  error  message.  It  a cal- 
culated time  step  is  larger  than  MAXSTEP,  then  the  time 
step  is  set  to  HAXSTEP  and  execution  continues. 
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Appendix  K 

MATERIAL  RESPONSE  ROUTINES 


There  is  a collection  of  subroutines  that  have  been 
devised  to  couple  with  the  SHTE  so  as  to  simultaneously 
calculate  selected  secondary  material  responses  of  the 
structures  being  analyzed.  These  secondary  responses 
include  thermal  stresses  for  cylindrical  geometries, 
boresight  errors  of  radar  through  a monolithic  wall 
radome,  and  stresses  in  a radome  due  to  aerodynamic 
load. 


SIGMA,  SIGMET 


These  two  routines  calculate  thermal  stresses  for 
cylindrical  geometries  according  to  a finite  difference 
technique  developed  by  Rivello  (Ref.  3,  and  "Tnermal 
Stress”  in  Section  1 of  Volume  1).  SIGMET  differs  from 
SIGMA  in  that  it  performs  the  additional  calculation  of 
the  electrical  thickness  change  of  the  wall  segment. 

The  thermal  stresses  are  calculated  by  making  a thick- 
wall  cylinder  assumption.  The  cylinder  can  be  made  up 
of  several  different  materials  and  may  be  assumed  either 
fully  restrained  from  axial  expansion  or  free  to  expand 
axially.  The  electrical  thickness  calculations  are 
based  on  calculations  of  the  wall's  physical  thickness 
change  and  the  temperature  coefficient  of  the  dielectric 
constant  versus  temperature.  The  complete  calculation 
of  boresight  errors  is  accomplished  with  the  help  of 
subroutine  MOBSER,  whose  use  will  also  be  described. 


SIGMA 


The  routine  SIGMA  calculates  thermal  stresses  in 
cylinders  in  the  hoop,  axial,  and  radial  directions.  It 
reguires  lists  of  temperature  versus  radius  ana  elastic 
property  data  versus  temperature  to  do  so.  These  data 
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can  be  read  in  and  made  available  to  SIGMA  via  routines 
READSG  and  READFH.  The  calling  sequence  is: 

CALL  SIGMA  (ID,  N REG , B,  F,  IND,  MAT,  TEMP,  XMOD , 
DELTAL,  XNU , EO , FN,  TNOSTRS,  ENDR, 
PLTCODE)  ; 

where : 

ID  is  an  identifying  number  (FIXED  BIN); 

N REG  is  the  number  of  regions  or  subdivisions 

that  will  be  made  through  the  wall  for  the 
finite  difference  solution  (FIXED  BIN)  ; 

b is  the  list  of  NREG+1  radius  values  that 

define  the  region  boundaries  (in.)  ((*) 

FLOAT  BIN)  ; 

R is  the  list  of  radius  values  at  which 

temperatures  will  be  specified  (in.)  ( (*) 

FLOAT  BIN)  ; 

IND  is  the  list  of  thermal  network  capacitor 

numbers  that  will  be  assumed  to  exist  at 
each  respective  radial  position  indicated 
in  R,  above  ((*)  FIXED  BIN); 

HAT  is  the  list  of  material  indexes  (property 

code  numbers)  that  will  identify  the 
material  properties  that  exist  at  each 
radial  position,  R,  above  (these  indexes 
are  keyed  to  identification  numbers  used 
by  READ  RK  and  READFH)  ( («■)  FIXED  BIN)  ; 

TEMP,  are  two-dimensional  arrays  of  temperature 

XMOD,  (°F)  , elastic  modulus  (psi)  , and  free 

DELTAL  thermal  expansion  (in. /in.)  for  the 

materials  used  in  the  cylindrical  section. 
The  first  subscript  position  of  each  of 
these  arrays  contains  the  material  iden- 
tification number  (specified  by  MAT)  and 
the  second  subscript  position  contains  the 
particular  property  data  (each  (*,*)  FLOAT 
BIN)  ; 

XNU  is  the  value  of  Poisson's  ratio  that  will 

be  used  in  the  calculation  (FLOAT  BIN)  ; 
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PO,  PN  are  the  internal  and  external  pressures, 

respectively,  that  may  exist  on  the  cylin- 
der (psi)  (FLOAT  BIN)  ; 

TNOSTBS  is  the  temperature  at  which  the  cylinder 
will  be  assumed  stress  free  (°F)  (FLOAT 
BIN)  ; 

ENDR  is  a flag  indicating  the  assumed  end 

restraint  for  the  cylinder;  e.g., 

1 rigidly  restrained  ends,  and 

0 ends  free  tc  expand  (FIXED  BIN)  ; 
and 

PLTCODE  is  a pseudo-rode  numner  that,  if  nonzero, 
will  be  used  to  store  the  maximum  tensile 
stress  that  exists  in  the  wall  at  each 
time  the  data  are  stored  (this  number  is 
used  as  an  index  to  the  array  of  thermal 
node  temperatures  so  that  if  plots  are 
used  (via  subroutine  PLOTTER)  , thermal 
stress  histories  may  oe  plotted,  if 
desired)  (FIXED  BIN), 

The  routine  SIGH  A may  or  may  not  be  used  with  the 
SKTP.  To  insure  the  proper  execution  in  either  case, 
routine  READSG  should  be  used  to  set  up  the  proper  data 
items.  When  used  with  the  SHTP,  SIGMA  will  execute  only 
at  the  times  listed  in  the  arguments  to  subroutine 
WRITE. 


SIGMET 


Subroutine  SIGMET  is  the  same  as  subroutine  SIGMA 
except  for  the  additional  calculation  of  the  electromag- 
netic (EM)  wave  thickness  change  of  the  wall  segment. 
This  calculation  regures  one  additional  piece  of  infor- 
mation, namely  the  incidence  angle  of  the  EM  wave.  This 
is  input  in  the  argument  list  as  the  variable  named  PHI 
in  the  example  below  of  the  calling  sequence: 
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CALL  SIGNET (ID  , NREG , B,  R,  IND  , HAT,  TEMP,  XMOD, 

DELTA,  XNU,  EO,  PN,  TNOSTRS,  ENDR,  PHI, 
PLTCODE)  ; 

where  all  the  variables  have  the  same  precise  defini- 
tions as  with  SIGMA  except: 

PHI  is  the  EM  wave  incidence  angle  (deg) 

(FLOAT  BIN) . 

The  calculation  of  electrical  thickness  change  is  essen- 
tial for  the  monolithic  boresight  error  slope  calcula- 
tions performed  in  subroutine  MOBSEE,  described  next. 

All  of  the  comments  pertinent  to  routine  SIGMA  apply  to 
routine  SIGMET  with  the  addition  that  the  material 
property  codes  now  apply  to  the  thermal  properties  (read 
in  by  KEADEK)  as  well  as  the  mechanical  properties  (read 
in  by  EEADFH)  . 


MOBSEE 


Subroutine  MOBSEE  uses  data  from  routine  SIGMET  and 
from  its  own  arguments  to  calculate  a boresight  error 
slope  (deg/deg)  for  a radar  beam  through  a monolithic 
radome  wall.  The  calculation  is  based  on  an  empirical 
relationship  among  the  parameters  mentioned  and  applies 
to  von  Karman  shapes.  The  routine  is  called  as  follows: 

CALL  MOBSEE  (ID,  FREQ,  SEP,  PLOTID)  ; 

where : 

ID  is  a case  identifying  number  and  keys  the 

routine  back  to  a previous  CALL  to  SIGMET 
wherein  the  electrical  thickness  change 
was  calculated  (FIXED  BIN)  ; 

FREQ  is  the  radar  frequency  of  the  incident 

beam  (H7)  (FLOAT  BIN)  ; 

SEP  is  the  effective  phase  center  separation 

of  the  radar  antenna  aperture  (in.)  (FLOAT 
BIN)  ; and 
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PLOTID  a psuedo-node  number  used  to  store  the 

boresight  error  slope  and  the  electrical 
thickness  change  values  at  each  time 
point.  This  number  and  the  next  highest 
integer  (-.e.,  PLOTID+1)  are  to  index  the 
array  of  thermal  node  temperatures  and 
therefore  must  not  be  the  numbers  of 
actual  network  nodes.  A value  of  zero  for 
PLOTID  will  cause  no  such  storage  to  take 
place.  Later,  if  the  routine  PLOTTER  is 
used  to  store  or  plot  the  nodal  temperat- 
ures, these  indexes  can  be  used  to 
generate  boresight  error  slope  and  elec- 
trical thickness  change  histories. 

The  variable  ID  joins  the  calculation  of  electrical 
thickness  change  to  the  other  parameters  and  routine 
SIGNET  must  be  called  before  routine  MOBSER. 


AERLOAD 


Subroutine  AERLOAD  is  used  principally  for  radome 
analyses,  wherein  the  stresses  due  to  aerodynamic  pres- 
sure and  maneuver  loads  at  the  missile-radome  attachment 
are  required.  The  routine  uses  the  radome  and  missile 
aerodynamic  force  coefficients  read  by  READAFC,  and  is 
keyed  via  the  argument  variable  TRAJID  to  a trajectory 
definition  read  in  by  a READTM  statement.  The  routine 
is  called  as  follows: 

CALL  AERLOAD  (ID,  L,  RB,  T,  CKSTA,  RHO  , WM,  YDPDOT, 
ALTHDOT,  ALTHCDT,  TRAJID,  PLOTID) ; 

where: 

ID  is  an  identifying  number  (FIXED  BIN)  ; 

L is  the  length  of  the  radome  from  tip  to 

attachment  point  (in.)  (FLOAT  BIN)  ; 

RB  is  the  radius  of  the  radcme  at  the  attach- 

ment point  (in.)  (FLOAT  BIN); 

T is  the  thickness  of  the  radome  wail  at  the 

attachment  point  (in.)  , (FLOAT  BIN)  ; 
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CHSTA  is  the  missile  body  station  that  locates 

the  missile’s  center  of  mass  (in.),  (FLOAT 
BIN)  ; 

EHO  is  the  density  of  the  radome  material 

(lbm/f t 3)  (FLOAT  BIN); 

WM  is  the  weight  of  the  missile  (lbf)  (FLOAT 

BIN)  ; 


YDDOT  is  the  maximum  possible  maneuver  accelera- 

tion normal  to  the  axis  of  the  missile 
(ft/s*)  (FLOAT  BIN)  ; 


ALTHDOT  is  the  maximum  possible  rate  of  change  of 
angle  of  attack  for  the  missile  (deg/s) 
(FLOAT  BIN) ; 


ALTHDDT  is  the  maximum  possible  acceleration  of 
angle  of  attack  for  the  missile  (deg/s2) 
(FLOAT  BIN)  ; 


TRAJID  is  an  identifying  number  for  the  trajec- 
tory information  neeaed  by  AERLOAD  (this 
number  must  match  an  identification  number 
for  a call  to  routine  READTM  elsewhere  in 
the  program)  (FIXED  BIN)  ; and 

PLOTID  is  a psuedo-node  number  used  for  storing 
the  maximum  te  stress  calculated  by  AEF- 
LOAD  at  each  time  point  (this  number  must 
not  correspond  to  an  active  thermal  node 
in  the  network;  a value  of  zero  will  cause 
no  storage  to  take  place). 


Part  of  the  data  that  will  be  expected  in  the 
trajectory  tables  is  the  maximum  permissible  angle  of 
attack.  These  data  are  allowed  to  be  variable  with  time 
and  is  therefore  implicitly  dependent  on  altitude  (c.f. 
Appendix  F,  subroutine  READTH)  . If  the  specified 
lateral  acceleration,  YDDOT  above,  requires  an  angle  of 
attack  greater  than  that  allowed,  then  the  maximum  angle 
of  attack  is  used  and  the  resultant  lateral  acceleration 
is  accepted.  The  AERLOAD  routine  executes  only  at  the 
times  listed  in  the  argument  to  subroutine  WRITE. 
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Appendix  L 
OUTPUT  ROUTINES 


There  are  two  kinds  of  output  available  from  the 
SHTP  --  printed  results  and  CalComp  plots.  The  printed 
results  come  mainly  from  the  individual  subroutines  that 
perform  the  respective  computations.  The  nodal  network 
temperatures  are  obtained  from  the  subroutine  WRITE 
which  keys  other  subroutines  to  print  their  results. 
Subroutine  WRITE  has  two  forms  --  one  that  prints  a well 
formatted,  compact  output  of  the  temperatures  and 
another  that  can  be  used  for  debugging  purposes  that 
prints  the  temperatures  and  other  important  program 
variables.  The  CalComp  plotting  routines  are  set  up  so 
that  plots  of  temperatures  versus  time  can  be  made 
either  during  a program  run  or  later,  after  the  run. 

The  specific  details  of  these  capabilities  are  discussed 
under  each  routine  name  below. 


WRITE 


Subroutine  WRITE  is  used  to  indicate  when  the 
printed  output  should  be  made.  The  calling  sequence  is: 

CALL  WRITE  (PRNTIHS)  ; 
where: 

PRNTHS  is  a list  of  times  (s)  at  which  the  prin- 
ted output  will  be  made  ((*)  FLOAT  BIN). 

If  PRNTIHS  is  dimensioned  to  1 then  the  single 
value  placed  therein  will  be  taken  as  an  increment  to  be 
successively  applied  to  the  starting  time  of  the  program 
to  achieve  an  evenly  spaced  seguence  of  printouts. 

If  the  value  of  PRNTIHS  (1)  is  negative,  then  the 
printed  output  will  be  made  after  every  time  step  in  the 
program.  This  option  should  be  carefully  used  only  for 
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debugging  because  considerable  unnecessary  output  can 
result. 

Within  routine  WRITE  is  a one-bit  flag  called 
PRINTIT  which  is  accessible  by  every  routine  in  the  SHTP 
that  generates  printed  output.  Each  time  routine  WRITE 
prints  its  own  output,  the  PRINTIT  bit  is  set  to  indi- 
cate to  the  other  routines  to  print  as  well.  In  this 
way,  the  schedule  of  ti*es  in  the  argument  PRINTIMS 
regulates  the  execution  and/or  printout  of  the  other 
computational  routines. 

Subroutine  WRITE  can  be  used  to  obtain  the  nodal 
temperature  data  in  a different  format  than  the  one 
described  and  shown  in  Appendix  C.  In  addition  to  the 
node  numbers  and  temperatures,  program  variables  C 
(nodal  capacitance) , H (nodal  conductance) , and  Q (net 
nodal  heat  flow)  are  printed.  Additionally,  if  there  are 
composite  conduction  paths  specified  in  the  program  then 
the  information  pertinent  to  these  connections  is  also 
printed.  The  routine  uses  more  paper  than  the  standard 
WRITE  routine,  but  can  be  helpful  in  debugging.  To 
obtain  this  alternate  output  feature,  simply  supply  the 
first  value  of  PRNTIHS  as  999. 


LIHIT 


Subroutine  LIHIT  was  devised  strictly  for  the 
radome  limitations  work  and  is  used  to  save  pertinent 
data  from  an  URLIM  run  for  plotting  by  the  routine 
LIHPRNT,  described  below.  The  calling  sequence  for 
LIHIT  is: 

CALL  LIHIT (ID,  FHHAT,  RKHAT , TIHEID,  HELTNODE, 
TITCARD)  ; 

where: 

ID  is  an  identifying  number  and  will  later 

show  up  as  the  limit  case  number  (FIXED 
BIN)  ; 

FHHAT  is  the  mechanical  property  data  number  for 
the  radcme  material  and  was  assigned  via 
an  earlier  call  to  routine  READFH  (FIXED 
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RKMAT 


TIHEID 


SIGID 


LODID 


HELTNODS 


TITCARD 


BIN)  ; 

is  the  thermal  property  ID  number  for  the 
radome  material  and  was  assigned  via  an 
earlier  call  to  routine  READRK  (FIXED 
BIN)  ; 

is  a list  of  trajectory  identification 
numbers  that  will  serve  to  define  the 
trajectories  over  which  the  limit  case 
will  be  performed;  these  identification 
numbers  must  have  been  previously  assigned 
via  calls  to  routine  READTM  ((*)  FIXED 
BIN)  ; 

is  a list  of  thermal  stress  case  numbers 
from  subroutine  SIGMET  that  will  be  used 
to  determine  thermal  stress  failure  points 
along  the  trajectories  listed  in  TIMEID  (a 
one-to-one  correspondence  of  the  values 
listed  in  SIGID  and  TIMEID  will  be  made  to 
determine  which  stress  case  goes  with 
which  trajectory;  the  thermal  stress 
identification  numbers  will  be  used  to 
identify  boresight  error  case  numbers  as 
well  since  the  electrical  thickness  change 
calculations  are  performed  by  the  same 
routine  — SIGMET)  ((*)  FIXED  BIN); 

is  a similar  one-to-one  list  of  aerodyna- 
mic load  calculation  indexes  that  identify 
previous  calls  to  routine  AERLOAD  ( (*) 
FIXED  BIN)  ; 

is  a two-dimensional  array  of  thermal  node 
numbers  that  will  select  the  nodes  to  be 
monitored  for  melting.  The  first  sub- 
script position  corresponds  one-to-one 
with  each  trajectory  in  TIMEID  and  the 
second  subscript  position  lists  all  of  the 
nodes  susceptible  to  melting  ((*,*)  FIXED 
BIN) ; and 

is  an  80-character  string  that  will  be 
used  to  identify  the  limits  plot  that  is 
made  by  LIM PENT  (CHARACTER  (80)). 
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Summary  Note 


For  each  trajectory  listed  in  TIMEID  above,  there 
will  be:  a stress  case  identification  number,  a load 
case  identification  number,  and  a list  of  node  numbers 
subject  to  melting.  The  boresight  error  case  numbers 
will  be  assumed  the  same  as  the  stress  case  numbers. 

Subroutine  LIMIT  uses  the  identification  number 
information  supplied  via  its  arguments  to  monitor  the 
level  of  the  radome  material's  responses  as  compared  to 
maximum  allowable  values.  In  the  case  of  thermal 
stresses  the  maximum  tensile  thermal  stress  for  each 
stress  case  is  stored  in  the  EXTERNAL  variable  MAXSTRS 
which  is  an  (*)  FLOAT  BIN  variable  and  is  indexed  with 
the  identification  numbers  in  either  routine  SIGMA  or 
SIGMET  (i.e.,  the  variable  SIGID  above).  The  LIMIT 
routine  compares  the  values  of  MAXSTES  (SIGID)  with  the 
maximum  permissible  tensile  stress  stored  by  the  user  in 
the  external  variable  MOR . MOP  is  an  (*)  FLOAT  BIN 
array  indexed  with  the  identification  numbers  from 
routine  READFH  (i.e.,  variable  FHMAT  above). 

The  boresight  error  rate  values  are  stored  in  the 
external  variable  BSE F which  is  an  (*)  FLOAT  BIN  array 
indexed  with  the  identification  numbers  from  routine 
SIGMET  (i.e.,  variable  SIGID  above).  The  values  of  bore 
sight  error  slope  that  are  calculated  during  the  program 
execution  are  compared  to  a value  stored  by  the  user  in 
the  EXTERNAL  variable  MAXBSER.  HAXBSEB  is  an  (*)  FLOAT 
BIN  variable  that  is  indexed  by  the  material  identifica- 
tion numbers  of  routine  READFK  (i.e.,  the  variable  RKMAT 
above)  . 

In  the  case  of  aerodynamic  and  maneuver  loads  at 
the  attachment  point  in  the  radome,  the  EXTERNAL 
variable  AERSIG  ((*,2)  FLOAT  BIN)  keeps  track  of  the  two 
principle  stresses,  one  at  the  upper  section  and  one  at 
the  lower  section  of  the  radome  (c.f.  Fig.  10).  The 
first  subscript  position  of  AERSIG  is  indexed  by  the 
identification  numbers  in  routine  AEFLOAD  (i.e., 
variable  LODID  above.) 

In  the  case  of  melting,  the  thermal  nodes  listed  in 
HELTNODS  are  monitored  with  respect  to  the  last  tern- 


188 


THE  JOHNS  HOPKINS  UNIVERSITY 

applied  physics  laboratory 

laurel  Maryland 


perature  entry  in  the  thermal  property  tables  read  in  by 
the  appropriate  READS K statement.  In  other  words,  the 
last  temperature  supplied  in  each  of  the  thermal  proper- 
ty tables  is  taken  to  be  the  melt  temperature  for  that 
material. 

The  LIMIT  routine  monitors  the  specified  cases  for 
each  trajectory  and  makes  note  of  the  time  at  which  the 
maximum  allowable  value  is  exceeded.  The  program  per- 
forms this  monitoring  in  two  modes  --  in  the  first  mode 
as  soon  as  any  one  of  the  four  limits  is  reached,  the 
time  is  noted  and  further  monitoring  along  that  trajec- 
tory is  halted;  in  the  second  mode,  all  of  the  limits 
are  monitored  all  of  the  time  up  to  the  end  of  the 
trajectory.  The  mode  selection  is  accomplished  by 
specifying  a control  value  for  the  external  variable 
MINRUN:  0 implies  all  the  limits  all  the  time,  and  1 
implies  only  the  first  limit. 

At  the  end  of  the  heat  transfer  program  run  (when 
the  time  variable  has  passed  the  STOPTIM  value  in 
subroutine  SET)  , LIMIT  writes  its  accumulated  limits 
data  on  two  files.  The  first  is  called  LIMPLT  and  is  a 
sequential  record  file  used  to  preserve  the  numerical 
results  for  later  use.  The  other  file  is  SYSPRINT  and 
the  data  is  presented  in  graphical  form.  This  later 
printed  output  is  generated  via  an  internal  call  to 
routine  LIHPRNT.  The  numerical  file  LIMPLT  must  have  a 
data  definition  card  describing  it,  in  the  JCL  tor  the 
execution  of  the  job.  Figure  6 contains  an  example  of  a 
suitable  JCL  card.  The  file  is  a record  file;  i.e.,  the 
data  are  written  in  a machine  format  that  has  not  been 
converted  to  a character  format.  The  data  organization 
and  a discussion  of  hew  to  generate  limits  plots  from 
the  saved  data  after  an  URLIM  run  are  given  in  Appendix 
0.  Examples  of  the  printer  plots  generated  by  LIHPRNT 
are  given  in  Figures  4 and  5. 

PLOTTER 


The  PLOTTER  routine  is  used  for  saving  and  plotting 
the  nodal  temperature  data  of  an  SHTP  run.  The  nodal 
temperatures  are  plotted  versus  time  and  can  be  arranged 
in  any  format,  drawn  to  any  scale  and  annotated  with  any 
phrase.  The  routine  is  designed  so  that  plots  may  be 
made  concurrent  with  a run  of  the  SHTP  or  separately, 
after  an  SHTP  run.  This  latter  feature  allows  for  plot 
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selection  after  the  printed  results  of  an  analysis  have 
been  obtained  and  reviewed.  The  routine  is  general  and 
allows  a degree  of  flexibility. 

The  calling  sequence  is: 

CALL  PLOTTER  (TIMES,  DDNAME,  FORMAT,  TMPSCAL , T1M- 
SCAL,  TMPINIT,  TIMINIT,  TITLE)  ; 


where: 

TIMES  is  a list  of  tines  (s)  during  the  run  of  a 

main  program  at  which  the  variables  will 
be  plotted  (if  it  is  required  to  have  the 
values  plotted  at  every  tine  step  in  the 
SHTP  then  the  value  of  TIMES  (1)  should  be 
negative;  if  the  plotting  frequency  is  to 
be  at  regular  intervals  of  time,  then 
TIMES  should  be  dimensioned  only  to  1 and 
the  value  of  the  plotting  interval  desired 
placed  in  TIMES  (1)  ) ((*)  FLOAT  BIN); 

DDNAME  is  a character  string  variable  of  type 

CHA RACT ER (7)  which  names  a DD  card  in  the 
JCL  of  the  job  being  run.  This  DD  card 
describes  the  file  onto  which  the  values 
to  be  plotted  will  be  written. 

If  the  plot  is  being  made  concurrent  with  the  run- 
ning of  the  SHTP  then  the  file  described  on  the  DD  card 
nay  be  a temporary  file  and  deleted  after  the  job  is 
completed.  In  general,  however,  the  plot  file  may  be 
saved  and  used  again  to  plot  other  variables  if  desired. 
The  sample  program  in  Appendix  A uses  a DD  card  named 
PLOT  to  identify  the  temperature  file.  For  a temporary 
file  the  DD  card  should  be: 

//DDNAME  DD  UNIT=TEMP, 

//  SPACE=  (3500,  (10, 10)  , RLSE)  , 

//  DISP=  (NEW, DELETE) 

To  save  the  plot  file  use: 

//DDNAME  DD  DSN*group. username. xxxx, 

//  UNIT-SAVE, 

//  DISP= (NEW,CATLG)  , 

//  SPACE=  (3500,  (10,10)  , RLSE) 
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FORMAT 


TMPSCAL 


is  a two-dimensional  array  of  integers 
((*,*)  FIXED  BIN)  that  indicates  what 
variables  are  to  be  plotted  on  which 
plots.  The  extent  of  the  first  dimension 
of  FORMAT  is  the  number  of  different  plots 
being  made,  and  the  extent  of  the  second 
dimension  is  the  maximum  number  of  varia- 
bles plotted  on  any  one  plot  (i.e.,  if 
FORMAT  was  dimensioned  FORMAT  (x,y),  there 
would  be  x different  plots  and  a maximum 
of  y variables  plotted  on  any  one  plot.) 
FORMAT  then  contains  x lists  (rows)  of 
numbers,  each  list  containing  y numbers. 
The  actual  numerical  entries  to  FORMAT 
specify  node  or  capacitor  numbers  to  be 
plotted;  i.e.,  the  numbers  are  used  to 
reference  the  temperature  vector  I in  the 
SHTP.  An  entry  of  zero  in  FORMAT  causes 
no  curve  to  be  drawn.  Furthermore,  an 
entry  of  zero  at  the  start  of  any  row  in 
FORMAT  will  cause  no  plotting  for  the 
entire  row.  if  an  entry  to  FORMAT  is 
negative,  then  a plotting  symbol  will  be 
plotted  with  each  data  point  on  that 
curve.  If  ABS  (FORMAT  (x,y))  = 999,  then  a 
2-  by  2-in.  grid  will  be  drawn  by  the 
routine.  If  this  is  used  with  plain  plot- 
ting paper,  then  the  resultant  plot  is 
suitable  for  formal  presentation.  An 
example  is  shown  in  Fig.  L- 1 . If  FORMAT 
(1,1)  = 0,  then  no  plotting  at  all  will  be 
done  and  the  plotting  values  will  be 
stored  on  the  file  DDNAME; 

is  a list  of  scaling  values  to  be  used  on 
the  vertical  axis  of  each  plot.  There 
should  be,  therefore,  at  least  x entries 
to  TMPSCAL.  The  numerical  entries  specify 
the  number  of  data  units  (ft,  deg,  psi , 
etc.)  per  inch  of  plotting  paper  on  the 
vertical  axis.  TMPSCAL  has  the  attributes 
(*)  FLOAT  BIN.  An  additional  feature  of 
the  plotting  capability  is  affected 
through  the  TMPSCAL  variable;  when  FORMAT 
(1,  1)  #0  and  any  entry  to  TMPSCAL  is 
negative,  then  a list  of  normalization 
factors  will  be  read  in  from  the  SYSIN 
file  in  list  format.  These  values  will  be 
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divided  into  each  value  for  each  curve  on 
the  plot  that  the  negative  value  of 
TMPSCAL  would  have  been  associated  with. 
Further,  the  scale  value  used  along  the 
vertical  axis  of  the  plot  will  be  in 
accordance  with  the  absolute  value  of  the 
value  of  TMPSCAL.  In  other  words,  if  a 
plot  (numbered  n)  is  to  be  made  and  the 
corresponding  value  for  the  vertical  scale 
(TMPSCAL(n))  is  negative,  then  the  SYSIN 
file  is  read  for  y numbers  in  LIST  format. 
Each  of  these  y numbers  will  be  used  to 
normalize  the  respective  thermal  node 
values  when  plotted.  The  ith  curve  on 
plot  number  n is  the  set  of  values  T (FOR- 
MAT (n,i))  divided  by  the  ith  normalization 
factor  read  in  because  of  the  negative 
value  of  TMPSCAL  (n)  ; 

TIMSCAL  is  a constant  value  specifying  the  scale 

to  be  used  for  the  horizontal  (time)  axis. 
It  is  a constant  and  therefore  the  same 
for  each  plot.  The  value  is  the  number  of 
seconds  per  inch  of  plotting  paper  to  be 
used  for  the  time  axis  (FLOAT  BIN)  ; 

TMPINIT  is  a list  of  initial  values  for  the  begin- 
ning of  each  of  the  vertical  axes.  There 
must  be  at  least  x entries  to  TMPINIT, 
which  has  attributes  (*)  FLOAT  BIN); 

TIMINIT  is  a constant  value  specifying  the  initial 
value  to  be  used  as  the  origin  for  the 
horizontal  (time)  axis  (FLOAT  BIN)  ; 

TITLE  is  a character  variable  into  which  may  be 
stored  commentary  phrases  for  labeling  and 
identifying  the  plots.  TITLE  is  an  array 
dimensioned  to  2 and  each  subscript  refers 
to  a character  string  no  longer  than  40; 
i.e.,  TITLE  has  the  attributes  TITLE  (2) 
CHARACTER  (40)  VARYING.  TITLE  (1)  appears 
at  the  beginning  of  the  plot  alo;ag  with 
constant  identifying  phrases  (described  in 
Fig.  L- 2)  and  at  the  top  of  each  plot. 
TITLE  (2)  appears  along  the  vertical  axis 
of  each  plot.  If  the  value  -999  is  sup- 
plied as  the  last  value  in  any  row  of 
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FORMAT  then  the  values  of  TITLE  will  be 
changed  by  reading  new  values  from  the 
SYSIN  file.  In  particular,  when  plots  are 
being  made,  (FORMAT  (1,1)  * 0)  and  FORMAT 
(n,y)  = -999  then  the  SYSIN  file  will  be 
read  for  two  character  strings,  in  LIST 
format,  each  not  longer  than  40  charac- 
ters. 

The  primary  purpose  of  routine  PLOTTER  is  to  store 
all  of  the  temperature  data  (the  T vector)  on  the  file 
DDNAME  at  each  time  pcint  requested.  If  FORMAT  (1,1)  is 
not  equal  to  zero  then  a plot  will  also  be  made;  for 
this  case,  in  addition  to  the  SHTP  subroutine  library, 
the  library  containing  the  CalComp  software  must  be  made 
available  and  a DD  card  of  the  following  type  must  be 
used : 


//PLOTAPE  DD  UNIT=  (TAPE9 ,, DEFER)  , 

//  VOL=SER=PLOTp,LABEL=  (,BLP)  , 

//  DCB-  (RECFM  = F, LRECL  = 480 , DEN  = 3) 

where  the  "pn  indicates  the  priority  desired  for  the 
off-line  running  of  the  plot;  MoM  is  for  overnight,  "N" 
is  normal,  "Q"  is  guick(i.e.,  twice  the  normal  charge), 
and  "E"  is  emergency  (i.e.,  three  times  the  normal 
charge) . 


PLOTEM 


The  routine  that  actually  performs  the  plotting 
functions  is  named  PLOTEM  and  is  called  by  PLOTTER  on 
the  condition  mentioned  above;  FORMAT  (1,1)  * 0.  When  a 
plot  is  required  after  an  SHTP  run,  then  a small  program 
calling  PLOTEM  must  be  executed.  This  miniprogram  need 
only  have  the  required  variable  declarations  and  a 
single  call  to  PLOTEM,  as  follows  (c.f.  Appendix  0); 

CALL  PLOTEM (DDNAME, FORMAT, THESCAL,TIMSCAL,  TIM- 
PINIT,TIMINIT, COUNT)  ; 


where: 

COUNT  is  the  number  of  time  points  at  which  the 
temperature  data  is  stored  during  the  SHTP 
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run.  This  nunber  is  printed  by  PLOTTER  at 
the  end  of  an  SHTP  run  as  follows: 

"PLOTTER  FILE  •DONABE'  HAS  BEEN  COM- 
PLETED. TEMPERATURES  AT  XXX  TIME 
POINTS  HERE  RECORDED. " , where  XXX  is 
the  value  required  for  COUNT  (FIXED 
BIN) ; and 

DDNABE  is  a CHARACTER  (7)  VARYING  string  that 

identifies  the  DD  card  (file)  that  descri- 
bes the  data  set  on  which  are  stored  the 
temperatures  that  were  saved  from  a 
previous  execution  of  the  SHTP.  In  Appen- 
dix 0 , DDNABE  has  the  value  'PLTFILE*. 

The  other  variables  in  the  calling  sequence  are  as 
described  above  for  PLOTTER. 

When  each  of  the  x plots  have  been  completed,  a 
message  is  printed  by  PLOTEM; 

"PLOT  NUMBER  n HAS  BEEN  COMPLETED." 
indicating  a normal  execution  cf  the  routine. 

Figures  L-1  and  L-2  show  sample  plots  in  which  six 
capacitor  numbers  (105,  109,  113,  117,  121,  and  125)  are 
plotted.  The  first  figure  (Fig.  1-1)  shows  the  iden- 
tifying header  and  the  date  the  data  were  stored.  Also 
TITLE (1)  is  printed.  The  plot  itself  shows  the  time 
axis  label  and  the  user  supplied  label  (TITLE  (2))  ad- 
jacent to  the  vertical  axis  along  with  the  various  scale 
annotations.  At  the  end  of  each  separate  curve  is  the 
node  number  for  the  particular  capacitor.  Figure  13 
shows  the  plot  resulting  from  the  formal  plotting  option 
ABS  (FORM  AT  (X,y))  = 999). 

Thermal  stresses  can  be  plotted  via  the  PLOTTER/ 
PLOTEM  routines  by  setting  the  PLTCODE  argument  in 
routine  SIGMA  (or  SIGNET)  appropriately  and  then  using 
that  sane  numerical  value  as  an  entry  to  the  FORMAT 
array. 

Sttnati-llais 

The  routine  PLOTTER  will  record  all  of  the  thermal 
node  and  psuedo-node  values  on  the  file  named  by  DDNAME. 
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If,  at  the  end  of  an  SHTP  run,  various  nodes  have  been 
selected  by  entries  in  FORMAT  (i.e.  FORMAT  (1,1)  * 0) 
then  plots  will  be  generated  via  an  internal  CALL  to 
PLOTEM.  The  routine  FLOTEM  proceeds  as  follows:  first, 
the  data  saved  on  the  file  naaed  by  DDNAHE  are  read  in; 
for  each  entry  to  TMPSCAL  that  is  negative,  a list  of  y 
normalization  values  are  read  in  LIST  format  from  SYSIN; 
after  all  such  values  are  read,  values  for  TITLE  will  be 
read  from  SYSIN  only  if  the  value  *999  is  input  as  the 
last  value  to  a row  in  FORMAT.  (Note  that  if  the  value 
of  the  last  number  in  the  last  row  of  FORMAT  (i.e., 
FORMAT  (z,y))  is  -999  then  formal  plots  will  be  made  and 
two  new  TITLE  values  will  be  read  for  the  x-th  plot; 
also,  the  values  ±999  occuring  in  FORMAT  will  not  cause 
a curve  to  be  drawn  for  node  999,  rather  the  value  has 
the  same  effect  as  a zero.)  The  reading  of  new  titles, 
if  indicated,  occurs  just  prior  to  the  generation  of  the 
plotting  commands  for  each  plot.  The  sequence  is  then: 

1 All  data  read  frcm  file  identified  by 
DDNAME, 

2 all  nor nalization  factors  read  in  as 
indicated  by  values  in  TMPSCAL, 

3 normalization  factors  applied  to  appro- 
priate values, 

4 new  values  of  TITLE  read  in  as  indicated 
by  the  last  entries  to  rows  in  FORMAT,  and 

5 Plotting  commands  written  out  to  PLOTAPE 
file. 


THE  .JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

Laurel  Maryland 


Appendix  N 
UTILITY  ROUTINES 


Many  of  the  routines  in  the  SHTP  deal  numerically 
with  functions  that  are  represented  in  tabular  form.  As 
a result,  considerable  use  is  wade  of  interpolation 
(table  look-up).  Several  varieties  of  interpolation  have 
been  developed  for  particular  types  of  data  and  are 
embodied  in  so-called  utility  routines.  The  usage  of 
these  routines  should  not  be  of  direct  concern  to  the 
user.  However,  it  could  prove  helpful  if  the  user  were 
aware  of  the  assumptions  and  techniques  employed  as 
errors  often  occur  due  to  inconsistencies  in  the  tabular 
data  used  by  these  routines.  A detailed  discussion  of 
the  interpolation  techniques  is  given  in  "Interpolation 
Techniques,"  Section  1 of  Volume  1. 


DECIDE 


The  routine  DECIDE  is  the  most  frequently  called 
routine  of  all  the  SHTP  subroutines.  It  is  used  to 
provide  values  of  thermal  conductivity  and  heat  capacity 
as  they  vary  with  temperature  for  each  node  in  a given 
problem  as  well  as  many  other  variables.  It  is  a linear 
interpolator  with  a serial  table  look-up  procedure; 
i.e. , it  searches  serially  through  the  table  of  indepen- 
dent variable  values  until  it  finds  two  entries  that 
bracket  the  given  value  of  independent  variable;  the 
returned  value  is  then  linearly  determined.  DECIDE  is  a 
function  subroutine  (i.e.,  it  may  appear  to  the  right  of 
an  equal  sign  in  any  algebraic  expression) ; for  example: 

Y * DECIDE  (KNOHNX , CEP,  INDEP)  ; 

where: 

Y is  the  variable  whose  value  will  be  deter- 

mined by  the  routine  DECIDE  (FLOAT  BIN)  ; 

KNOHNX  is  the  known  value  of  the  independent 
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variable  (FLOAT  EIN); 

INDEP  is  the  tabular  list  of  independent 

variable  values  uhose  first  entry  (INDEP 
(1))  is  the  number  of  dependent  variable 
values  (#N)  ; INDEP  must  therefore  be 
dimensioned  at  least  to  (#N*1)  ((*)  FLOAT 

BIN)  ; 

DEP  is  the  tabular  list  of  dependent  variable 

values  ((*)  FLOAT  BIN)  that  correspond 
with  the  data  values  INDEP (2)  through 
INDEP  (#  N+1 ) . 

If  INDEP  is  dimensioned  only  to  1,  then  the  corres- 
ponding single  value  stored  in  DEP  Hill  be  the  value 
returned.  This  option  is  made  available  to  facilitate 
the  inclusion  of  constants  where  functional  values  are 
normally  required. 


PIF1 


The  PIF1  routine  is  a linear  interpolator  similar 
to  DECIDE  except  it  expects  the  data  in  its  arguments  in 
a slightly  different  format: 

Y = PIF1  (KNOWN!,  INDEP,  N,  DEP)  ; 

where  the  variables  are  as  previously  defined  for  DECIDE 
and  N has  the  same  value  as  #N  previously  used  and  has 
the  FIXED  BIN  attribute.  PIF1  does  not  reguire  the 
number  of  entries  in  the  table  (#N)  to  be  stored  as  the 
first  number  in  INDEP  so  that  the  values  in  INDEP  cor- 
respond 1 to  1 with  the  values  in  DEP  for  this  routine. 


AIT 


Subroutine  AIT  uses  a modification  of  Aitken's 
method  to  interpolate  tabular  data.  The  routine  is 
called  as  follows: 

Y = AIT  (KNONNX , INDEP,  N,  DEP); 
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and  is  used  exactly  as  PIF1.  AIT  is  used  primarily  by 
the  SIGBA/SIGBET  routines  to  interpolate  temperature 
versus  radius  tables  and  produce  smooth  curve  fits. 


PIF1D 


PIF1D  is  a linear  interpolator  that  also  return.*- 
the  derivative  of  the  dependent  variable  with  respect  to 
the  independent  variable.  The  routine  is  used  as  fol- 
lows: 


CALL  PIF1D  (DEP,  1NDEP,  KNOHNX,  Y , DYOX)  ; 

with  DEP  and  INDEP  being  the  dependent  and  independent 
variable  tables  ( (*)  FLOAT  BIN);  KNOHNX  being  the  known 
value  of  the  independent  variable  and  Y being  the  re- 
quired dependent  variable  so  that  Y =£  (KNOHNX)  . DYDX 
will  be  the  slope  of  the  function  at  the  value  KNOHNX; 
i .e.  , DYDX  = f»  (KNOHNX).  Y and  DYDX  are  FLOAT  BIN 
variables.  As  with  DECIDE,  PIFID  requires  the  indepen- 
dent variable  list  to  begin  with  the  number  of  data 
entries  and,  therefore,  be  dimensioned  to  at  least  #N*1. 


LINLOG 


The  LINLOG  routine  is  a linear-logarithmic  inter- 
polator that  takes  the  logarithm  of  the  dependent 
variable  tabular  values,  interpolates  linearly  between 
these  logarithmic  values  and  then  returns  the  antilog  of 
the  interpolated  value.  This  procedure  might  also  be 
represented  as: 

Y = EXP  (PIF 1 (KNOHNX,  INEEP,  N,  LOG  (DEP))); 

The  routine  is  called  as  follows: 

Y = LINLOG  (INDEP,  DEP,  KNOHNX,  N)  ; 

where  the  variables  have  the  definitions  as  previously 
stated.  Note  that  LINLOG  does  not  require  the  number  of 
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data  table  entries  (N)  to  be  stored  as  the  first  word  in 
INDEP. 

The  algorithm  in  LINLOG  makes  use  of  some  assump- 
tions that  the  user  should  be  aware  of.  If  the  inter- 
polation requires  the  use  of  a tabular  value  of  DEP 
between  *1  and  -1  then  a strictly  linear  interpolation 
is  used.  For  DEP  table  values  less  than  -1,  the  loga- 
rithmic interpolation  is  carried  out  using  absolute 
values  and  the  returned  value  will  have  the  correct 
sign.  The  routine  automatically  limits  the  absolute 
values  of  entries  in  DEP  to  less  than  1043.  It  is 
generally  good  practice  to  include  values  in  DEP  and 
INDEP  that  cover  DEP  values  near  +1,  -1,  and  0 if  the 
function  represented  by  DEP  is  both  positive  and  nega- 
tive. 


BIVLID 


BIVLID  is  a linear  interpolator  of  dependent  varia- 
bles that  vary  with  two  independent  variables.  The 
routine,  according  to  a control  value,  will  also  return 
the  partial  derivatives.  (The  name  is  an  acronym  formed 
from:  Bivariate  Linear  Interpolator  with  Derivatives.) 

In  this  case  y = f(X1,X2)  and  the  derivatives  are  either 
with  respect  to  XI  or  X2.  The  routine  can  also  perform 
an  inverse  interpolation;  i.e.,  given  a value  of  the 
dependent  variable  and  a value  of  one  of  the  independent 
variables,  the  corresponding  value  of  the  other  indepen- 
dent variable  will  be  returned.  The  routine  is  called 
as  follows: 

CALL  BIVLID (FOFX S , XI,  X2,  X1VAL,  X2VAL,  Y,  DYDX , 
FLAG)  ; 

with: 

FOFXS  a two-dimensional  array  of  dependent 
variable  values  ((*,*)  FLOAT  BIN), 

XI  a vector  of  independent  variable  values 

((♦)  FLOAT  BIN)  , 

X2  the  other  vector  of  independent  variable 

values  ((*)  FLOAT  BIN)  , 
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X1VAL  a given  value  of  the  independent  variable 
represented  by  XI  (FLOAT  BIN), 

X2VAL  a given  value  of  the  independent  variable 
represented  by  X2  (FLOAT  BIN), 

X the  required  interpolated  value,  (FLOAT 

BIN)  , 

DYDX  the  returned  value  of  the  partial  deriva- 

tive as  required  by  the  value  of  FLAG 
(FLOAT  BIN)  , and 

FLAG  a control  variable  to  indicate  which 

partial  derivative  is  required  and  whether 
the  interpolation  is  inverse  or  not  (FIXED 
BIN).  FLAG  can  have  the  values  ♦!,  -1, 

-2,  or  +2. 

If  FLAG  is  negative  then  inverse 
interpolation  is  performed  and  the 
returned  value  will  be:  X1VAL  if  FLAG 
was  -1,  or  X2VAL  if  FLAG  was  -2. 

If  FLAG  is  positive  then  normal  inter- 
polation is  performed  and  Y is  retur- 
ned. In  each  case,  the  partial 
derivative  returned  is  determined  by 
the  absolute  value  of  FLAG:  If  ABS 
(FLAG)  = 1 then  DYDX  is  the  partial 
with  respect  to  XI  If  ABS (FLAG)  = 2 
then  DYDX  is  the  partial  with  respect 
to  X2. 

It  aust  be  noted  that  the  first  subscript  of  FOFXS  has 
corresponding  entries  in  XI  and  the  second  subscript 
position  of  FOFXS  has  corresponding  values  in  X2;  i.e., 
FOFXS (I, J)  is  the  value  of  the  dependent  variable  at 
XI  (I)  and  X2  (J) . 


BIVLLID 


BIVLLID  is  exactly  the  sane  as  BIVLID  except  that 
the  linear-logarithmic  interpolation  nethod  described 
above  for  LINLOG  is  used  (i.e.,  the  routine  is  a 
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Bivariate  Linear-Logarithmic  Interpolator  with  Deriva- 
tive). All  of  the  argunents  and  options  described  for 
BIVLID  are  the  same  fcr  BIVLLIE;  the  calling  sequence  is 
identical,  only  the  interpolation  algorithm  used  by  the 
routine  is  different. 


TRAP 


The  routines  SIGNET  and  SIGHA  sake  use  of  the 
routine  TRAP  to  perform  integrations  of  elastic  modulus 
and  thermal  expansion  data  as  functions  of  radius.  TRAP 
uses  a trapezoidal  integration  algorithm  and  is  called 
internal  to  SIGH A and  SIGNET  as  follows: 

Y = TRAP  (ORDER,  FOFT)  ; 

where 

ORDER  indicates  the  weighting  factor  in  the 
equation  shown  below  (FIXED  BIN),  and 

FOFT  is  the  temperature  dependent  vananle  ( (*) 

FLOAT  BIN)  . 

If  FOFT  is  a function  of  temperature  (FOFT  * f(T))  and 
the  temperature  (T)  is  a function  of  radius  (T  = g(R)), 
then  the  returned  value  from  TRAP  is: 


b 

y • f Rnf  (g  (B) ) dR 
a 

where  the  limits  of  integration,  a and  b,  and  the  other 
functional  tables  required  are  all  variables  that  are 
internal  to  the  SIGHA/SIGNET  routines  and  do  not  enter 
into  the  calling  sequence. 
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CROSIM 


The  solution  to  the  finite  difference  equations 
used  in  SIGMA/SIGMET  are  obtained  in  subroutine  CROSIM. 
The  routine  is  called,  internal  to  SIGMA/SIGMET  as 

CALL  CROSIM  (MAT,  N)  ; 

where 

MAT  is  a coefficient  matrix  for  the  set  of 

equations  being  solved  and  is  dimensioned 
to  (N*2)  by  (N*  2)  ( (*)  FLOAT  BIN  (53)), 
and 

N is  equal  to  2 tines  the  number  of  regions 

the  cylindrical  wall  segment  is  divided 
into  (FIXED  BIN)  . 

The  solution  vector  is  returned  in  the  (N*1)th  column  of 
the  matrix,  MAT.  The  solutions  obtained  in  the  CkOSlM 
routine  are  done  in  double  precision  to  eliminate  sig- 
nificant round-off  error  that  can  occur  otherwise. 


1 
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Appendix  N 


STEADY-STATE  TEMPERATURE  ROUTINES 


A series  of  subroutines  has  been  developed  that 
allow  the  determination  of  the  steady-state  temperatures 
in  bodies  modeled  by  the  SHTP  subroutine  library.  The 
routines  employ  a Gaussian  elimination  technique  that  is 
iterated  to  achieve  the  nonlinear  solution  to  the 
steady-state  heat  balance  equations  for  each  node  in  the 
network  being  considered.  Details  of  the  solution 
technique  are  given  in  Section  2 of  Volume  1. 


STEADY 


The  routines  involved  with  the  steady-state  solu- 
tions have  been  purposely  designed  to  be  completely 
compatible  with  the  existing  SHTP  structure.  The  user 
who  wants  a steady-state  analysis  instead  of  a transient 
analysis  must  take  three  basic  steps: 

1.  Change  the  TYPRUN  argument  in  Subroutine 

STORE  to  a 4,  (c.f.  Appendix  H)  , 

2.  Replace  the  "CALI  STEP"  statement  to  a 
"CALL  STEADY"  statement  as  shown  below. 

3.  Include  the  appropriate  subroutines  during 
the  LINK-EDIT  step  of  the  job  assembly, 
also  described  below. 

In  the  main  program  used  to  describe  the  particular 
problem  being  analyzed,  the  CALL  STEP  statement  is 
replaced  by  a CALL  STEADY  statement  of  the  following 
syntax: 

CALL  STEADY (f ITERS,  RELAX,  T TOLER,  QTOLER,  PUNCH); 

where: 

r- 

IITERS  is  the  maximum  number  of  iterations  to  be 
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allowed  in  attempting  to  solve  the  set  of 
simultaneous  equations  (FIXED  BIN)  ; 

RELAX  is  a constant  factor  which  is  used  in  the 
correction  term  of  an  alternate  solution 
scheme,  described  later  in  this  appendix 
(the  argument  is  unused  in  the  present 
version  of  STEADY)  (FLOAT  BIN)  ; 

QTOLER,  are  the  convergence  criteria  for  the  net 

TTOLEF  heat  flew  (Btu/h)  and  temperatures  (°R) 
(i.e.r  the  steady  state  is  achieved  if 

| Ti  - Ti'l  < TTOLEH  and 

| Qi  | < QTOLEE, 

with  Ti  - Ti'  being  the  difference  in 
temperature  of  nede  i between  iterations 
and  Qi  being  the  net  heat  flow  to  node  i)  ; 
these  criteria  must  be  satisfied  for  each 
node  in  the  network;  and 

PUNCH  is  a code  for  having  the  values  of  tem- 
perature saved  on  an  output  file  named 
SAVE.  If  PUNCH  is  not  equal  to  1,  the 
temperatures  will  not  be  saved;  if  PUNCH 
is  equal  to  1 then  the  temperatures  will 
be  written  on  the  file  SAVE  in  the  follow- 
ing format: 

LASCAP  * N; 

T (1)  - xxxxx 


T (N)  = xxxxx; 

In  order  to  have  the  temperatures  saved,  PUNCH 
must  egual  1 and  a JCL  card  of  the  following  type  must 
be  included  at  execution  time: 
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//SAVE  DD  DSN=group. userid. dsname, 

//  DISP= ( NEW , CATLG ) , 

//  DCB= (RECFH=FB, LRECL=80 , BLKSIZE=800)  , 

//  SPACE=  (3200,  (2,2)  ,RLSE)  , 

//  UNIT=SAVE 

Note  that  LASCAP  is  the  uame  of  the  program  variable 
into  which  is  stored  the  value  N and  "group. userid, 
dsname"  is  an  arbitrary  name  fcr  the  data  set  that  will 
contain  the  saved  temperatures.  N will  be  the  largest 
node  number  in  the  network  and  will  be  equal  to  the 
# CAPS  argument  to  subroutine  STORE  (c.f.  Appendix  H)  . 

In  the  linkage-editing  step  of  the  job  being  run,  the 
following  input  must  be  supplied: 

//L.SYSLIB  DD  D SN =B EE. F RAZE R , B A SIC02 ,DI SP=SHfi 
//L.SYSIN  DD  * 

INCLUDE  SYSLIB  (STEADYH) 

INCLUDE  SYSLIB  (CCMCONM) 

INCLUDE (CAPCONH) 

/* 

If  direct  calls  to  subroutine  RAD,  CON,  FORCER,  or 
FLUX  are  made  in  the  main  program  then  the  following 
INCLUDE  cards  must  also  be  used  after  the  "//L.SYSIN  DD 
*"  card  above. 

INCLUDE  SYSLIB  (P  ADM) 

INCLUDE  SYSLIB  (CONN) 

INCLUDE  SYSLIB  (FCRCERM) 

INCLUDE  SYSLIB  (FLUXM) 

Note  that  no  subroutine  names  need  to  be  changed  in  any 
of  the  user-supplied  code  that  was  or  could  be  used  for 
a transient  analysis.  It  is  important,  however,  to 
include  the  required  input  cards  to  the  link-editor  at 
job  assembly  so  the  proper  subroutine  modules  are  used. 

The  order  in  which  subroutines  are  called  within 
the  SHTP  main  program  loop  is  important  to  the  steady 
state  analysis.  In  particular  the  routines  FORCER,  FLUX 
or  SAERO  must  be  called  after  all  calls  to  CON,  COHCON, 
or  CAPCON  and  before  the  call  to  STEADY.  If  this  simple 
convention  is  followed,  correct  results  will  be  ob- 
tained. 
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The  STEADY  routine  prints  a message  indicating 
convergence,  when  such  is  the  case,  or  nonconvergence  if 
the  criteria  for  convergence  are  not  net  in  illERS 
attempts.  Refer  to  Appendix  E for  a description  of  this 
message. 

During  the  early  developaent  of  the  steady-state 
routines  another  solution  technique  was  employed  and 
later  found  inferior  in  terms  of  machine  tine  usage  and 
accuracy  to  the  one  described  above.  This  earlier 
technique  used  a Newton-Raphson  iteration  schema  for 
arriving  at  the  steady  state  field. 

The  Newton-Raphson  solution  technique  can  be 
utilized  in  exactly  the  sane  way  as  described  above  for 
the  Gaussian  elimination  technique  with  the  following 
changes: 


1.  The  TYPRON  argument  to  routine  STORE 
should  be  set  to  a value  of  2; 

2.  The  following  input  should  be  provided  to 
the  link-editor  in  place  of  that  showr. 
earlier: 

//L. SYSLIB  DD  DSN=B BE. FRAZER ,B ASIC02,DISP=SHR 
//L.SYSIN  DD  * 

INCLUDE  SYSLIB  (STEADYS) 

INCLUDE  SYSLIB  (COMCONS) 

INCLUDE  SYSLIB  (CAPCONS) 

/* 

If  direct  calls  are  made  in  the  main  program  to  the 
subroutines  RAD  or  CON  then  the  following  cards  should 
also  be  included  after  the  "//I.  SYSIN  DD  *"  card: 

INCLUDE  SYSLIB  (RADS) 

INCLUDE  SYSLIB  (CONS) 

3.  A value  needs  to  be  assigned  to  the 
argument  RELAX  in  the  CALI  statement  to 
STEADY.  This  value  is  used  to  modify  the 
successive  corrections  to  the  nodal  tem- 
peratures and  is  normally  1.  Refer  to 
'•Steady-State  Temperature  Fields"  (Section 
2 of  Volume  1)  for  a complete  description 
of  the  relaxation  factor. 
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I 


No  other  changes  are  required  and  the  output  aessages 
generated  are  the  saae  as  previously  mentioned. 

I 

I 

y 

i 


Appendix  0 

GEN EG  ATI MG  LIMITS  PLOTS 
FFOM  PREVIOUSLY  STORED  DATA 


During  an  URLIM  run  the  pertinent  radome  Imita- 
tions data  are  saved  automatically  on  the  peripheral 
storage  device  specified  by  a data  definition  (DD)  card 
for  the  file  LIMPLT.  These  data  embody  all  of  the 
velocity  limits  data  that  were  determined  during  the 
URLIM  run  and  can  be  accessed  for  display.  The  follow* 
ing  will  shov  how  a PI/I  program  can  be  written  to  read 
the  data  stored  on  the  LIMPLT  file  and  generate  the 
appropriate  limits  plcts.  The  data  stored 'on  the  file 
LIMPLT  were  put  there  by  a call  to  routine  LIMIT  whose 
usage  is  discussed  in  Appendix  L. 

The  data  that  are  reguired  tc  produce  a velocity 
limits  plot  include  the  trajectory  (velocity  versus 
time)/  maximum  material  capabilities,  and  the  appropri- 
ate failure  times.  All  of  these  data  along  with  appro- 
priate control  values  are  read  by  the  program  shown 
below  from  the  LIMPLT  file.  The  program's  main  function 
is  to  re-establish  the  proper  crder  of  storage  for  the 
data  within  addressable  memory  and  then  invoke  the 
routine  that  creates  the  limits  plot,  namely  LIMPRNT . 

The  significant  points  to  notice  about  the  program  shown 
below  and  the  use  of  subroutine  LIMPRNT  are  as  follows: 

1.  For  each  limit  plot  called  for  (iCALLS) 
there  will  be  a limits  structure  (LINSIR) 
allocated  and  read  in  from  the  LIMPLT 
file; 


2.  The  external  E01NTER  variables  LIMPNT  and 
THPNTR  must  have  the  addresses  of  the 
limits  structure  and  the  trajectory  data 
properly  assigned;  and 

3.  The  largest  values  of  time  and  velocity 
from  the  trajectory  tables  (MAXTIHE  and 
NAXV)  must  be  calculated  and  provided  to 
the  LIMPRNT  routine. 


1, 
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Appendix  P 

OBTAINING  CALCCMP  PLOTS 
FROM  PREVIOUSLY  STORED  DATA 


The  PLOTTER  routine  (described  in  Appendix  L)  can 
be  used  to  generate  tine  history  plots  of  program  data 
(temperatures,  stresses,  etc.)  either  at  the  time  of  the 
heat  transfer  run  or  later  by  using  the  data  as  stored 
on  a peripheral  device.  The  following  short  program 
shows  how  the  data  may  be  read  from  the  appropriate  data 
files  and  the  plots  made  subsequent  to  an  SHTP  run.  The 
program  shown  simply  provides  a means  for  conveniently 
reading  in  the  desired  control  values  and  calling  the 
plotting  routine,  PLOTEM.  The  variables  are  the  same  as 
defined  in  Appendix  L.  Note  that  the  program  will 
require: 

1.  Access  to  the  load  module  libraries  con* 
taining  the  routine  PLCTEM  and  the  CalComp 
routines  PLOT,  PLOTS,  LINE,  NUMBER,  SYM- 
BOL, and  WHERE;  (at  APL  this  access  is 
automatic  via  system  convention) 

2.  The  file  PLOT  to  be  defined  and  contain 
the  stored  temperature  data; 

3.  A PLOTAPE  file  to  be  defined  for  the  Cal- 
Comp  generated  data.  These  items  and 
files  are  discussed  fully  in  the  PLOTTER 
discussion  in  Appendix  L. 
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The  program  just  listed  was  executed  with  the  cards 
shown  next  to  generate  the  plot  shown  in  Fig  L-2.  Close 
observation  will  show  the  use  of  the  formal  grid  format, 
normalization  of  plots,  and  title  changing  (c.f.  Appen- 
dix L,  PLOTTER  and  PLCTEM) . 
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Appendix  Q 


External  Variables 


Throughout  the  SNTP,  communication  of  calculated 
data  between  subroutines  is  accomplished  via  the  class 
of  variable  known  as  EXTERNAL.  Since  external  variables 
must  have  unique  names  the  following  table  is  provided 
to  show  the  external  variable  names  now  in  use.  The 
data  type  FLOAT  means  float  binary,  single  precision, 
FIXED  means  fixed  binary  and  is  for  integers.  Asterisks 
for  the  extent  of  an  array  indicate  that  the  variable  is 
controlled,  and  has  its  storage  dynamically  allocated. 
Also  when  variable  names  occur  as  the  extent  of  an 
array,  the  variable  is  controlled  and  the  extent  is  the 
value  of  that  variable.  Further,  these  variables  are  as 
defined  in  Table  1 of  this  report. 


Variable  Nane 

Data 

Value  or 

and  Dimension 

Type 

Use 

ALPHA  (1 ALPHA) 

FLOAT 

Angle  of  attack  (deg) 

AERSIG  (CALLLOD,  2) 

FLOAT 

Aerodynamic  load  stresses 
(psi) 

BSER  (CALLBSE) 

FLOAT 

Boresight  error  slope 

Cl# 

FIXED 

number  of  capacitor  with 
smallest  allowable  time 
step 

CALPHA  (2, 5) 

FLOAT 

array  of  coefficients  used 
in  heat  transfer  correla- 
tions 

CTl  (*) , CT2 (*) 

FLOAT 

lists  of  contact  temperat- 
ures in  composite  conduc- 
tances (°B) 

C(#NODES) 

FLOAT 

heat  capacity  values  for 
each  node  (Btu/ft*-°B) 
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Variable  Name 
and  Dimension 

Data 

Type 

Value  or 

Use 

CNN (# ALPH A, IMACH) 

FLOAT 

normal  force  coefficients 
for  missile 

CNR  (IALPHA,  IMACH) 

FLOAT 

normal  force  coefficient  of 
the  radome 

CAR  (IALPHA,  IMACH) 

FLOAT 

axial  force  coefficient  of 
the  radome 

DCHI (CALLBSE) 

FLOAT 

list  of  electrical  thick- 
ness changes  for  radome 
vails 

DLTHS 

FLOAT 

the  time  step  (s) 

EQ(*) 

FLOAT 

values  of  derivatives  of 
the  net  heating  rate  of  a 
node  with  respect  to  the 
node's  temperature  (Btu/h- 
°R) 

END  TIH 

FLOAT 

time  for  termination  of  the 
run  (h) 

FIRST 

FIXED 

loop  counter 

FILENAME 

RECORD 

INPUT 

FILE 

used  by  PLOTEM  for  reading 
in  data  to  be  plotted 

FLPNTR  (*,*) 

POINTER 

addresses  of  tables  of 
local  flow  ratios  versus 

Mach  number 

GAHHA  (*,*) 

FLOAT 

values  of  specific  heat 
ratio  of  air  versus  tem- 
perature and  pressure 

HAIR  (*,*) 

FLOAT 

values  of  enthalpy 

(Etu/lbm)  for  air  versus 
temperature  and  pressure 


H HAT (INODES, INODES)  FLOAT  conductance  matrix  for 

steady-state  solutions 
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Variable  Name 

Data 

Value  or 

and  Dimension 

Type 

Use 

H 4# NODES) 

FLOAT 

thermal  conductance  values 
fcr  each  node  (Btu/h-°R) 

IA  (*)  t IB(*) 

FIXED 

lists  of  nodes  that  are 
joined  by  composite  conduc- 
tances 

ICOBCN 

FIXED 

number  of  composite  conduc- 
tances 

IASEQ  (*) 

FIXED 

list  ot  indexes 

LASCAP 

FIXED 

number  of  theraal  nodes 

L1MPNT  (CALLLI H) 

POINTER 

list  of  addresses  for 

LIMSlfi  structure 

LASTIH 

FLOAT 

value  of  a previous  time 
point  (s) 

BACH#  (#BACH) 

FLOAT 

Bach  number  values 

BAXPRNT 

PBINT 

FILE 

optional  file  for  extended 
output 

BINPRNT 

FIXED 

Flag  for  print  option 

BAXSTRS  (CALLSIG) 

FLOAT 

Baximum  tensile  thermal 
stress  (psi) 

BAXBSEB  (TUBOPS) 

FLOAT 

maximum  allowable  boresight 
error  values  for  thernal 
materials 

BOB (HCPBOPS) 

FLOAT 

flodulus  of  Rupture  values 
for  mechanical  materials 
(psi) 

BIMBON 

FIXED 

flag  for  limiting  execution 
in  UBLIfl  runs 

BOAIB  (*,*) 

FLOAT 

values  of  viscosity  for  air 

(lba/ft-s)  versus  tempera- 
ture  and  pressure 
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Variable  Name 

Data 

Value  or 

and  Dimension 

Type 

Use 

PALT (*) 

FLOAT 

static  pressure  values 
(function  of  Z,  lbf/ft2) 

PP (45) 

POINTER 

list  of  theriual  property 
table  addresses 

PLTFILE 

REnnrrr- 

CUTPUT 

FILE 

-Jrsenl  fry  PLOTTER  to  store 

data  for  plotting 

PAIR  (*) 

FLOAT 

logarithm  of  pressure 
variable  lor  air  property 
values 

PRAIR  (*,*) 

FLOAT 

values  of  Prandtl  number 
for  air  versus  pressure  aud 
temperature 

PRINTIT 

BIT  (1) 

flag  indicating  output 

Q (♦NODES) 

FLOAT 

S 

beating  rate  values  for 
each  node  (Btu/h) 

QE  (ANODES) 

FLOAT 

values  of  net  heat  flows 
fer  steady-state  solutions 
(Etu/h) 

SAVSIG (CALLLOD, 2) 

FLOAT 

aerodynamic  load  stresses 
at  previous  time  step  (psi) 

SOUND  (*) 

FLOAT 

sonic  velocity  (ft/s) 
function  of  Z 

SYSPRINT 

PRINT 

FILE 

standard  system  output  file 

SYSIN 

STREAK 

INPUT 

FILE 

standard  system  input  file 

SAVBSER  (CALL BSE) 

FLOAT 

boresight  error  slope  at 
previous  time  step 
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Variable  Name 
and  Diaension 

Data 

Type 

Value  or 

Ose 

SAVSTRS  (CALLSIG) 

FLOAT 

aaxiaua  tensile  thermal 
stress  at  previous  time 
point  (psi) 

STOP 

BIT  (1) 

flag  for  stopping  execution 

SAVT (tNODES) 

FLOAT 

ncdal  temperatures  from 
previous  time  point  (°R) 

SAVE 

STREAM 

CUTPU1 

FILE 

output  file  tor  saving 
nodal  temperature  data  from 
steady  state  solutions 

T (#NO  DES) 

FLOAT 

temperature  values  for  each 
node  (° R ) 

TIMS 

FLOAT 

time  variable  (s) 

TMPNTR (CALLTIN, 5) 

FOIHTLF 

address  table  for  time- 
dependent  variables 

TEMPIN 

FLOAT 

initial  temperature  for  all 
nodes  at  start  of  run  (°R) 

TAIR  (*) 

FLOAT 

values  of  temperature  for 
air  property  tables  (°R) 

TALT  (*) 

FLOAT 

values  of  temperature 
versus  altitude  (®e) 

TSP  ACE  (*) 

FLOAT 

values  of  the  space  tem- 
perature versus  altitude 
<°R) 

TIM 

FLOAT 

current  simulation  time  (h) 

VI(*) 

FLOAT 

a list  of  nodal  volumes 
(ft*) 

XCPL (# ALPHA, iMACH) 

FLOAT 

center  of  pressure  values 
for  radome 

XCAP (XCAPLIM , 7) 

FLOAT 

nodal  interconnection 
matrix 
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Variable  Name 

Data 

Value  or 

and  Dimension 

Type 

Use 

z<*) 

HO  AT 

altitude 

(ft) 

ZAIR  (*,*) 

FLOAT 

values  of 

compressibility 

for  air  versus  pressure  and 
temperature 
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